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The Aerodynamics of Golf Balls 


Joun M. Davies 
Akron, Ohio 


(Received May 23, 1949) 


The aerodynamic forces on golf balls were studied by dropping 
spinning balls through the horizontal wind stream of the B. F. 
Goodrich wind tunnel. The lift, Z, and drag, D, were calculated 
from the drift of the balls, rotating at speeds, V, up to 8000 r.p.m. 
while falling through a wind stream having a velocity of 105 feet 
per second. For a standard dimple or mesh golf ball weighing 
0.1 Ib., the lift varied with the rotational speed as L=0.064 
X[1—exp(—0.00026N)], with a maximum observed value of 
0.055 Ib. or more than half the weight of the ball. The drag in- 
creased nearly linearly from about 0.06 lb. for no spin to about 0.1 
Ib. at 8000 r.p.m. 


For a smooth ball the lift was negative at all rotational speeds 
below 5000 r.p.m. Above this speed, the lift was positive but was 
less than for the standard ball. The drag for these balls was nearly 
constant at about 0.08 lb. Balls with shallower dimples than 
standard gave intermediate results. Driving tests were consistent 
with the wind tunnel results. 

These results explain why a golfer cannot obtain long drives 
with a ball having a smooth surface and why the standard dimple 
or mesh surface gives him greater distance and better control 
of the ball. 





INTRODUCTION 


HE force acting on a sphere or cylinder rotating 

in a wind stream about an axis perpendicular to 

the stream is generally not in the direction of the wind. 
It may be resolved into two components, one in the 
direction of the wind and one perpendicular to the 
wind and to the axis of rotation. Magnus! and Newton! 
discussed the transverse force and gave the accepted 
explanation for its existence many years ago but ap- 
parently only a limited number of experimental studies 
have been made. Several measurements*-* have been 


' (a) Magnus, Poggendorf’s Ann. d. Physik Chemie 88, 1 (1853) ; 
(b) Newton, Phil. Trans. 6, 3078 (1672); (c) Rayleigh, Messenger 
of Math. 7, 14 (1877); (d) Thomson, Nature 85, 251 (1910), 
reviewed in Sci. Am. p. 136 (February 11, 1911); (e) Goldstein, 
Modern Developments in Fluid Dynamics (Oxford University Press, 
New York, 1938), Vol. 1, p. 83; (f) Prandtl, Tietjens, and den 
Hartog, Applied Hydro and Aeromechanics (McGraw-Hill Book 
Company, Inc., New York, 1934), first edition, pp. 82-85. 

? (a) Lafay, Rev. Mécanique 30, 417 (1912); (b) Int. Crit. Tab. 
1, 408 (1926). 

3(a) Reid, N.A.C.A. Tech. Note 209 (1924); (b) see reference 
2(b); (c) see reference 1(e), Vol. 2, p. 545. 

*(a) Betz, Zeits. des Vereines deutscher Ingenieure 69, 11 
(1925); (b) see reference 1(e), Vol. 2, p. 545. 

* Thom, Brit. Aero. Res. Comm. Rep. and Memo. No. 1018, 
Annual Reports 1925-26, p. 82 (1927). 

* Acheret, Zeits. f. Fiugtech. Motorluftschiffahrt 16, 45 (1925). 

? Flettner, Zeits. f. Flugtech. Motorluftschiffahrt 16, 52 (1928). 

* Thom, Brit. Aero. Res. Comm. Rep. and Memo. No. 1082, 
Annual Reports 1926-27, p. 66 (1928). 


made of the forces acting on rotating cylinders in air 
and photographs have been taken of the flow around 
cylinders in air® and in water.*!° The Flettner rotor 
ship’ was an attempt to make a practical application of 
these forces. Robins" observed the effect of such forces 
on bullets fired from gun barrels curved in such a way 
as to produce rotation of the bullet. Maccoll? measured 
the forces on a 6-inch diameter sphere rotating at 
speeds up to about 1800 r.p.m. in a wind stream having 
velocities up to 34 feet per second. 

These forces are important in many ball games, for 
example, baseball, tennis, and golf, but little definite 
information seems to be available. Measurements on 
non-rotating golf balls have been made at the Bureau 
of Standards® and by Derieux,“ but in both cases 
relatively low velocities were used. Thomson! measured 
the difference in pressure between opposite sides of a 
golf ball rotating in a wind stream and showed qualita- 
tively how the forces would affect the flight of the ball. 


® See reference 1(f) pp. 281-287. 

10 Relf and Lavender, Brit. Aero. Res. Comm. Rep. and Memo. 
No. 1009, Annual Reports 1925-26, p. 74 (1927). 

1 Robins, “New principles of gunnery,” p. 206 (1805) ; presented 
before the Royal Society in 1747, discussed in reference 15(d). 

12 (a) Maccoll, J. Roy. Aeronautical Soc. 32, 777 (1928); (b) see 
reference 1(e), Vol. 2, p. 504. 

48 Unpublished results by Dr. H. L. Dryden, now at N.A.C.A. 

4 Derieux, J. Elisha Mitchell Soc. 51, 207 (1935). 
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Tait'® calculated the forces on a golf ball in flight by 
observing the trajectory and the time of flight and then 
attempted to explain the trajectory in terms of these 
forces. However, he had no convenient way of determin- 
ing the ball velocities and, in general, the values he 
assumed were too high and the air forces correspond- 








e 


Vertical Distance 











Horizonte! Distance 
Fic. 1. Forces on a golf ball in flight. 


ingly high but he was able to show that the transverse 
force had a significant effect on the trajectory. 

A typical golf ball trajectory is shown in Fig. 1, 
with a convenient way of defining the aerodynamic 
forces. The total force, R, is resolved into two com- 
ponents, a drag, D, in the line of and opposite to the 
motion, and a lift, L, perpendicular to the motion. The 
lift is upward for the usual case of under spin and down- 
ward for over spin. Frequently the axis of rotation is 
not exactly horizontal and then these forces have hori- 
zontal components which cause “hooking” and “slicing.”’ 
In play velocities vary over a wide range; translational 
velocities of more than 225 feet per second and rota- 
tional velocities up to 8000 r.p.m. have been observed 
at the A. G. Spalding and Brothers laboratory’ with 
Edgerton’s high speed photographic technique." To 
learn more about the lift and drag under flight condi- 


tions, the motions of different types of golf balls spinning 
in a wind stream were studied in the B. F. Goodrich 
wind tunnel. 


METHOD OF MEASURING LIFT AND DRAG 


The spinning device is shown in Figs. 2 and 3. The 
ball was held between shallow cups and rotated about a 
horizontal axis perpendicular to the wind stream with a 
variable speed motor whose speed was measured with a 
tachometer. When the trigger was released, the springs 
pulled the cups apart, allowing the ball to fall freely 
through the wind stream. 

The wind tunnel'* is of the horizontal, open type, 
with the throat about 18 inches high and 36 inches wide. 
It was modified somewhat for these tests by adding the 
top and bottom sections, A and B, and a screen, as 
indicated in Fig. 4. The ball was rotated at C and the 
landing spot was marked on waxed paper at B. In some 
cases the ball struck the screen before falling to B and 
for those tests a raised platform, § inch thick, was used 
at E. The various heights of fall were from 0.67 to 1.3 
feet. The wind speed was 105 feet per second. 

The motion of the ball under the action of gravity 
and the aerodynamic forces is indicated in Fig. 5. 
With no spin, the ball is pulled down by its weight and 
horizontally by the drag, so that it travels a horizontal 
distance x» while dropping a vertical distance y. With 
the ball spinning, the value of the drag is different and 
there is an additional vertical force, ordinarily up for 
clockwise rotation and down for counterclockwise rota- 
tion. This increases or decreases the rate of fall, allow- 
ing the ball to travel a horizontal distance x; or x2 while 
falling the vertical distance y in a time /; or ¢2. Neglect- 
ing any variation in lift and drag arising from the in- 
crease in ball velocity during fall and the variation in 
wind velocity throughout the testing space, the equa- 





Fic. 2. The device for rotating and dropping a ball in the wind tunnel. 


6 Tait, (a) Nature 42, 420 (1890); (b) ibid. 44, 497 (1891); (c) ibid. 48, 202 (1893); and (d) Trans. Roy. Soc. Edinburgh 37, 


427 (1893). 


6 (a) Ind. Eng. Chem., News Edition, 16, 673 (1893); (b) other results by Mr. George H. Temple, not published. 
17 Edgerton and Germeshausen, American Golfer, 17 (November 1933). 


18 Geer, Aero. Eng. 4, 33 (1932). 
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tions of motion are: 


Clockwise spin 


Counterclockwise spin 


(W/g)(@x/d?) =D, (1) (W/g)(@x/d?) =D, (5) 
(W/g)(@y/d?)=—-W+L, (2) (W/g)(@y/d?)=—W—-L, (6) 
xi= Dgt;?/2W, (3) xo= Degt.?/2W, (7) 

y= —(W—L)gi?/2W, (4) y= —(W+L)gt.*/2W, (8) 

D= —2Wxx2/y(x1+%2), (9) 


L= W (x1 —2%2)/(a1+22). 


By spinning the ball first in one direction and then the 
other, values of the drag and lift can be obtained. 


RESULTS 


Some typical results are given in Fig. 6 for standard 
golf balls with dimple and mesh markings and for 
smooth balls. Each result is the average of three tests 
of each of two dimple or mesh balls or of four smooth 
balls. The dimples were spherical segments about 0.013 
inch deep and inch in diameter and the mesh surface 
consisted mostly of approximate squares about 0.011 


inch deep and } inch wide. The smooth balls were 


standard with the exception that the covers were ap- 
plied in a highly polished smooth mold. All the balls 
were very close to the standard diameter, 1.68 inches, 
and weight, 1.62 ounces. 

The values for the lift and drag for the standard balls 
show the effects of dimple and mesh markings to be very 
similar. This might be expected as a result of experience 
in play. The values for the lift fall close to the curve 
represented by the equation 


L=0.064[1—exp(—0.00026N)], (11) 


where J is the rotational speed in r.p.m. and L is the 
lift in pounds. At the higher rotational speeds, the lift 
was just about half the weight of the ball. The drag 
increased nearly linearly from about 0.06 lb. at zero 
spin to just about the weight of the ball at 8000 r.p.m. 

The most surprising result is the lift curve for the 
smooth ball. Up to about 5000 r.p.m., or an equatorial 
surface speed of 37 feet per second, the lift is negative, 
just opposite in direction to that expected from the 
Magnus effect. Above that speed it assumes the normal 
direction. The lift is much smaller than for the standard 
balls but there is an indication that at still higher rota- 
tional speeds it will be at least nearly as great. The drag 
of the smooth ball changes very little with change in 
rotational speed ; above 4000 r.p.m. it is less than that 
of the standard balls. 

In Fig. 7 smooth and standard balls are compared 
with a ball with intermediate marking, i.e., shallower 
dimples with rounded edges. Each result is the average 
of five tests on one ball. The results on the smooth ball 
agree roughly with those of the previous test but are not 
as consistent. For the standard dimple ball, the lift 
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(10) 





curve is nearly identical with the previous curve, except 
possibly at 1000 r.p.m., but the drag curve is somewhat 
lower at the higher speeds. These standard balls were 
similar to those of Fig. 6 but the markings were slightly 
different. For the ball with the shallower dimples, the 
lift was intermediate between that of the standard and 
smooth balls, with the drag somewhat lower than for 
either. It is interesting that for a bramble ball, with 
reversed or protruding dimples, the lift and drag curves 


are nearly identical with those for the normal or recessed 
dimples. 





Fic. 3. Photograph of a falling ball just after release. 


As a check on the validity of the method, photo- 
graphs were obtained of the ball paths. The camera 
shutter was set for an exposure of one-fifth to one-half 
second and was opened by a solenoid tripped by the 
dropping mechanism. Figures 8, 9, and 10 show com- 
parisons of the paths for standard and smooth balls. 
The paths are essentially straight and the directions are 
in agreement with the other tests. 


DRIVING TESTS 


Driving tests on smooth balls are consistent with the 
wind tunnel results. The B. F. Goodrich driving ma- 
chine consists essentially of a mallet-type club, rotated 
about a horizontal axis by a motor. The club speed and 
face angle can be varied. The results plotted in Fig. 11 
were obtained with the smooth and standard balls for 
which the wind tunnel results are given in Fig. 7. 
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Each result is the average of six dimple balls or of three 
smooth balls. The distance is the carry or the hori- 
zontal distance from the tee to the spot where the ball 
struck the ground. The roll was not measured. The 
measurement of trajectory height was rather arbitrary, 
the figures representing the apparent maximum angle 
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Fic. 4. Diagram of the wind tunnel, showing the 
location of equipment. 


of the ball observed from a point about ten feet directly 
above the tee. 

For the standard balls, the carry and trajectory 
height increased with the club speed and varied with 
the club face angle in the usual manner. With the 
smooth balls, however, as the speed and face angle 
increased, the balls started out at a higher angle but 
curved downward more rapidly. Consequently, neither 
the carry nor the trajectory height changed appreciably 
with variation in club speed or face angle. The magni- 
tude of the negative lift seemed to increase rapidly 
with linear speed of the ball. 

Driving tests at the A. G. Spalding and Brothers 
laboratory on similar balls gave results consistent with 
these tests.'® 


DISCUSSION OF ERRORS 


The free fall method of studying lift and drag was 
selected because it eliminated the uncertain corrections 
for the effect of the ball supports. However, the results 
were subject to some other errors. 

The wind stream had some turbulence and was 
affected by the equipment in the tunnel. Pitot tube 
measurements showed that the variation in wind ve- 
locity throughout the testing region was, at the most, 
about +4 feet per second. 

Near the end of the path, the ball velocity was ap- 
preciable and this changed the magnitude and direction 
of the drag and lift. Corrections for this effect are small 


- because of the short time in this part of the path. 


Calculations show that for counterclockwise spin the 
errors tend to cancel and the path is practically straight. 

Also, in its initial position, the ball is close to the 
top of the tunnel and the wind velocity between the 
ball and the top of the tunnel is increased, with a 
resulting vertical force on the ball. 


1% Unpublished results by Dr. J. B. Dickson at the A. G. 
Spalding Brothers, Inc., Laboratory. 
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The dropping mechanism had a slight effect on the 
direction of fall in some cases but this was not serious. 
In tests with dimple balls at 5000 r.p.m. with no wind, 
the average path deviated from the vertical by no 
more than +0.007 feet with individual variations of no 
more than +0.014 feet, for a drop of 1.33 feet. 

A combination of these effects, especially the first 
three, causes some curvature of the ball paths, as shown 
in the photographs. However, the photographs and 
also tests made with the landing platform at various 
heights show that these effects are not large. The extent 
of the over-all error is probably less than 10 percent. 


COMPARISON WITH OTHER TESTS 


For the smooth ball with no spin, the Reynolds’ 
number is about 9X10‘ and the drag coefficient, Cp, 
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Fic. 5. Motion of a ball in the wind tunnel. 


about 0.44, where Cp is defined by the equation 
Cp=D/(p/2AV"), (12) 


where V is the velocity, p the air density, A the pro- 
jected area, and D the drag. Tests in wind tunnels 
having low turbulence give slightly higher values for 
Cp, for example, Weiselberger®® at Géttingen found 
Cp=0.48 at this Reynolds’ number. In a more turbu- 
lent air stteam, the critical region where the drag 
coefficient decreases rapidly is reached at lower air 
velocities and this is a reasonable explanation for the 
above difference. 

Apparently the only quantitative data on rotating 
spheres are those of Maccoll" at the National Physical 
Laboratory. Over the same range of rotational speeds, 
his drag and lift curves resemble those in Fig. 6, 
especially in showing a negative lift at low rotational 
speeds. The lift becomes normal in direction at equa- 
torial speeds of about half the wind speed or less. 

Negative lift has been observed also for smooth 
cylinders. Lafay* found a negative lift at surface speeds 
of less than half the wind speed and Reid* and Thom! 


» See reference 1(f), p. 100. 
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indicated a slight negative lift at a slightly lower speed. 
In testing models to study the ground handling forces 


_on airships, Knoblock and Troller! at the Guggenheim 


Airship Institute and Silverstein and Gulick” at 
Langley Field found that for certain conditions, the 
force on a model near a flat surface was negative. 

In Derieux’s experiment, a non-rotating gelf ball was 
suspended by a thread in a horizontal wind stream at 
velocities of 24 to 45 feet per second and the drag 
determined from the angle of deflection. His data at low 
velocities and the present results at 105 feet per second 
fit the equation 

D=23X10-*V!-, (13) 


This agreement may not be significant because un- 
published results from the Bureau of Standards” show 
that the critical region for golf balls is reached at 
velocities of 75 to 90 feet per second and extrapolation 
through the critical region is not justified. It is difficult 
to compare the present results with the measurements at 
the Bureau of Standards. A velocity of 105 feet is in 
the critical region and in this region the drag seems to 
depend very much on the surface of the ball. There is 
a small difference in size, the ball diameter being 1.68 


leer 
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Fic. 6. Drag and lift for smooth and standard dimple and mesh 
balls. Wind velocity about 105 feet per second. 


inches in the present case and 1.62 inches for their 
experiments. Interpolating their results at 105 feet per 
second and making allowance for this difference in size, 
their values of D for two different balls were 0.063 and 
0.050 Ib., one a little higher and the other a little lower 
than the present value of 0.06 Ib. 

It is also difficult to make a comparison with Tait’s’® 
results. He attempted to measure the initial ball speed 
and the trajectory. By assuming a particular depend- 
ence of D and L on translational speed, D« V? and 
L« V, and of L on rotational speed, L« N, he was able 
to calculate both forces. He was uncertain about both 
the translational and rotational speeds and for the most 
part gave his results in terms of assumed values of V. 


* Knoblock and Troller, Guggenheim Airship Inst. Pub. No. 2, 
p. 53 (1935). 


* Silverstein and Gulick, N.A.C.A. Rep. No. 556 (1936). 
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Generally these assumed values were too high and the 
values of D and L correspondingly high. Calculated for 
a velocity of 105 feet per second, his values of D are 
between 0.10 and 0.19 Ib. and the values of L between 
0.06 and 0.12 Ib. Considering the possible differences in 
size and surface condition of the balls, the agreement 
is reasonably good and undoubtedly better agreement 
would be found by substituting accepted values for 
the velocities. 


DISCUSSION OF NEGATIVE LIFT 


Bernoulli’s theorem states that for ideal frictionless 
flow, along any stream line the total energy is constant, 
i.€., 

p+ pV?/2+ pgh=a constant. (14) 


A qualitative explanation of the Magnus effect can be 
given on this basis. With the sphere rotating clockwise 
as in Fig. 12, the velocity at a on stream line A is higher 
than at 6 on stream line B. Since the kinetic energy is 
higher at a than at 5, the pressure is lower at @ than at b 
and there is a resultant upward force. The wake is 
deflected downward. This accounts for the experience 
with a rough sphere such as a standard golf ball and 
with a smooth sphere at high rotational speeds but not 
for the opposite effect with a smooth sphere at low 
rotational speeds. 

At low velocities this situation is disturbed by fric- 
tion and turbulence. With a non-rotating spHere in a 
wind stream, at low velocities the flow follows the sur- 
face for a short distance and then breaks away. As the 
velocity increases, the laminar boundary layer changes 
gradually into a turbulent layer and since such a 
turbulent layer is able to follow the surface better than 
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Fic. 7. Drag and lift for balls with reduced dimple depth com- 


pared with smooth and standard dimple balls. Wind velocity about 
105 feet per second. 


a laminar one, the breakaway point is moved in the 
direction of the wind stream. Relatively, the inertia 
forces increase faster than the frictional forces. Condi- 
tions approach those of ideal flow and there is a decrease 
in the drag coefficient during this transition. 

These effects are illustrated by the pressure distri- 
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bution around a sphere, for example, the measurements 
by Fage* on a 6-inch diameter sphere at various ve- 
locities, as shown in Fig. 13. The critical range was 
from velocities of 50 to 95 feet per second, the curve 
for Vo=35 feet per second corresponding to nearly 
laminar flow and that for Vo=135 feet per second to 








Fic. 8. Photograph of ball paths in the wind tunnel. Standard 
dimple ball, 8000 r.p.m. A. Clockwise. B. Counterclockwise. 








Fic. 9. Photograph of ball paths in the wind tunnel, 8000 r.p.m. 
clockwise. A. Dimple ball. B. Smooth ball. 


nearly completely turbulent flow. Another curve has 
been added, calculated for the case of frictionless flow. 
This is based on the Bernoulli law and the assumption 
that the velocity at the surface is given by 


V=$Vosing, (15) 


where V> is the undisturbed velocity and ¢ the angle 
measured from the stagnation point.™ Below the critical 
region the pressures deviate considerably from the ideal 


% (a) Fage, Brit. Aero. Res. Comm. Rep. and Memo. No. 1766 
(1937); (b) see reference 1(e), Vol. 2, p. 497; (c) similar curves 
by Eiffel in Int. Crit. Tab. 1, 407 (1926). 

*“Milne-Thomson, Theoretical Hydrodynamics (MacMillan 
Company, Ltd., London, 1928), p. 417. 
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curve but above the critical region the experimental 
curves approach the ideal case rather closely. 

To a reasonable approximation, the conditions on the 
opposite sides of a rotating sphere may be represented 
by the pressure curves for a stationary sphere, using 
the appropriate velocities for each side. On the side 
rotating against the wind stream, the effective Rey- 
nolds’ number is high and on the side rotating with the 
wind stream it is low. For low rotational speeds, the 
conditions at @ in Fig. 12 might correspond roughly to 
the curve for Vo=35 feet per second in Fig. 13 and the 
conditions at 5 to the curve for Vo= 135 feet per second. 
This analogy does not hold for the whole surface but 
fits fairly well in the region near ¢=90° which is of 
most importance in determining the lift since the effec- 








Fic. 10. Photograph of ball paths in the wind tunnel, smooth 
ball, 1000 r.p.m. A. Clockwise. B. Counterclockwise—rest of 
path cut off by short exposure. 


tive vertical component of the normal pressure is high 
in this region. On this basis it is evident that the pres- 
sure at b is lower than at a and the transverse force will 
be downward or the lift is negative. 

Some care is needed in arriving at this conclusion as 
the ordinate in Fig. 13 is the dimensionless quantity 
(p— po)/$pV.* and the absolute pressure depends on the 
velocity. The velocity at 5 will be somewhat lower than 
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Fic. 11. Driving distances (Carry) and trajectory heights for 
dimple and smooth balls. 
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at a. The lift will be negative if the equivalent velocity 
for the lower side is less than 0.6 that for the upper side, 
positive if it is more than that. 

Obviously conditions are rather critical and it would 
not take a very great change in the state of turbulence 
or velocities to reverse the direction of the resultant 
force. Presumably this is just what happens when the 
rotational surface speed becomes about half the wind 
speed. As the surface speed becomes higher, there may 
be a tendency to approach the ideal curve on both 
sides of the spinning sphere, but for different reasons. 
On the side rotating against the wind, the effective 
Reynolds’ number becomes high; on the side rotating 
with the wind, with the ball surface moving faster than 
the wind stream, the flow follows the sphere very well 
and the pressures may approach the ideal curve even 
with a laminar boundary layer. In such a case, the 
direction of the resultant force would be determined 
by the usual Magnus explanation. . 

This does not explain the nature of the flow and the 
resulting pressure distribution but accepting that dis- 
tribution, it does account for both the positive and 

















Fic. 12. Flow around a rotating sphere. 


negative lift for a smooth sphere. Undoubtedly the 
situation would be considerably clearer if more pressure 
distribution data were available, for example, for a non- 
rotating rough sphere and for a rotating smooth sphere. 


CALCULATED TRAJECTORIES 


Trajectories for given typical conditions would be 
interesting but without additional information the 
calculations cannot be made with any certainty. 
Neither the manner in which LZ and D vary with trans- 
lational speed nor the time rate of change of rotational 
speed are known. Consequently, only a few calculations 
have been made, for the conditions assumed in Table I, 
and the results are shown in Fig. 14. 


The equations of motion are: 
dx/d?=(gD/W) cos@—(gL/W) sind (16) 
d*y/d? = —g—(gD/W) sind+(gL/W) cos@, (17) 


where the various terms have the meaning indicated in 
Fig. 1. The values of L, D and 6@[tan@= (dy/dt)/(dx/dt) ] 
for the various cases were substituted in these equations 
and the values of x and y obtained by Moulton’s method 
of successive approximations.® 


*® Moulton, New Methods in Exterior Ballistics (University of 
Chicago Press, Chicago, 1926). 
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Fic. 13. Pressure distribution around a sphere in a wind streaam— 
experimental] curves by Fage. 


Data for a direct check on these trajectories are not 
available but the results are reasonable. For all the cases, 
air resistance reduces the maximum theoretical carry 
by at least one-third. For an initial angle of 20°, the 
lift increases the carry appreciably even with increased 
drag. The effect of loss of spin is surprisingly small. 
The carry for the smooth balls seems to be too high, 
very likely because the magnitude of the negative lift 
increases more rapidly with linear speed than was 
assumed. 

The forces acting on a ball at the beginning, near the 
middle and at the end of the flight, calculated for Case 6, 
are given in Table II. It is apparent that after it leaves 
the club, a ball is acted upon by relatively large forces. 
In this case, the aerodynamic forces are from one-third 
to more than twice the weight of the ball and with 
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Fic. 14. Calculated golf ball trajectories. 
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TaBLeE I. Conditions for calculated trajectories. 


Initial velocity, Vo= 200 feet per second. Initial angle with horizontal @)= 20°. 




















No. Description Lift Drag 
1 Ideal case L=0 D=0 
2 Dimple ball, no spin L=0 D=5.2X10-*V? 
=0.057 lb. at 105 feet per second 
3 Smooth ball, no spin or hypo- L=0 D=7.5X10-*V? 
thetical case of dimple ball, no =0.083 Ib. at 105 feet per second 
spin, but same drag as ball spin- 
ning 4000 r.p.m. 
4 Dimple ball, spinning at con- L=4xX10~V D=7.5X10-*V? 
stant velocity of 4000 r.p.m. = (0.042 Ib. at 105 feet per second =0.083 Ib. at 105 feet per second 
5 Dimple ball, spinning. Start at L=6.2X10"V D=9.5X10-5V? 
r.p.m., decrease at rate of < {1—exp[—2.6x 10-* 4000 exp(—0.1/) ]} -X [0.057+-6.5 X 4000 exp(—0.1#)] 
10 percent per second. =(0.042 lb. at 105 feet per second and =0.083 lb. at 105 feet per second 
4000 r.p.m. 
6 Dimple ball. Similar to Case 4 L=4X10“V D=3.54X 10-*V5/8 
except Da V5/3 = (0.042 lb. at 105 feet per second =0.083 Ib. at 105 feet per second 
7 Smooth ball, spinning at con- L=—0.75X10-V D=7.5X10*V? 


stant velocity, 4000 r.p.m. = —0.0079 lb. at 105 feet per second =0.083 lb. at 105 feet per second 








higher translational and rotational speeds they would 
be still higher. Sufficient data are not available to 
calculate accurately the effects of a transverse force 
but such a force of 0.04 Ib., resulting, for example, from 
rotation about an axis at an angle to the horizontal, 
would cause a hook or slice of roughly 50 yards in a 
200 yard drive. A following, head, or transverse wind 
changes the trajectory appreciably. The magnitude of 
these lift and drag forces illustrates the extreme pre- 


TABLE II. Lift and drag for a calculated trajectory. 











Time Distance from 
in Velocity, ft./sec. origin, ft. Lift Drag 
seconds y x v y x Ib. Ib. 
0 68 188 200 0 0 0.080 0.242 
2.6 0 83 83 82 320 0.033 0.056 
5.5 —53 55 83 0 512 0.033 0.056 
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cision necessary for satisfactory control of a golf ball in 
play and, at the same time, the degree of control that 
can be exerted by an expert player with precise swinging. 
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The Projection of Crystal Structure Factors on the Reciprocal Lattice 





Dan McLACHLAN, Jr. 
Department of Metallurgy, University of Utah, Salt Lake City, Utah 
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C. L. Curist 
Stamford Research Laboratories, American Cyanamid Company, Stamford, Connecticut 


(Received August 3, 1948) 


A discussion of the use of the wave addition method for making projections of the crystal structure 
factor, Fixo, upon the reciprocal lattice, in a fashion similar to that used in making electron density pro- 
jections in real space, through the use of the “sand machine” and devices operating on similar principles, 


has been developed. 





INTRODUCTION 


| be recent years, several devices more or less specifi- 
cally designed for the computation of electron 
densities in crystal structure analysis have been re- 
ported in the literature. These machines, all based 
fundamentally upon the idea of wave addition originally 
proposed by W. L. Bragg,' include the photographic 
method* * the sand deposition method,‘ and the oscillo- 
graphic method.® With these devices the function which 
is computed is of the type! 


p(x, y)= La Lee | Faro] cosl2a(hx+ky)—anro]. (1) 


The Fourier equation used for the computation of 
the crystal structure factor is quite analogous to that 
used for the computation of electron density, and has, 
for the centro-symmetric cell containing M atoms, the 
form 


M 


Fiko= >, fi cos2a(hx;+ky;). (2) 


i=1 


By virtue of this analogy, it would seem to be of interest 
to inquire into the feasibility of using the machines 
referred to abovef for the purpose of calculating struc- 
ture factors, rather than electron densities. This paper 
presents an analysis of this problem. 


ELECTRON DENSITY COMPUTATIONS 


A rather general description of the nature of the 
wave addition involved in such devices as the “sand 
machine” has been given.* We may write, for example, 
for the projected electron density p(x, y)’ 


i) @ 
p(x, y)’= DL (Sa)ato, (3) 
h=—o k=—00 

'W. L. Bragg, Zeits. f. Krist. A70, 475 (1929). 

*M. L. Huggins, J. Am. Chem. Soc. 63, 66 (1941). 

*M. L. Huggins, J. Chem. Phys. 12, 520 (1944). 

*D. McLachlan, Jr. and E. F. Champaygne, J. App. Phys. 17 
1006 (1946). 

*R. Pepinsky, J. App. Phys. 18, 601 (1947). 

t Although, by means of certain modifications, the discussion 
which follows may be made to apply to some of the numbered 
strip methods, these will not be given consideration. 
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where 
(Szy)ato= | Frio | cos[ 24(hx+ ky) oe Oinko |. (4) 


The wave (S.y)rzo is shown in projection in Fig. 1. Its 
wave-length Axxo is precisely the interplanar spacing 
dixo, Which may be computed by Bragg’s Law, or from 
the unit cell constants of the crystal under considera- 
tion. Using the latter, dy. has, for the general case, 
the form 

sin8 


diko= ‘ 
[(h/a)?+ (k/b)?— (2hk cosB/ab) }} 


The direction of the wave front may be referred to the 
angle Waxo, which may be calculated from the equation 


Waro= arc sin(hdpxo/ a). (6) 


The quantities dro and Wrro may be readily determined 
graphically by a construction such as is given in Fig. 1. 
The constants dio and Wako, the experimentally deter- 
mined |Fixo|, and the phase angles arko0, obtained by 
other means, then completely define (Szy)rto- 





(5) 


RECIPROCAL LATTICE PROJECTIONS FROM 
CENTRO-SYMMETRIC STRUCTURES 


For structures having a center of symmetry at the 
origin, and including those structures having a center 
of symmetry in projection, we may write _ 


M/2 


Fixo= 2 > f; cos2a(hx;+ky;), . (7) 
i=1 


taking account explicitly of those atoms related by the 
center. The atomic coordinates x; and y; are now to be 
considered as constants fixed by the structure for which 
the Fro are to be computed, while the Miller indices 
h and k are the variables. Similarly to Eqs. (3) and (4), 
we may write 


M/2 
Frro= 2 2 (Fnko) is (8) 
and 
(Faro) i= fi cos2x(hx:;+ ky;). (9) 


At this point, however, the analogy between Eqs. (4) 
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and (9) is not complete. While | Fy.0; is a fixed constant 
independent of the real plotting area defined by x and y, 
fi is not similarly independent of the reciprocal plotting 
area defined by h and k. This difficulty may be resolved, 
if the assumption is made that the atomic scattering 
factor curves, for the particular set of atoms involved, 
all have the same shape as that of some conveniently 
chosen reference atom. Choosing some reference atom r 
whose atomic scattering curve constitutes the best 
average of those involved in the computation at hand, 
we may write 


fAhko) = Z,g(hko). - (10) 


Here, the left-hand member is the atomic scattering 
factor written explicitly as a function of hko, the func- 
tion g(hko) defines the shape of the scattering curve, 
and Z, is the atomic number of atom r. Proceeding 
with the assumption that the scattering curves for the i 
atorns all have the same shape as that of atom r, we 
may write 


fi(hko) = Zg(hko), (11) 
and Eq. (9) becomes 
(Fino) i/fe(hko)=Z,/Z, cos2e(hxi+ky:). (12) 




















Fic. 1. Projection of the wave (Szy)aco, showing the wave crests 
only. The symbols used have the usual crystallographic meanings: 
a and b are the unit cell lengths, 8 the interedge angle of the unit 
cell, dixo the interplanar spacing and / and & the Miller indices. 


- Wako gives the direction of the wave front. For the wave shown 


Chko = 0. 


Making the definition 


(Gro) i= Zs/Z, cos2(hxi+ ky), (13) 
then 
M/2 
Fano/2f(hko)= ¥ (Gite): (14) 
830 


The wave (Gis); is shown in projection in Fig. 2, 
Comparison of this with Fig. 1 for (S.,)aze shows that a 
one-to-one correspondence between the two waves 
exists, except that the anio of (Szy)axo has no counterpart 
in (Grio)i. The wave-length \,* for (Gio); is given for 
the general case by 

sin p* 
a (15) 


[ (x,/a*)?+ (y,/b*)?— (2xy,cosB*/a*b*) }t 
and the direction of the wave front by 


y*=arc sin(x,d,*/a*). (16) 








Again, the type of construction given in Fig. 2 may be 


profitably employed for the graphical determination of 








- = 


} ANS / 
|_. iN 


Fic. 2. Projection of the wave (Gro)i, showing the wave crests 
only. The symbols used have the usual crystallographic meanings: 
a* and b* are the reciprocal cell lengths, 8* the interedge angle of 














' the reciprocal cell, d;* the wave-length of (Gixo):, and & and k the 


Miller indices. 


d;* and y,*. These latter quantities, together with Z,/Z, 
then completely define (Gnio). 

If the computation defined by Eq. (14) is carried out 
according to the scheme outlined above, i.e., with h 
and k as the continuous variables in reciprocal space, 
using the constants x; and y; for one unit cell only, such 
a projection will exhibit measurably large values of Fito 
at fractional, as well as integral values of / and &. 
This projection will have the same distribution of Fic 
in reciprocal space as would be given by the x-ray 
diffraction pattern of the crystal containing one unit 
cell only, if this were experimentally realizable. How- 
ever, since only the values of Fy. at integral values of 
h and k are usually of concern only these need be 
measured. 

It is at least of academic interest, however, and 
possibly of practical interest in considerations of very 
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small crystals, to investigate briefly the method for 
obtaining the F),./2f,(hko) projection as a function of 
the number of unit cells considered. Figure 3 shows a 
calculation of Fhoo/fzn for n=1, 2, 3, and N (very large) 
unit cells of 8B— ZnS, as a function of h. It will be noted 
from this figure that although the function Fhoo/ fzn has, 
for m small, non-zero values at non-integral values of h, 
the values of Fioo/fzn at integral values of h are inde- 
pendent of n. 

Taking the above facts into consideration, we may 
recast Eq. (14) in the form 


nM /2 


F ko, ‘2f-(hko) = 1/n »e (Ginko) i- (1 7) 


The summation limits are then taken to mean that 
the atomic coordinates of all the atoms (not related by 
the center of symmetry) of the » unit cells are included 
in the summation. An example will make this clear: 
consider an atom in the unit cell containing the origin, 
the first unit cell; its coordinate in the x direction is xp, 
while the equivalent atom in the mth unit cell along 
the positive x direction will have the atomic coordinate 
(n—1)+<2,y. Thus, the coordinates of all the atoms in- 
volved are obtained by adding the appropriate integers 
(negative or positive) to the coordinates of the atoms 
in the first unit cell. 


RECIPROCAL LATTICE PROJECTIONS FROM 
STRUCTURES NOT HAVING A 
CENTER OF SYMMETRY 


The computation of the crystal structure factor for 
the centro-symmetric case has been given in some detail 
above. For those structures not having such a center, 
even in projection, a quite analogous treatment may be 
made. The crystal structure factor for the general case 
is given, in the present notation, by 


M 
F ro/ fe(hko) => (Z;/Z,)e2*irzittve | (18) 
i=1 


or 


M 
Fiko/f(hko) => {(Z;/Z,) cos2a(ha;+ky,) 


i=1 


+ i(Z;/Z,) sindw(hx;+ky;:)}. (19) 
We may rewrite this as 
Frko/ f-(hko) =C(hko)+-1S(hko). (20) 


The sums C(hko) and S(hko) may be plotted separately. 
Simple vector addition, on the Argand diagram, of the 
real and imaginary contributions to each Fyxo will then 
yield | Frio| / fe(hko) and hko-» 
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GENERAL CONSIDERATIONS 


So far in the discussion presented, we have neg- 
lected the fact that several of the devices under con- 
sideration+*** give for each wave a sinusoidal distri- 
bution which follows the law! (neglecting phase shifts 
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Three unit cells ———— 

\ N unit cells 
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Fic. 3. Fhoo/fzn vs. h for B— ZnS, as a function of the number 
of unit cells used in the calculation. 


for the electron density case) 


sin 
lanl + lanl lo, (21) 
cos 
rather than 
sin 
onl] | (22) 
cos 


For the projections under consideration, the effect of 
the term >~,|a@,| is to introduce a constant background 
which may be quantitatively subtracted from the result- 
ing projection. 

The factors which affect the accuracy of projections 
of the type p’(x, y), made with the “sand machine,” 
have been discussed previously.‘ These same considera- 
tions apply to projections of the type presently con- 
sidered, when made with this device. - 

It should be pointed out that, in general, the param- 
eters x; and y; are irrational numbers. It is quite easy 
to prepare wave-length templates appropriate to these, 
within the necessary accuracy, for the “sand machine.” 
It would appear that this would be less convenient in 
the case of the other published methods, cited above. 
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Antennas and Open-Wire Lines. Part I. Theory and 
Summary of Measurements* 


RONOLD KING 
Cruft Laboratory, Harvard University, Cambridge, Massachusetts 
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Theoretical models of antennas driven by hypothetical genera- 
tors are contrasted with antennas driven from transmission lines. 
The coupling between an antenna and a two-wire line driving it, 
as well as end-effects on the line, are considered in their relation- 
ship to the definition of an independent impedance for the an- 
tenna. It is concluded that a physically significant impedance can 
be defined only outside a terminal zone near the junction of load 
and line. The derivation and approximate solution of generalized 
transmission-line equations that take account of coupling-effects 
and end-effects are given for (1) antenna as end-load in the plane 
of a two-wire line with and without a high impedance-stub 
support; (2) antenna as center-load in the plane of a two-wire line 
symmetrically driven at each end; (3) antenna as end-load in the 


plane perpendicular to the line with stub support. Theoretical 
predictions of measured impedances defined outside a terminal 
zone on a feeding line in these several arrangements may differ 
greatly from one another and from the ratio of scalar potential 
difference to current for an isolated antenna with gap. These 
differences are due primarily to capacitive end-effects and coup- 
ling-effects of which account may be taken by a lumped positive 
or negative capacitance across the load terminals of an ideal trans- 
mission line. Inductive effects are also involved but are relatively 
small. Methods for evaluating the elements in a terminal-zone 
network are given. The theoretical results are verified experi- 
mentally with good agreement in all cases. 





1. INTRODUCTION 


HE essential quality of a radiating circuit is to 
maintain an electromagnetic field at distant 
points. Such a field is due to currents strictly in all 
parts of the circuit, but the contributions due to cur- 
rents in the several parts may vary widely. In some 
circuits, the currents in all parts contribute significantly 
so that the entire circuit is both antenna and complete 
transmitting system. In others, the contributions by 
currents in most of the circuit may cancel and be neg- 
ligible compared with the contributions from currents 
in a specified part, that is designated as the antenna; 
obviously, this part is not a complete transmitting 
system. 
The electromagnetic energy-transfer function (de- 
fined as the integral, 


T= f (n’, S)de, 
Z (closed) 


over a completely closed surface of the normal com- 
ponent of the real part of the Poynting vector) cannot 
be used to define or even locate an antenna. For, if the 
surface of integration, 2, in (1.1) is progressively con- 
tracted and moved to maintain the same value of T 
as when = is a great sphere, it is always the generator 
and never the antenna that is located. The value of T 
when > encloses an antenna but not a generator (as 
when the antenna is driven from a two-wire line) is the 
net energy dissipated in the antenna. 

Theoretical studies of the impedance of straight 
antennas and of its practical significance have been 
confused by a failure to distinguish between (1) the 


(1.1) 


* Presented, Joint Meeting, U.R.S.I-I.R.E., Washington, 
D. C., May 4, 1948. Abridged from R. King, “Theory of antennas 
driven from a two-wire line,” Cruft Lab. Tech. Rep. No. 41 (June 
5, 1948). This research was supported in part by the Navy De- 
partment (ONR) and the Signal Corps, U. S. Army under Con- 
tract NSori-76, T.O. 1. 
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physical problem of an antenna as an actual circuit 
element in a complete system that also includes a 
transmission line and a generator with its coupling 
network, and (2) the analytically convenient abstraction 
of a rotationally symmetrical conductor that is itself a 
complete transmitting system by virtue of a fictitious 
source of e.m.f. at its center. 

If an antenna is the load on a long transmission line 
(Fig. 1), the physically significant and practically useful 
impedance is measured between two rings, C and D, 
around the surfaces of the two conductors at a cross 
section far from both the generator and the load. Each 
of the rings is an equipotential contour, and on a 
balanced line the cross sections they bound are traversed 
by equal and opposite currents. If C and D are far from 
both load and generator, and the transmission line is 
straight, there is no significant transverse component, 
A,, of the vector potential, A, along any contour ¢ in the 
transverse plane containing C and D, so that the trans- 
verse electric field is given by 


E.= = 06/dt— jw: = —d¢/dl, (1.2) 
where ¢ is the scalar potential. With (1.2), the voltage 
across CD is defined by 


c 
Voo= oo—be= f Edt. 


D 


(1.3) 


Only when the line integral of the electric field is sensi- 
bly equal to the scalar potential difference, as in (1.3), 
is it possible to introduce a unique voltage independent 
of the path of integration. 


Zpoc= Voe/Ic=Voc/Ip. (1.4) 


This impedance is independent of the distant generator, 
so that the actual oscillator with its coupling network 
may be replaced by an impedance in series with a fic- 
titious, localized e.m.f. that maintains the same voltage 
across CD, 
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If the antenna is assumed to be driven at its center 
by a fictitious, self-contained generator, the im- 
pedance of the antenna must be defined with respect 
to terminals that are a part of the generator and involve 
its structure and dimensions. A reasonable mathe- 
matical problem is obtained only if the properties of 
the generator are idealized to yield convenient boundary 
conditions. Two types of fictitious generators have been 
used: The “‘conducting-slice” generator shown in Fig. 2 
and the “dielectric-gap”’ generator of Fig. 3. 

The “conducting-slice” generator is an e.m.f. concen- 
trated at the center of a continuous, straight conductor 
or at the base of half of this conductor erected vertically 
on an infinite, perfectly conducting plane. It is the 
impedanceless and dimensionless generator familiar to 
electric circuit analysis. Since the conductor is un- 
broken, the e.m.f. must be maintained by forces acting 
from the outside, much as in a conventional alternator in 
which part of a continuous conductor is in a varying 
magnetic field. In practice, such forces can be main- 
tained only along finite sections of the conductor, and 
it is in assuming a potential difference maintained across 
an infinitely thin ring or belt that the physically possible 
is violated to achieve mathematical simplicity. Im- 
pedances for cylindrical antennas driven essentially by 
this type of fictitious generator are given in the litera- 
ture.'-* It is shown later that the impedance seen by a 
“conducting-slice” generator can be determined ex- 
perimentally, and hence given operational significance, 
by measuring the physically meaningful impedance 
defined in (1.4) for a number of different spacings of a 
two-wire line and then extrapolating these results to 
zero line-spacing. 

In the “dielectric-gap” generator the fictitious e.m.f. 
is across a dielectric gap between the halves of the 
antenna or between a vertical antenna and an infinite, 
perfectly-conducting, horizontal plane. This voltage is 
assumed to be maintained by rotationally symmetrical 
forces within the gap. In practice, a voltage can be 
maintained across a dielectric gap only by an externally 
attached transmission line, since generators built inside 
the antenna must be connected to its terminals by 
current-carrying conductors which cross the gap. The 
voltage across the gap may be defined either as the line 
integral of the electric field along a specified path across 
the dielectric (Fig. 3a), or as the scalar potential differ- 
ence (Fig. 3b). Since the vector potential necessarily 
has a non-vanishing tangential component, the two 
definitions are not equivalent. Infeld‘ has discussed the 
former definition in terms of an assumed, rotationally 





1 E. Hallén, “Theoretical investigations into the transmitting 
and receiving qualities of antennas,” Nova Acta (Upsala) 11, 
(November, 1938). 

*R. King and F. G. Blake, Jr., “The self-impedance of a sym- 
metrical antenna,” Proc. I.R.E. 30, 335-349 (1942). 

*R. King and D. Middleton, “The cylindrical antenna, current 
and impedance,” Quart. App. Math. 3, 302-335 (January, 1946); 
pp. 200-201 (July, 1946). 

* LL. Infeld, “Influence of width of gap upon theory of antennas,” 
Quart. App. Math., 5, 113-132 (July, 1947). 
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symmetrical field across a gap that is a physical ab- 
straction. King and Winternitz® have used the latter 
definition in terms of scalar potential difference. Note 
that the ratio of voltage across the gap (using either 
definition) to the current into the antenna is not a 
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Fic. 1. Antenna driven from two-wire line. 


physically meaningful or experimentally measurable 
impedance, so long as the current and voltage are due 
to generators with fictitious properties. However, 
since the ratio of scalar potential difference to current 
is not limited to the particular electric field and path’of 
integration appropriate to a fictitious generator in the 
gap (as is the line integral of the electric field) it may 
be used as a boundary condition at the load point of an 
open-wire transmission line actually driving the antenna 
by maintaining this potential difference and current. 
This is a step in the approximate determination of the 
measurable impedance defined in (1.4) as carried out 
in the following analysis. 


2. TWO-WIRE TRANSMISSION-LINE THEORY; 
GENERALIZED EQUATIONS 


The derivation of the transmission-line equations 
from fundamental electromagnetic theory involves 
assumptions and approximations that are often over- 
looked, but that are of great importance in the defini- 
tion, the theoretical determination, and the experi- 
mental measurement of impedances terminating the 
line, especially if these are antennas. One method of 
obtaining the transmission line equations is given in the 


A 


Fic. 2. Antenna with fic- 
titious, “conducting-slice” +¢¢ conducting 
mite. eed a Slice 
generator maintaining a@ . 
scalar potential difference. 








U 


literature® for a two-wire line of finite conductivity 
immersed in an imperfect dielectric. This method does 


5R. King and T. W. Winternitz, “Cylindrical antenna with 
gap,”’ Quart. App. Math. 5, 403-416 (January, 1948). 

®*R. W. P. King, Electromagnetic Engineering, I, (McGraw-Hill 
Book Company, Inc., New York, 1945). Chapter VI. 
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not separate effects due to non-uniformity of current on 
the line and radiation since both effects are ultimately 
neglected. In order to analyze conditions near the end 
of a line terminated in an antenna, it is necessary to 
take account of the non-uniformity of current. 

The electromagnetic foundations of the transmission- 
line equations are the defining equation for the scalar 
potential ¢, and the equation of continuity for poten- 
tials. If the two-wire line is parallel to the x axis in the 
xz plane, with the x axis midway between its two con- 
ductors, the following forms are involved: 


E,=2'l,= —0/dx— jwAz, (2.1) 
8A,/dx+0A,/dy+0A,/d2+ jB2¢/w=0, (2.2) 


where z' is the internal impedance per unit length and 
Bo=2x/Xo. The components A,, and A,, and their 
derivatives 0A,/dy, 0A,/0z, are zero everywhere on the 
two conductors if the line is infinitely long or has open 
ends. A, and A, and their derivatives may differ from 
zero near the end of a line terminated in a coil, a straight 
wire, an antenna, etc., in which currents are perpen- 
dicular to the x axis. Hence, 0A,/dy and 0A,/dz may 
not be omitted from (2.3) unless shown to be zero or 
negligible, if the resulting equations are to be valid 
near the ends of the line. 

Let (2.1) and (2.2) be written successively with sub- 
scripts 1 and 2 for the two conductors of the line and 


A 


S 


dielectric 
gop Fic. 3. Antenna with fic- 
E Foust Pe 4 es = titious “dielectric-gap” gen- 


erator maintaining in (a), 
a rotationally symmetrical 
electric field, in (b), a 
scalar potential difference. 


mm 
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a, 
(a) (b) 


subtracted to give: 
— (0/0x)(i— 2) = 21'T12—22'Toz+ jw(Aiz—A2z); (2.3) 
aA, aA, aA, aA, 

Tet Gd 

, aA, aA, 
-| dz ) -( Oz 

If the potentials and their derivatives in (2.3) and (2.4) 

are defined on the equipotential surfaces of identical 

conductors 1 and 2 with centers at (x, 0, 6/2), (x, 0, 


—b/2), and the impedances terminating both ends of 
the line are symmetrical with respect to the xy plane, 

















Be? 
) |- j—(¢i—¢2]. (2.4) 
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it follows that 
gii=giezi/2, I.=—l,=—1,, (2.5) 


where z’ is the internal impedance per loop unit length 
of the two-conductor line. The x component of the 
vector potential is an odd function of z calculated from 





1 e~ 1Bok 
A,(x)=A,= f 1,,——dx’ , (2.6) 
R 


4rvo 


where R=[(x—2x’)?+(y—y'?+(c—2’)?]! and = 
=10'/4mr meters/henry. If there is a 2 component of 
current in the symmetrical termination, it must be an 
even function of z, and hence the z component of the 
vector potential must also be an even function. It 
follows that, (0A,/0z),;— (0A ,/02)2=2(9A,/dz);. Note 
that (0A,/dz),-0= —(0A,/dz)0=0. 

For terminations that are symmetrical with re- 
spect to the xz plane as are here required, (0A,/dy); 
=(0A,/dy)2. For the antenna of Fig. 1, A, and its 
derivatives are zero everywhere. With the notation** 
V = ¢1—¢2= 261; We= Aiz— Aor= 2A 1, (2.3) and (2.4) 
become: 


— (dV /dx)=2'l.+joW;; (2.7a) 
— (OW ,,/dx)=2(0A,/d2):+j(Bo?/w)V. (2.7b) 


Differentiation of (2.7a,b) and separation of variables 
gives 


(PW ,/dx*)+BeW. 
= j(Bo?/w)z‘Iz—2(8/dx)(0A./0z), (2.8a) 


(PV /dx*)+ BeV = —2'(0I,/0x)+ j2w(GA 2/92):, (2.8b) 
where 


W,=W,(x) =2A1,(x) 


1 8 e1BoR. g—iBoR, 
= f 1e')(“ — )ax (2.9) 
2p 9 R Ry 


Ra=[(x’—x)*+a?}!;  Ro=[(x’—2x)?+8]!. (2.10) 


The length of the line is s. In expressing W, as in (2.9) 
it is implied that the only significant x components of 
current are the currents in the transmission line itself. 
If there are x components of current in the termina- 
tions, these are assumed to be so oriented and arranged 
that the resultant contribution to A, along the conduc- 
tors of the transmission line is zero. 

The z component of the vector potential in (2.8) is 
defined at the surface of line conductor 1. It is 


e—iBoR 
J I,(2’)——dz’, (2.11) 
R 


** Conventional transmission-line theory arbitrarily assigns 4 
positive direction to its current that is the opposite of that 
logically derived from the field equations. Hence, it is necessary 
to set V=¢1—¢» instead of V=¢2— qi, as is usual. 








A,(x)=A,= 
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where R is the distance between a point (x, 0, 5/2) 
where the vector potential is evaluated to the element 
dz’ at a point (x’, y’,2’) in the termination where 





I,(2’)#0. 
V=V(x)=Vi(«)+V r(x) (2.12) 
iw OW, 
ve inate (x) 
Bo? Ox 








1 8 e iBoR. e iBoR, 
“ f ae)( ~ )ae’ (2.13) 
2r€o 0 R, R 


b 


j2w OA, 1 e 1BoR 
V r(x) =— ‘)-—f q(z’)———dz’, 
Be X\ dz 27re€o R 





(2.14) 


where g(x’) and q(z’) are the charge per unit length in 
the line and in the termination, respectively. With this 
notation, (2.8a,b) may be written as follows: 














PW (x .. Bo : OV r(x ) 
BeW (x)= is ‘T,(x)+ | (2.15a) 
P Ox 
dV (x) 01 ,(x) 
+ BV (x)= —2* + BePV r(x). (2.15b) 
s ox 


V(x) is the potential difference between the wires of the 
line at the cross section specified by the coordinate x, 
due to all the charges in the line and in the termina- 
tions; V r(x) is that part of V(x) due to the charges in 
the terminations alone. For a line with conductors that 
are sufficiently closely spaced, V r(x) is significant only 
over short distances near each termination. Over most 
of the line V 7(x)=0; V(x)=Vz(z). 

Equations (2.15a,b) are integrodifferential equations 
for determining the current J,(x). The scalar potential 
difference V(x) must be a solution of (2.15b), but if 
(2.15a) involves J, in the terminations through A(x), it 
is clear that an equation similar to (2.15a) in A,(z) and 
I,(z) must be obtained from the boundary conditions 
of the termination and solved simultaneously with 
(2.15a) for J,(x) in the line and J,(z) in the termination. 
Such an equation has been derived in the literature for 
an antenna.*® Clearly, an exact solution of such a pair 
of simultaneous integrodifferential equations is involved. 

For perfect conductors (for which 2‘=0) and at 
distances from the termination at which the terms in 
A,({x) are negligible, W(x) and V(x) individually satisfy 
the homogeneous wave equation with no restriction to 
make radiation negligible. Unlike the potential functions, 
T, (x) satisfies the homogeneous wave equation only if 
radiation is made negligible by imposing the following 
condition on the distance 6 between centers of the con- 
ductors of the transmission line: 


(2%b/ do)? = Bereb’<1. (2.16a) 


This restriction is sufficiently severe to make radiation 
negligible in determining the amplitude and the dis- 
tribution of voltage and current along the line even 
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- Fic. 4. Lumped networks for use with a fictitious uniform line 
in the terminal zone to simulate an actual line with variable 
parameters. 


when not loaded. For a line with a load, and even for 
determining the approximate distribution of voltage or 
current on an unloaded line, (2.16a) is more severe than 
necessary. It is sufficient to require the average radiation 
resistance per loop unit length to be small compared 
with the inductive reactance per loop unit length. The 
condition is 








B2b?<«4x In(b/a). (2.16b) 
tty 
walle 
2c, s Zy 
: +z 
seid ‘ 
ite 


Fic. 5. Like Fig. 4b, but arranged for use with image-plane line. 


Let the variables w=s— x, w’=s—x’ be introduced in 
(2.9) and (2.13). Since the contributions to W.(w) and 
V(w) at a point w, due to currents and charges at w’, 
are significant only for distances (w’—w) that are not 
large compared with 3), it is a good approximation to 
expand the current J,(w’) in terms of J,(w), and the 
charge per unit length g(w’) in terms of g(w), and to 
retain only the first two terms. Thus, 


I,(w’)=I,(w)+(w’—w)(01,(w’)/dw’)wrew (2.17) 


q(w’) = q(w)+ (w’—w)(dq(w’)/dw")wraw. (2.18) 
Using the equation of continuity in the form 
jwq(w) = 01,(w)/dw, (2.19) 


the following expressions may be derived 


1 
W.(w)= —— 2(w)ko(w)+ jq(w)ki(w)/Bo] (2.20) 











” 1 \ba(to)-+ ki(w) dJ,(w) (2.21) 
Hu) = —| aw o(w i: a =f 
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where 


“7 (Ct 
ka(u)= f (~~~) 





4 eBook. — 1 e~ 1BoR,—] 
+f ( - ) (2.22a) 
0 Ra Ry 





A 
db 








Fic. 6. Measurable impedance must be defined at A’B’, not at AB. 


s 1 1 8 
ka(w)= bo f (w’—w)(—-— )du'+ sof (w’—w) 
0 b 0 


a 


c-iB Re] eit] 
x( _ ) a’ (2.22b) 
R, Ry 





The second integrals (2.22a) and (2.22b) take account 
of radiation. They are negligible when (w’—w) is large 
compared with 6 since then R,= R». On the other hand, 
when (w’—w) is of the order of magnitude of small 
multiples of 6 or less, BoR. and BR, are necessarily 
small compared with unity if (2.16a) is satisfied. If the 
exponentials are expanded in series, all terms below 
those of third power cancel. Hence, subject to (2.16a), 
these integrals, and with them radiation, are always 
negligible compared with the first integrals in (2.22a) 
and (2.22b). These latter are readily evaluated as 
follows: 


w w x x 
ko(w) = sinh—— sinh~'—-+ sinh—"—-— sinh (2.23) 
a b a b 
ki(w) = Bol (w*+ b*)!— (w*+-a*)! 
— (x°-+ 8?)!+ (a+a7)#]. (2.24) 


The greatest value of k:(w)/ko(w) is for w=0. Subject 
to (2.16a) or (2.16b), 


os kii(w) k?(0) Bo(b—a)? 








*(w) = = ~ - (2.25) 
’ ke(w) k(0) [lnb/aP 
Let the following symbolism be introduced: 
ko(w) 
lo*(w)= 5 
2rVv9 
1 
1o*=[lo°(w) use ~=— In- (2.26) 


Vo 
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( ) 2reo 
Co(w)= ; 
: ko(w) 
=[ca(w)] — (2.27) 
ia eatiiiaaal ~ In(b/a) 
Iy*(w)}! fo 
Rew) =| =—k)(w) ; 
Co(w) 2 
lt\? fo 5b 
R=(—) Pl wal (2.28) 
c r a 
Note that {9= 1/(vo¢o)!= 376.7 ohms and that 
Vo 4 1 1 
w=(“*) = = 29) 
€0 [Zo*(w)co(w) }* [lo*co]}! 
WwW 
Bo=—= of Lo*(w)co(w) = wf lo*co }}. (2.30) 


Vo 


With (2.23)-(2.27), and after eliminating q(w), (2.20) 
and (2.21) lead to 


OV (w)/dw+ Bop(w) Vi(w) 
I,(w)= , : 
‘+ julo*(w) 
a 





(2.31) 












































Fic. 7. Distances involved in definition of potential functions on 
the surface of one line conductor at P; end-loaded line. 
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The generalized transmission-line equations for the 
voltage V(w) and the current J,(w) are (2.31) together 
with (2.15b) with w replacing x. Since 2‘ for good con- 
ductors is small, only an insignificant error is involved 
if the correction term with the factor p(w) in (2.31) is 
neglected in evaluating 2‘(07,(w)/dw), and *(w) is 
treated as constant in the differentiation. These ap- 
proximations involve errors only very near the ends, 
and they are tolerated only in the small term in 2’. 
Thus, 

Den (w) | z' aV (w) 


Ow ~ sit jlo! (w) Ow" 





(2.32) 


Let the complex propagation constant be defined as 
follows: 


— Bo?Lz'+ julo*(w) ] 


= = jwco(w)[2*+ julo*(w)]. (2.33) 
julo*(w) 


2 





y 


Since 2‘ is small, 

















3'+ julo*(w) = [z'+ julo® lL ko(w)/ko(~)]=2(w), (2.34) 
so that _ 
2= jweo(z'+ jlo’), (2.35) 
which is independent of w. 
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Fic. 8. Plot of function involved in determining the equivalent 
lumped capacitance of terminal zone with end-loaded line. 


Using (2.32)—(2.35) in (2.15b) with w=s—x replacing 
x, this becomes 





eV (w) 
——— 7’¢i(w) V (w) =0, (2.36) 
Ow? 
where 
V(w)—V r(w) V L(w) Vi(w) - 
¢1(w)= (2.37) 
V(w) Vw) ~ Vi(w)+ V r(w) 


is the ratio of the scalar potential difference across the 
line due to the charges on the line to the scalar potential 
differences across the line due to the charges on the line 
and on the termination. $,(w) reduces to unity except 
near the ends of the line. Moreover, near these ends, 
Vr(w) is essentially in phase with Vz,(w) since the 
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charges that. contribute significanily to Vr(w) and 
V1(w) cannot be separated more than a fraction of a 
wave-length. Accordingly, ¢:(w) is essentially real. 
It is regarded conveniently as a variable scale factor 
for the capacitance co(w) that appears in the expression 
(2.33) for y*. Thus, let 


cr(w)=co(w)di(w) (2.38) 

1(w)= 7’ oi(w) = jwcr(w)2(w), (2.39) 
so that (2.36) becomes 

[eV (w)/dw* ]— y?(w) V(w) =0. (2.40) 


With (2.34) and (2.37), (2.31) reduces to 
T,(w) = 1/2(w) {LOV (w)/dw]+ Bop(w)i(w)V(w)} (2.41) 


Equations (2.40) and (2.41) are the final equations from 
which the scalar potential difference and the current 
in the two-wire line must be determined. 

At distances from the ends of the line specified by 
the conditions, w*>>b’, 2°>>6b, it is readily verified that 
(2.40) and (2.41) reduce to the well-known formulas of 
conventional transmission-line theory. Also, 


W(w)=I,(w)lo*; Vi(w)=V(w)=q(w)/co. (2.42) 


Evidently, the distributions of current and axial vector 
potential difference are the same as are the distributions 
of charge per unit length and scalar potential difference, 
except near the ends of the line where these simple pro- 
portionalities are not true. 

The results thus far obtained imply a*<<6*. They are 
readily generalized to permit unrestricted values of 
a/b by substituting b*= (b/2)[1+(1—(2a/b)*)*] for 5.’ 
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Fic. 9. Plot of equivalent lumped capacitance Cr, for end- 
loaded line for two lengths d of the terminal zone as a function of 
b/a. Broken line curve is Cr, computed from approximate formula. 


7 See reference 6, p. 440. 
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Conventional line theory assumes the consequences of 
w >}, 22>0 to be valid for all values of w and x. That 
this leads to serious errors whenever {ob is not extremely 
small is shown in the following sections. 


3. APPROXIMATE SOLUTION OF THE GENERALIZED 
TRANSMISSION-LINE EQUATIONS 


The generalized transmission-line equations that are 
to be solved for the scalar potential difference V and 
the current J, are (2.40) and (2.41). Their solution 
cannot be carried out directly since y*(w) is not a 
constant parameter independent of w in the terminal 
zones defined by 


0Sx=d~10b; 0=wSd~10b, (3.1) 


but is a complicated function of x or w, different for 
each termination. However, ?(w) reduces to the con- 
stant +’ at all points outside the terminal zones. For 
distances b between the conductors of the line that 
satisfy B,*b*<1, the length d of the terminal zone does 
not exceed small multiples of 6. Accordingly, if the solu- 
tion of (2.40) along the principal line outside the 
terminal zone is expressed as follows, 


V(w)=A coshyw+B sinhyw, (3.2) 


A, B, and y are constants outside the terminal zones. 
In the terminal zones a solution in the form (3.2) is 
not helpful, since A, B, and y would have to be com- 
plicated functions of the variable w (or x). This suggests 
the following method for determining an approximate 
solution at all points outside the terminal zones. The 
physical basis of the method is straightforward. It in- 
volves reducing (2.40) and (2.41) to conventional form 
at all points along the line by replacing *(w), z(w), and 
and p(w) by the constant values these functions have 
outside the terminal zone, viz. y?, z, and zero. This 
substitution is exact outside the terminal zone; it is 
increasingly in error as the terminal zone is entered and 
the terminated end of the line is approached. This error 
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Fic. 10. Antenna as center load on two-wire line symmetrically 
driven from both ends. 


occurs over a distance short compared with the wave- 
length if bd<X so that it is a good approximation to 
compensate for the error with appropriately defined 
lumped elements connected in parallel or series with the 
actual terminal impedance Z;. This impedance is de- 
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Fic. 11. Distances involved in definition of potential functions on 
the surface of one line conductor at P; center-loaded line. 


fined as the ratio of scalar potential difference to cur- 
rent; (the discussion is confined to the load at w=0 or 
x=s with the understanding that a parallel procedure 
applies at x=0). In other words, the analytical method 
is equivalent to replacing the actual terminal zone with 
its terminating impedance Z; by an equal length of 
conventional line with constant parameters and a ficti- 
tious, apparent load Z,, that is a combination of 
lumped elements in parallel or series with Z;. (Note that 
Z.a is the impedance actually measured on the line at 
do/2 from the load, if the line is lossless.) These lumped 
elements, C7 and Z7= jwL7 are defined as follows: 


d d 
Cr= f [er(w)—co |\dw= f [co(w)oi(w)—co dw. (3.3) 


That is, a constant-distributed capacitance per unit 
length, co, and a lumped capacitance, Cr, at w=0 are 
made approximately equivalent to a distributed capaci- 
tance per unit length cr(w). Since Cr takes account of 
the coupling between the termination and the line, it 
can be evaluated only if the structure of the termination 
is given. 


d 
Zr= f [s(w)—z ]dw. (3.4a) 


Since a small error in 2‘ is insignificant if confined to the 
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terminal zone, (3.4a) is practically equivalent to 


Lr= f [1o*(w) —Lo* |dw. (3.4b) 


Thus, a lumped inductance Z7 must be assumed con- 
nected in series with Z; at w=0 if the constant param- 
eter z is substituted for the actual, variable parameter 
2(w). 

The lumped series inductance defined in (3.4b) is 
readily evaluated. Near the end w=0 or x=s, x is very 
great compared with 6; furthermore, d>0*. With 
(2.23) the result is 


Lr=—(b—a)/2». (3.5) 


Since the lumped network of Cr and Lr is designed 
to make the terminated transmission line like an 
infinite line, it follows that p(w) must assume the ap- 
propriate value, viz. zero. Accordingly, for determining 
currents and voltages at points outside the terminal 
zone, the actual line with variable parameters may be 
replaced by a fictitious line with constant parameters 
and a modified termination. The two lines are the same 
outside the terminal zone, and they maintain the same 
ratio of scalar potential difference to current across the 
load. The equations for the fictitious line are the con- 
ventional ones, viz., 


[2?V(w)/dw* |—y?V(w)=0 (3.6) 
T(w)=(1/2)LdV(w)/dw]. (3.7) 


It may be concluded that for determining the scalar 
potential difference and current for a line outside the 
terminal zone, or for determining the theoretical load 
admittance Y,, it is possible to treat the terminal zone 
as though conventional line theory with constant 
parameters applied, provided Y; is treated as though 
connected in parallel with a lumped capacitance Cr and 
in series with a lumped inductant Zr as shown in 
Fig. 4a or 4b. For use in conjunction with an “image 
line” the circuit of Fig. 5 is convenient. It is equivalent 
to Fig. 4b but with elements so arranged that they are 
symmetrical with respect to a plane midway between 
the wires of the line. 

The lumped series inductance L7 has been shown to 
be negative. It compensates for the use of a constant 
inductance per unit length /)* in the terminal zone in 
place of the variable )*(w). The lumped parallel 
capacitance Cr may be positive or negative, depending 
upon the structure of the termination. For an antenna 
it is negative. It corrects for the use of a constant 
capacitance per unit length co in the terminal zone 
instead of the variable cr(w). The function of cr(w) 
includes the variation in capacitance per unit length 
near the end of the line as contained in co(w) and, in 
addition, the effect of coupling between the line and 
the termination. (Note that there is no contribution to 
k*(w) due to coupling only because it has been assumed 
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that the termination is so designed that currents in it 
are either perpendicular to the line or include equal and 
opposite elements with canceling effects.) 


4. ANTENNA TERMINATING TWO-WIRE LINE 


The rigorous one-dimensional analysis of the circuit 
consisting of a completely balanced two-wire line 
terminated in a symmetrical center-driven antenna 
(Fig. 6) requires the solution of simultaneous integro- 
differential equations for the distributions of current 
both in the antenna and in the two-wire line. Since the 
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Fic. 12. Like Fig. 8 but for center-loaded line. 


antenna and the line are coupled over distances along 
the atenna and along the line comparable with small 
multiples of the finite separation 5 of the conductors of 
the line, the definition of an input impedance Zaza 
(Fig. 6), that is characteristic of the antenna alone 
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Fic. 13. Like Fig. 9 but for center-loaded line. 


and independent of the feeding line, is not possible. 
An impedance in the sense of conventional transmission- 
line theory can be defined only for a sufficiently long 
section of line terminated in the antenna. Thus, Zp-4’ 
can be defined if the distance A’A in Fig. 6 is large com- 
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Fic. 14. Equivalent terminal-zone networks for center-loaded line. 


pared with the spacing 6. Note, however, that the im- 
pedance Zz-4 of the section of line to the right of 
B’A’ is not independent of the presence of the line to 
the left of B’A’. The assumption is implicit in the 
definition of every impedance on a transmission line 
that the uniform line continues in both directions from 
the points of definition B’A’. 

In order to obtain an approximate expression for the 
impedance of an antenna terminating a two-wire line, 
the method described in the preceding section may be 
used. In effect, this considers the characteristics of the 
section of line to the right of B’A’ (Fig. 6) in three 
parts that depend primarily upon (1) the transmission 
line, (2) the antenna, and (3) the coupling between 
them. Since significant coupling is confined to distances 
from the junction of antenna and line that are of order 
of magnitude 5b, the condition {,*b*<1 assures essen- 
tially near-zone coupling. However, conditions of sym- 
metry and perpendicularity eliminate inductive coupling 
between the currents in the line and in the antenna, 
so that the principal approximation is to replace the 
distributed capacitive coupling between line and antenna 
by a concentrated coupling in the form of a lumped 
capacitance in parallel with the antenna at its junction 
with the line. This is included in Cr. In this way, the 
transmission-line problem (1) and the antenna problem 
(2) are made effectively independent of each other ex- 
cept for the common potential difference across the 
terminals at the junction points and the requirement 
of continuity for the current at these terminals, 
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The determination of the lumped capacitance Cr 
defined by (3.3) involves the prior evaluation of the 


function 


[co(w)p(w)—co]=[gr(w)/V(w)]—co (4.1) 
V(w)=Vi(w)+Vr(w), (4.2) 


where V(w) is the contribution to the scalar potential 
difference V(w) across the line due to the distribution of 
charge gz(w) on the line alone, and where V 7(w) is the 
contribution due to the distribution of charges on the 
termination—in this case the antenna. The function 
Vz(w) is defined by (2.13) with w substituted for x; 
(Fig. 7). The function V r(w) is defined, in general, as 
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Fic. 15. (a) Stub-supported line; (b) equivalent network for 
terminal zone; (c) approximate, more convenient replacement 
of (b). 
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Fic. 16. End-loaded line with antenna perpendicular to the plane 
of the line; equivalent terminal-zone network. 
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and where gq(z’) is the charge per unit length on the 
upper half of the antenna; q(z’) is an odd function of 2’. 


The exact evaluation of Vz(w) and Vr(w) is intri- . 


cate. Fortunately, subject to the condition 6,°b’<1, 
significant contributions to the integrals in (4.3) and 
(4.5) come only from the integrands when evaluated 
at distances from the junction that are not large com- 
pared with b. Accordingly, the exponentials may be 
replaced by their leading terms of unity, and the upper 
limit of integration may be made infinite. 








Le a 

V.(w)=— ——— ibe’, (4.5) 
2reo 0 LR, Ry 
q of 1 

Vr(w)= ~f _——— |e (4.6) 
2mreowr2lRis Res 


The quantity gz is the average charge per unit length 
on the upper conductor of the line in the terminal zone. 
Similarly, gr is the average charge per unit length on 
the upper half of the antenna in a short region near the 
junction point. The approximations made in obtaining 
(4.5) and (4.6) are essentially the same as those implied 
in the derivation of the conventional line equations, 
but since V7(w) and V,(w) are involved here only in 
correction terms, these approximations are equivalent 
to higher order errors. The integrals may be evaluated 








[co(w)o1(w) —co/co= — 


sinh (w/a)—sinh—!(w/b)+ k, sinhb/(w?+ a”)! ]+-1n(b/a) 


Rot 











Fic. 17. Like Fig. 16 but with dielectric and stub supports 
adjusted to give maximum voltage across antenna terminals when 
the antenna is removed. 
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Fic. 18. Image-plane equivalent of Fig. 10. 
directly. 
qi dr | 
Vi (w)= ; Vr(w)= sinh~'————--_ (4.7) 
Cow 2r€o (w?+ a*)} 


In order to evaluate [¢o(w)¢1(w)—co] it is necessary 
to substitute (4.7) in (4.1). For convenience, let the 
charge-ratio factor k, be defined by k,=@7/gx. Then, 





(4.8) 


sinh (w/a)—sinh—(w/b)+k, sinh[b/ (w+ a?)*]—In(b/a) 


It is readily argued that k,=1 is a satisfactory ap- 
proximation. 

The function [co(w)¢i(w)—co]/co, as defined in 
(4.8), is shown plotted in Fig. 8 with b/a as parameter 
and w/b as independent variable in the range from 
w/b=0 to w/b=20 with k,=1. The area under each 
of these curves from w=0 to w=d multiplied by 
—cob= —bre,/Inb/a is the desired lumped capacitance 
Cr. (subscript e for end-loaded line) for each value of 
b/a and an appropriate choice of d. These areas, ob- 
tained by numerical methods, are shown in Fig. 9 as a 
function of b/a and for d=10b and d=20b. Since the 
values of —C r,/bco for these two choices of the length d 
of the terminal zone differ appreciably, it is clear that 
significant capacitive coupling between the antenna 
and the line exists for distances from the antenna 
exceeding 10 times the line spacing 6. This means 
that the representation of this distributed effect by an 
equivalent lumped capacitance at w=0 is necessarily 
only a fair approximation unless d can be taken quite 
large compared with 6 and yet be a small fraction of a 
wave-length. In order to make it possible to treat end- 
effects on lines at frequencies for which b/o=0.02, it 
IS necessary to restrict d to at most 10 or 0.2X9. Even 
for this length the assumption of equal average charge 
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concentrations per unit length on antenna and line 
(k,=1) is justified only near antiresonance. 

For the rapid evaluation of C7,, a closed formula is 
desirable. The shapes of the curves obtained in Fig. 8 
suggest 


@=«Kb d 
Cr=—cof ———dw=— 6 bK sinh—- (4.9) 
o (w?+k*b*)! kb 


K and k are determined by fitting — Kb/(w*®+ kb")! to 
[co(w)¢1(w)—co]/co in the range of w from zero to 
d=10b. This is done by requiring (4.9) to be satisfied 
exactly at w=a and w= 36. 

Curves showing the approximate function —C7./bco 
as evaluated from (4.8) are shown in Fig. 9 together 
with the same function determined from (4.9) by 
numerical integration. The agreement is seen to be fair, 
with the approximate curve falling between the curves 
for terminal zones of length d=10) and d=20b. Since 
the proposed representation of a distributed effect by a 
lumped capacitance is in any case far from exact, the 
values of Cr, taken from the simple integrated formula 
(4.9) are quite satisfactory. Note that the order of 
magnitude of Cr, is that of a section of uniform line 
of length equal to the line spacing b. 
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Summarizing the analysis of the two-wire line termi- 
nated in a center-driven antenna, the following may be 
said: the apparent impedance Z,, terminating the two- 


wire line—which is the actual impedance looking toward © 


the antenna from a point on the line (assumed lossless) 
one-half wave-length from the antenna—may be de- 
termined using the circuit of Fig. 4a or 4b if the ratio 
of scalar potential difference to current at the terminals 
of the antenna, viz., Z;= V3/I;s, is known. In the figures 


Lr=Lr.= —(b—a)/2rv (4.12) 


d 
Cr=Cri= f [co(w)oi(w) —co dw 
, = —Keob sinh-(d/kb). (4.13) 


The subscript e identifies the values for an end-loaded 
line with antenna in the plane of the line. 


5. ANTENNA AS MIDPOINT LOAD FOR 
SYMMETRICALLY DRIVEN LINE 


The circuit of Fig. 6, in which a center-driven an- 
tenna terminates a balanced two-wire line, has several 
practical shortcomings if interest lies in the accurate 
measurement of the impedance of the antenna. One 
of the problems is to support the antenna and line at 
or near their junction. A common method uses a piece 
of dielectric material. However, even the theoretical or 
experimental determination of the actual added capaci- 
tance due to a very thin dielectric support and not 
merely the apparent added capacitance is complicated. 
In usual practice only the apparent combined impedance 
terminating the line is of interest, and this is readily 
measured directly in each case. On the other hand, a 
complete, even though perhaps quantitatively only 
approximate understanding of the factors that con- 
tribute to the apparent impedance of a termination re- 
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Fic. 19. Possible method of driving antenna to approximate slice 
generator with coaxial line in (a) and image equivalent in (b). 
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quires a separation of the effect of a dielectric support 
from those effects characteristic of the configuration of 
conductors. 

A very flexible structure which provides both a means 
of support for the line at its junction with the antenna 
and complete symmetry is the circuit shown in Fig. 10. 
The antenna is attached to the center of a symmetrical 
line that is driven from each end by identical generators. 
The currents in the two driving lines are equal in mag- 
nitude and opposite in direction at corresponding 
points measured from the antenna. The charge per 
unit length at these points are the same in magnitude 
and sign. The currents J,(w=0) at the junction point 
w=0 combine into the antenna current /,(4) ; the charge 
per unit length on the antenna near the junction point 
is equal to the charges per unit length on each of the 
two line conductors very near the junction. 

If the impedance Z; of the antenna is defined formally 
as the ratio of the scalar potential difference V; at the 
junction point to the current /,(4), 


- V3 _ o(8)—9(—8) 
18) 1,(8) 





, (5.1) 


it follows that, 
Vs 


Z,=———_, (5.2) 
21, (w=0) 


where J, is the current in conductor 1 of the line. 








Fic. 20. Conventional 
method of driving antenna 
above image plane from 
coaxial line. 





Hence, if complete symmetry is maintained, each line 
may be treated as though terminated in an impedance 


Vs V; 


Zs= = : (5.3) 
1.(w=0) —I,(w=0) 





Therefore, the entire analysis of Sections 2-4 applies, 
except that the lumped elements L7 and C7 must be re- 
evaluated since the two parts of the line are coupled 
not only to the antenna but also to each other. 
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It is readily verified that 


Lre= ae (b—a)/mvo. (5.4) 


Note that Lc in (5.4) is the lumped inductance to be 
connected in series with each half of the symmetrical 
line with uniform constants. The subscript C identifies 
the center-loaded line. 

It remains to evaluate C rc. This is carried out as for 
the end-loaded line but with the following modified 
expressions for Vz;(w) and V7(w) which must now in- 
clude the effect of the additional feeding line (Fig. 11). 








wr 6 qg 
Vi(w)= an Fi —-— |e’ he (5.7) 
2re_wLR, Rs Teg a 
or i 1 
Vr(w)= tn f ———|d 
2re€o 2 Ris Roa 
q 
=—— sinkt—_——- (5.8) 
27é (w*+-a?)! 
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As before, kg=4,/,, where 9, is the average charge 
per unit length on the antenna near the junction with 
the line, and q, is the average charge per unit length on 
one conductor on each of the two lines. 

Substitution of (5.7) and (5.8) in (4.1) gives 


[co(w)gi(w) —co]/co 

k, sinh[b/(w?+ a*)*] 
2 In(b/a)-+-k, sinh"[1/(w?-+-a2)*]) 
Curves of [co(w)¢(w)—co ]/co for kg=1 are in Fig. 12. 





(5.9) 


Comparison with Fig. 8 shows that for the center- 
loaded line, [co(w)¢(w)—co]/co is approximately equal 
to its value for the end-loaded line except at very small 
values of w/b. The lumped capacitance Cre is given 
by (3.3). 

Curves showing Cr¢ as a function of b/a are in Fig. 13 
as evaluated by numerical methods with d=10b and 
d=20b. Comparison with Fig. 13 for Cr. shows that 
they are roughly equal. 
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Since the general shapes of the curves of [co(w)¢(w) 
— co |/co for the center-loaded line is the same as for the 
end-loaded line except near w/b=0, they may be 
evaluated as before by setting 


[co(w)o1(w) — co ]/co= —K/[(w/b)?+k}#. (5.10) 


K and are obtained by evaluating (5.10) at w=0 and 
w= 3b. Curves of —Cre/bco are in Fig. 13. 

To conclude the discussion of the symmetrically 
driven line, attention is called to the circuit of Fig. 14a 
which must replace the antenna as center-load if con- 
ventional transmission-line theory is to be a satisfac- 
tory approximation. In actual use the circuit of Fig. 
14b, which is equivalent to Fig. 14a, is convenient. 
Since no current crosses the central axis of symmetry, 
the two connections between the two loads 2Z; may be 
removed and each side treated separately. Alternatively, 
the circuit may be arranged as in Fig. 14c which is 
useful when an “image line” is used, as discussed in 
Section 6 and Part III. 

Since Lrc is exactly twice Lr, but is connected in 
series with 2Z; in Fig. 14b, the inductive corrections are 
the same for the center-loaded and end-loaded lines. 
Since Crc is approximately the same as Cr., Cre in 
parallel with 2Z, for the center-loaded line is effectively 
double Cr, in parallel with Z; for the end-loaded line. 
The apparent terminal impedance Z,, of a particular 
antenna must, therefore, be different when determined 
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for the antenna as end-load from its value for the same 
antenna as center-load unless Cr is negligible. 

Instead of a dielectric support, a so-called “stub” 
support consisting of a high impedance, antiresonant 
section of transmission line may be used (Fig. 15). In 
order to achieve a sufficiently high impedance, losses 
in the stub must be very small and, hence, the adjust- 
ment to antiresonance very critical. A practical pro- 
cedure is to determine accurately with a sensitive probe 
the location of the voltage maximum with the antenna 
absent, and then to attach the antenna exactly at this 
point. If the stub is so constructed that the impedance 
of the stub is very great compared with that of the 
antenna, the current in the stub near the junction with 
antenna and line is very small, so that the inductance 
per unit length of the transmission line differs only 
negligibly from that for the unsupported line and an- 
tenna. Hence, Lr=Lr,. However, the concentrations 
of charge on the sections of the stub near the junction 
point are by no means negligible compared with those 
on adjacent sections of the antenna and the line proper. 
When there is a maximum voltage at the terminals of 
the antenna, the distribution of voltage and charge 
along the transmission line including the stub is essen- 
tially the same as the distribution for the symmetrical, 
center-loaded line. Accordingly, Cr=Crc. However, 
in this instance the impedance terminating the line is 
not 2Z; but simply Z; so that the equivalent terminal- 
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Fic. 22. Theoretical resistance of antenna in parallel with lumped positive and negative capacitances; 2= 10. 
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zone network is in Fig. 15b rather than in Fig. 14b. 
Note also that Lr, replaces Lrc. No inductances are 
required on the right, since the current into the high 
impedance stub is always small. Since Lr, is small, 
the two capacitances may be combined as in Fig. 15c 
without significant error. The resulting circuit is the 
same as for the center-loaded antenna in Fig. 14b. 
Note that the capacitance C rc on the right in Fig. 15b 
is due entirely to coupling between the antenna and the 
stub, so that it can mot be compensated by adjusting 
the stub with the antenna removed. For this reason it 
must be included in the circuit even if the stub is 
adjusted experimentally to put a voltage maximum at 
the load terminals before the antenna is connected. 


6. ANTENNA CENTER-DRIVEN FROM A TWO-WIRE 
LINE WITH MINIMUM COUPLING; IMAGE 
METHOD OF DRIVING AN ANTENNA 
USING OPEN-WIRE LINE; ANTENNA 
PERPENDICULAR TO IMAGE 

PLANE : 


In order to eliminate all coupling between antenna 
and transmission line it is necessary to orient the an- 
tenna so that the vector and scalar potentials due to 
currents and charges on the transmission line vanish at 
all points along the antenna. In the end-loaded and 
center-loaded lines the tangential vector potential due 
to currents in the line is zero on the surface of the 
antenna since antenna and line are mutually perpendicu- 


500 eyes 


lar; however, the scalar potential on the surface of the 
antenna due to charges on the line does not vanish. 
In order to make this zero on the antenna, opposite 
points on the two conductors of the line must be made 
equidistant from each point on the axis of the antenna. 
If this is done, the contribution to the scalar potential 
at any cross section of the antenna by the charges in one 
conductor of the line is canceled by the scalar potential 
due to an identical distribution of charge of opposite 
sign on the second conductor of the line. One method 
of accomplishing this cancelation is shown in Fig. 16a 
for an antenna driven as an end load. (In this circuit 
there is no inductive coupling between the antenna and 
the line or the antenna and the short bridges connecting 
the antenna to the line so that the inductive end-effect 
is due exclusively to the absence of an infinite line just 
as in the more conventional circuit in Fig. 1. Therefore, 
the same value of Lr may be used.) 

In the conventional end-loaded line the correcting 
capacitance Cr combines the effect of non-uniformity 
in the capacitance per unit length of line with coupling 
between line and load. In the arrangement of Fig. 16 
this latter effect vanishes so that Cr is simplified to the 
extent that direct coupling between the two parallel 
conductors of the line and the antenna is eliminated. 
If the two short conductors connecting the ends of the 
line to the terminals of the antenna were exactly per- 
pendicular to the antenna, the scalar potentials on the 
antenna due to the equal charges of opposite sign on 
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Fic. 23. Theoretical reactance of antenna in parallel with lumped positive and negative capacitances; 2= 10. 
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them would cancel exactly. Evidently this is possible 
only in the idealized case of zero gap. For a finite gap 
the potentials due to the charges on the short connect- 
ing sections do not quite cancel, so that a small capaci- 
tive coupling between 4ntenna and the connecting 
conductors exists. Since each of the conductors is es- 
sentially of length 6/2 the magnitude of the coupling 
is very small and approximately equivalent to a small 
lumped capacitance connected across the terminals of 
the antenna. 

In practice the arrangement of Fig. 16 is not con- 
venient without a dielectric or stub support, of which 
account must be taken. Fortunately, and unlike the 
corresponding situation in Fig. 15, the presence of a 
stub in the form of a continuation of the line as shown 
in Fig. 17 involves no coupling to the antenna, although 
Lr and Cr are affected due to modification of end- 
effect on the line. However, since the antenna is not 
directly involved, it is possible to use a dielectric sup- 
port in conjunction with a stub so adjusted that in the 
absence of the antenna the impedance of the dielectric 
in parallel with the stub is very great compared with the 
maximum impedance of the antenna. Experimentally, 
this is readily accomplished by adjusting the length 
of the stub in the absence of the antenna until a voltage 
maximum occurs exactly at the dielectric support as 
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Fic. 24(a). Apparent terminal resistances of antennas in the 
plane of the two-wire line. Curve A, second-order King-Middleton 
theory for antenna with slice generator; curve E, theory cor- 
rected for end-load; curve S, C theory corrected for center-load 
or stub-supported end-load; curve M,, experimental data of K. 
Tomiyasu for stub-supported end-load. 
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determined from the main line. In this way end-effect 
and the presence of the dielectric and the stub are all 
combined in a manner that virtually eliminates them. 
If the antenna is attached after this adjustment has 
been made, the impedance measured for the termination 
is essentially the impedance Z; of the antenna. That is, 
Z.a= Zs, if the series inductance of the half-bridges is 
negligible, as is true when these are of sufficiently large 
cross section. 

The experimental study of antennas center-driven 
from a balanced two-wire line may be facilitated in some 
instances if use is made of the theorem of images* to 
provide equivalent structures which provide a shield 
for the observer. For the conventional end-load line 
(Fig. 6) or the center-loaded iine (Fig. 10) (but not for 
the circuits of Figs. 16 and 17), the plane z=0, per- 
pendicular to and bisecting the rectangle formed by the 
line, is a plane of symmetry for the entire circuit. 
On this plane the tangential component of the electric 
field is zero everywhere, so that nothing is changed 
electromagnetically if a perfectly conducting sheet is 
placed to coincide with the plane z=0. Such a sheet 
completely isolates the halves of the circuit and the 
fields associated with their charges and currents, so 
that one-half may be removed without disturbing the 
other. Each then consists of one conductor of the orig- 
inal two-wire line parallel to a perfectly conducting 
sheet at a distance }/2 and one-half of the original 
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Fic. 24b. Like Fig. 24a but for reactances instead of resistances. 


8 See reference 6, Section IV, p. 20. 
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antenna perpendicular to the sheet as shown in Fig. 18 
for the center-loaded line. 

Instead of using half of a two-wire line over a con- 
ducting plane, it is often convenient to use a coaxial 
line in the arrangement shown in Fig. 19. In so far as 
the electromagnetic field and distributions of current 
and charge in the half-space z=0 are concerned, the 
theorem of images states that the circuit of Fig. 19a 
is equivalent to the physically awkward, symmetrical, 
center-driven circuit of Fig. 19b. In order that the 
scalar potential difference across the double gap in the 


1800 


















































VARIATIONS IN APPARENT ANTENNA IMPEDANCE 7 
WITH LINE SPACING O/d5 3.9710 A 
b/r» 
need 0.0000 KING-MIDDLETON SECONO-ORDER \ 
— —— 0.0159 (CORRECTED) — v} 
Re |----— 0.0238 (CORRECTED) | Fs \\ 
1400 b ¥ae — 00556 \++ 
01032 |] \ | 
i i ‘| \ 
1200 A; T 
4 ‘* 1 | 
/ \ hi 
1000 | i | 
/ \ 
} \ 
200 AALA | ht \ 
\ 
\ 
600 a + 
¥\ | \\ 
‘ 
© oe Yh NI .\ \ 
oe | 4) 
: | | XA WN] 
Z \ N 
— VA lw 2 ol Sian 
‘ | x ~ Z 
| \: : | J TA 
¢/ 4 ; 
\ } 4 
° es Se f- 





' f7 ) ; | %y 
| Aer 











OT Ta 
aw | 7A WN AA 
of {| At 

| | lV i 
¥ é 



































“1000 L 





Fic. 25a. Apparent terminal impedances of antennas as center 
load obtained experimentally by P. Conley on image-plane line 
together with theoretical curves for antenna driven by slice 
generator. 


symmetrical circuit of Fig. 19b be twice the scalar 
potential difference across the gap in the circuit of 
Fig. 19a, the gap must be sufficiently short to make the 
following condition a good approximation : 


V;/2=4(6)—(0)= f " Bds> ju f "Adz. (6.1) 


For a sufficiently small gap it is reasonable to assume 
that little error is introduced if the gap is in the con- 
ventional form of a hole in the conducting plane as 
shown in Fig. 20 instead of a slit around the surface of 
the antenna. The requirement now is 


Vs/2=4(b)—4(a)= f * pdr > je f "Addr. (6.2) 
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Fic. 25b. Resonant resistances and lengths as determined ex- 
perimentally by P. Conley with curves extrapolated to theoretical 
values for slice generator. 


A more detailed investigation of the antenna driven 
from a coaxial line as shown in Fig. 20 is not con- 
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Fic. 25c. Antiresonant resistances and lengths as determined 
experimentally by P. Conley with curves extrapolated to theo- 
retical values for slice generator. 
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templated here. The present discussion is intended tances due to dielectric supports or binding posts, to tc 
merely to permit comparison of this arrangement as a_ correspondingly large negative capacitances encoun- ta 
more usual alternative to the open-wire line in measur-_ tered in the terminal zone networks of antennas as - 
ing impedances for cylindrical antennas. end-loads, both with and without stub supports, and ae 
antennas as center-loads. Since the difference between “ 
7. THEORETICAL AND EXPERIMENTAL the impedance Z; of the isolated antenna with gap ; 
IMPEDANCE CURVES 25=6 and the impedance Zp» of the same antenna with fi 
It has been demonstrated theoretically that the ap-- 5=9 is small, present purposes are adequately served a 
parent impedance of an antenna as a load on a trans-__ if the capacitance C is connected in parallel with Zo as t 
mission line may be determined approximately from obtained from the second-order King-Middleton ex- ” 
the theoretical impedance (defined as the ratio of scalar pansion.’ , 
In Figs. 21-23 are shown curves of Bo, Ro, and Xo in = 
eet parallel with a lumped capacitance, 
~ am fz * Rig* jx, | Cs= KBoh/w= Kh/t», (7.1) of 
“= where K is a constant parameter to which a range of al 
— eaten positive and negative values is assigned. Gp is un- gi 
Debts asethaetitee affected by a capacitance in parallel. The independent of 
_ variable is Bok, the radian half-length of the center- tv 
- driven antenna. If C is assumed to be a fixed constant, a 
hi 
aia 10,000 cl 
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Fic. 26. Theoretical and experimental impedance curves. Zo is A. 
; the King-Middleton second-order impedance for antenna with 100 “= 
slice generator; Zs; is the theoretical impedance for antenna with 80 co 
gap; Zsa is the apparent terminal impedance determined theo- fo 
retically using Zs; with terminal-zone network and measured e : 
experimentally by P. Conley on image-plane line. ie tr 
in 
potential difference to current) of an isolated antenna ™ 
with gap if account is taken of coupling- and end-effects se 
° ° ° . ° Qs 22tn 22 m 
near the junction points. This may be done with ap- e 
propriately determined inductances in series with the e “ 
= . eng . 05 1.0 1.5 2.0 2.5 3.0 3.5 4.0 m 
antenna and capacitances in parallel with it. Inductive A 
effects usually are quite small, capacitive effects may be a , . 
; taal A ding! h Saad f Fic. 27. Apparent terminal impedance of antenna in plane cu 
| surprisingly large. /Accordingty, the generat nature Ol perpendicular to two-wire line. Theoretical curves are King- So 
the possible differences between the apparent impedance Middleton second-order impedance, experimental points deter- , 
of antennas driven from various transmission lines ™imed by K. Tomiyasu using dielectric and stub supports ad- = 
justed to give maximum voltage at terminals of antenna when this 





using several methods of connection may be investi- was removed. 
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. gated by determining the effect on the admittance and N. 
impedance of the isolated antenna of a lumped capaci- {oh=wh/v is varied by changing the frequency in lin 
’ tance in parallel. By allowing the capacitance to assume w=2xf. The curves marked K=O are for the antenna (te 
4 a range of positive and negative values, a true picture alone. Note that a positive capacitance in parallel, as Cr 
is given of the measurable admittance and impedance with a dielectric support, leads to reduced magnitudes 7 

. under conditions ranging from large positive capaci- of resistance and reactance with shifts of antiresonance ma 
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to smaller values of Boh. Conversely, a negative capaci- 
tance in parallel increases the magnitudes of resistance 
and reactance very appreciably if K is not too great, 
and produces a shift of antiresonance to much greater 
values of Bok. It is clear that the nature of the ad- 
mittance or impedance of an antenna is greatly modi- 
fied by a capacitance in parallel. For example, the 
curves for K=0 would be quite inadequate to describe 
the apparent impedance of antennas and terminal 
zones for which K = —0.0003 or K=+0.0003. 

In Figs. 24a and 24b are shown the calculated re- 
sistance and reactance curves and one experimental 
pair® for a particular antenna and line. Curves A give 
the real and imaginary parts, R; and X;, of the ratio 
of scalar potential difference to current for an isolated 
antenna driven by a fictitious generator. Curves E 
give the apparent terminal impedance Zsa= RyatjX ea 
of the same antenna used as end-load in the plane of the 
two-wire line (Fig. 6). Curves S are for the same 
antenna and transmission line as E but with an added 
high impedance-stub support (Fig. 15). Similarly, 
curves C (which are the same as curves SS) are for the 
same antenna when connected as center-load on a line 
symmetrically driven at both ends, as in Fig. 10. 
Note that the curves are shifted to longer electrical 
lengths for antiresonance and greater magnitudes 
of impedances in proportion to the negative shunt 
capacitances of the appropriate terminal-zone networks. 

The dotted curves M, were obtained experimentally 
with the particular antenna for which the theoretical 
curves were evaluated ;? a high impedance-stub support 
was used as in Fig. 15. Note that agreement between 
the experimental curves M, and the curves A for the 
isolated antenna with fictitious generator is very poor. 
On the other hand, the agreement with curves S, ob- 
tained using the values of curves A in conjunction with 
the appropriate terminal zone network, is quite good. 
As predicted, the attempt to correct for a distributed 
coupling- and end-effect by lumped elements evaluated 
for antiresonant conditions (k,=1) is not quantita- 
tively exact. But whereas curves A are almost useless 
in attempting to predict the experimental curves M,, 
curves S are adequate for many practical purposes. 

In Fig. 25a are shown curves of R,q and X,q deter- 
mined experimentally for an antenna connected as 
center-load in the plane of a transmission line sym- 
metrically driven at both ends (Fig. 10) for four values 
of the ratio of line spacing to wave-length.” Theoretical 
curves of Ry and Xo (King-Middleton second-order) 
for the same antenna driven by a “slice” generator are 
also shown. Owing to a slight bow in the image plane 

°K. Tomiyasu, “Problems of measurement on two-wire lines 
with application to antenna impedances,” Cruft Lab. Tech. Rep. 


No. 48 (June 15, 1948). K. Tomiyasu, “Antennas and open-wire 
lines, part II, measurements on two-wire lines,” J. App. Phys. 
(to be published). 

* P. Conley, “Impedance measurements on open-wire lines,” 
Cruft Lab. Tech. Rep. No. 35 (March 18, 1948). P. Conley, 
Antennas and open-wire lines, part III, image-line measure- 
ments,” J. App. Phys. (to be published). 
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the two sets of experimental curves for the smaller 
values of b/d» involve considerable error for which an 
approximate correction was introduced.” In Fig. 25b 
are resonant and antiresonant lengths and resistances. 
The effect of the line spacing on the apparent impedance 
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Fic. 28. Theoretical antenna impedances for slice generator with 
coaxial-line measurements by D. D. King. 


is great. In Fig. 26 one of the pairs of curves for R,, and 
Xsq in Fig. 25a is reproduced together with the corre- 
sponding theoretical values of Ro and Xo, Rs and X; 
for the antenna with gap 26, and of R,. and X,q using 
the appropriate terminal zone network with elements 
computed from the constants of the particular line. 
The experimental curves are seen to be in very much 
better agreement with the theoretical curves of Ria 
and X,, than with Ry and X» or R; and X;. Note that 
the maxima of the experimental curves of R and X 
are relatively too high as seen from Fig. 25a. Actually, 
experimental and theoretical maxima are almost 
coincident if interpolated experimental values are used. 
Unfortunately, the theoretical values of R,g and X,q for 
the two larger values of b/d» cannot be determined due 
to the limitation of the theory to small values of b/Xo. 

As a final test of the new theory, Fig. 27 shows the- 
oretical curves of R; and X; for an isolated antenna with 
gap 26 and fictitious generator together with experi- 
mental points’ for an antenna connected as end-load 
pependicular to and in the neutral plane of the trans- 
mission line, as in Fig. 17. The stub and dielectric 
supports were adjusted to give maximum voltage at the 
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load terminals before the antenna was attached. The 
connecting half-bridges were of large cross section so 
that their series inductance was negligible. 

As predicted for this special case in which there is 
virtually no coupling between antenna and line, and in 
which line end-effects are compensated experimentally, 
the apparent terminal impedance Z,, is the same as the 
theoretical ratio of scalar potential difference to cur- 
rent, viz., Zs for the isolated antenna. 

For all other types of connection previously dis- 
cussed such an equivalence can exist only in the theo- 
retical limit as the line spacing, and with it the length 
of the terminal zone, approach zero. This limit, while 
physically unattainable, may be defined by graphical 
extrapolation as in Fig. 25b and so serve as an opera- 
tional definition of the idealized slice generator. It has 
already been shown" that for an antenna driven from a 
coaxial line in the arrangement of Fig. 20, the apparent 
terminal impedance Z,, is closely approximated by Zo 
in parallel with the lumped positive capacitance due to a 
thin dielectric wafer at the junction between antenna 
and line, provided the inner radius of the outer con- 
ductor of the coaxial line at this junction is a very 
small fraction of a wave-length. This is illustrated in 
Fig. 28 where compensated critical data obtained by 
King" are compared with theoretical values. 


8. CONCLUSION 


An approximate theory has been derived and verified 
for calculating the measurable, apparent impedance of 


1D. D. King, ““Measured impedance of cylindrical dipoles,” 
Cruft Lab. Tech. Rep. No. 2 (July 19, 1946); J. App. Phys. 16, 
(October, 1946). 





an antenna loading an open-wire line in several arrange- 
ments. Use is made of conventional transmission-line 
theory and antenna theory together with a properly 
designed terminal-zone network that takes account 
of transmission-line end-effects and coupling between 
antenna and line. In this way the King-Middleton 
second-order impedances as generalized to a finite gap 
by King and Winternitz are shown to lead to good 
agreement with impedances measured on open-wire 
and coaxial lines using widely different connections. 
These measured impedances differ greatly. 
Significantly, the new theory is not limited to an- 
tennas, but applies equally to all types of load. It con- 
stitutes an essential supplement to conventional trans- 
mission-line theory whenever an increase in the length 
of the line by an amount equal to its spacing produces 
more than a negligible change in impedance. For ex- 
ample, the new theory has been applied to open-end 
and bridged-end transmission lines to determine end- 
corrections. It may also be applied to discontinuities. 
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Potential Functions for a Thermionic Vacuum Tube 


E. C. OKREss 
Westinghouse Electric Corporation, Lamp Division, Bloomfield, New Jersey 
(Received February 21, 1949) 


Space charge free potential functions for a particular type of thermionic vacuum tube are derived and 
represented graphically for four sets of boundary conditions. The geometry under consideration comprises 
two spaced parallel equipotential metallic planes between which are three symmetrically interposed planes 
of spaced metallic strands. The strands are mutually parallel in their respective planes, though the strands 
in planes adjacent to the equipotential surfaces are both orthogonal to those in the central plane. This 
structure is divided principally into two planes, due to the orthogonal relationship between the strands in 
the central and adjacent planes. These two planes are each conformally transformed such that the analysis 
thereof yields the desired potential functions for the original structure, neglecting edge effects. Finally, 
the theoretical results are applied to a typical commercial thermionic vacuum tube. 


I. INTRODUCTION 


HE object of this exposition is primarily the 

derivation of potential distribution functions by 

the application of conformal transformations for a type 
of geometry illustrated essentially by Fig. 1. 

Although, either (a) the experimental method, in- 

volving electrolytic tank field plotting,’ or the (b) re- 


1L. M. Myers, Electron Optics (D. Van Nostrand and Company, 
Inc., New York, 1939), pp. 129-142. 
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laxation method,’ involving the solution of Laplace’s 
equation by first reducing it to a system of difference 
or algebraic equations which are solved by an arith- 
metical technique due to R. V. Southwell,? or the 
(c) conformal transformation method*-* may be used 


2R. V. Southwell, Relaxation Methods in Theoretical Physics 
(Oxford University Press, London, 1946). 
: ; rT Julia, Representation Conforme (Gauthier-Villars, Paris, 
931). 
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to derive the potential functions desired. Often, as in 
the present case, the last method is the most ap- 
propriate. 

The general expression for the potential function 
involving the particular geometry under consideration 
is derived and represented by expression (8). This 
expression is then evaluated for each of four boundary 
conditions, 


(a) between the grid and filament strands, 
(b) between the grid and through the filament strands, 
(c) between the filament and through the grid strands, 
(d) through the grid and filament strands, 


encountered with a typical thermionic vacuum tube. 
The resulting potential functions are illustrated in 
graphical form by Figs. 6 through 9. 

As illustrated by Fig. 1 the geometry under considera- 
tion comprises two spaced parallel equipotential planes 
1,1 between which are symmetrically interposed planes 
2,2 and 3 of isolated metallic strands. The strands in 
planes 2,2 are mutually parallel to those in plane 3 also, 
but strands in planes 2,2 are orthogonal to those in 
plane 3 as shown. Planes 1,1 are at a prescribed uni- 
potential value, E,. These planes are referred to as the 
anode or plate. Strands in planes 2,2 are also at a 
prescribed potential Eg with respect to potential Ec of 
the strands in assemblage 3. Planes 2,2 are referred to 
as the grids and plane 3 as the cathode or filaments. 

This structure is next represented on a composite 
plane commonly referred to as the W plane. This plane 
is comprised of two corresponding parts of the structure 
illustrated by Fig. 1 and will be referred to as the W, 
and W, planes. The former of these two planes pertains 
to the filament as center, as illustrated by Fig. 2a, 
whereas the other pertains to the grid as center, as 
illustrated by Fig. 2b. Finally, in order to derive the 
potential functions in the original structure, Fig. 1, it 
is expedient to transform conformally the composite W 
plane into a composite Z plane. The Z plane is neces- 
sarily composite and constitutes transformations of 
the two corresponding parts of Fig. 2. As illustrated by 
Fig. 3a and Fig. 3b these are referred to as ZX and YX 
transformations. The type of conformal transformation 
resorted to may be represented by the complex function 
W=InZ. From the theory of functions of a complex 
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Fic. 1. Structure under 
consideration, comprising 
metallic planes 1,1 with 
metallic strands in planes 

--4 2,2 orthogonal to those in 
plane 3. 
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*C. Caratheodory, “Conformal Representation,” Cambridge 
Tract 28 (Cambridge University Press, London, 1932). 

*W. G. Dow, Fundamentals of Engineering Electronics (John 
Wiley and Sons, Inc., New York, 1937), Chapter 2. 
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Fic. 2a. The W, plane with a strand in plane 3 of Fig. 1, 
e.g. the filament as center of coordinates. 
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Fic. 2b. The W:2 plane with a strand in plane 2 of Fig. 1, 
e.g. the grid as center of coordinates. 


variable it is shown that equipotential surfaces and 
field lines retain their original relationship by such 
transformation. Hence, the self and mutual capacitances 
and potentials between elements remain invariant as a 
result of such conformal transformations. With the aid 
of these transformations the potential function is de- 
rived and evaluated for each of the four important 
boundary conditions previously outlined. 

The physics of thermionic vacuum tubes involves an 
important well-known constant referred to as the am- 
plification factor, u. It represents a measure of the 
relative effectiveness of the control grid and pla‘ 
potentials in the determination of the electric fiewl 
strength at an infinitesimal distance from the filament 
or cathode surface. This constant is derivable with the 
aid of the potential function (8). Explicitly, the ampli- 
fication factor is expressed as a ratio of potential 
coefficients, (9) and (10), from relation (15). It is also 
expressible in terms of geometric factors, by expressions 
(13) or (14). Incidentally, the application of conformal 
transformation for the evaluation of the amplification 
factor was made early by King,® Vogdes and Elder,’ 
and others. 

The derived potential functions are applicable with 
a fair degree of approximation to structures of the type 
illustrated by Fig. 1 and under restrictions set forth in 
the analysis. Aside from the neglect of space charge 
and edge effects, the treatment is approximate to the 
extent that the derivation of exact electron trajectories 


®R. W. King, Phys. Rev. 15, 256 (1920). 
7 F. B. Vogdes and F. R. Elder, Phys. Rev. 24, 683 (1924). 
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is not strictly derivable with the aid of the present 
potential functions. An appropriate application of the 
theoretical results to a typical commercial thermionic 
vacuum tube terminates the work. 


Il. THE DERIVATION OF THE POTENTIAL FUNCTION 
BY CONFORMAL TRANSFORMATION, METHOD C 


In the light of the geometry revealed by Fig. 1, 
W, and W, half planes are constructed as shown in 
Figs. 2a and 2b. This follows as a consequence of the 


ent CENTER GRIO,G, 











Fic. 3a. The XZ conformal transformation with a strand in 
plane 3 of Figs. 1 and 2a, e.g. the filament as center. 


relation that the structure on either side of the YZ 
plane in Fig. 1 is identical. 

First, the W, plane is transformed into a ZX plane 
as shown in Fig. 3a. In the W; plane, the origin, i.e. 
the position of the filament center, Y, is represented by 
the coordinates 

¢' =2xncZ=0, 
and 
r= erncXc, 
in the ZX transformation. The size of the filament 
strand is represented by the coordinates 
¢ =0, 
ra? = ernc(—Re) | 
and 
r’ p = e2ncRc. 


“he position of the grid plane is represented by the 
coordinates 
¢'=0, 

and 

7 qr = ence, 
Note that, in the ZX transformation, all elements to 
the left of the filament C in the W, plane makes r’-0 
rapidly, but not negative. Hence, in the minute domain 
about the polar origin exists the left-hand grid and 
plate. Furthermore, no passage is permitted through the 
grid plane as this is assumed to be the immediate 
boundary by virtue of the fact that the grid strands are 
orthogonal to the filament strands. 

Next, the W, plane is transformed into a YX trans- 
formation, Fig. 3b. In the W2 plane the origin, i.e. the 
position of the grid center, is represented by the co- 
ordinates 

¢” = 2amngY = 0, 
and 
7’ = eeXG=1, 
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in the YX transformation. The size of the grid strand 
is represented by the coordinates 

¢”=0, 

rn ermaCRa) 

and 

9!” Br on emake 
The position of the plate is given by the coordinates 

¢”=0 
and 
9” on = Pena, 

Now, note that in the YX transformation, all ele- 
ments to the left of the grid strand, G, in the W2 plane 
makes r’”’-—0 rapidly. Hence, in the minute domain 
about the polar origin exists the filament, the left hand 
grid and plate. In other words, the complete W; plane 
or ZX transformation exists in the minute domain 
about the polar origin. 

Consider the ZX transformation again and select 
any point on the filament surface, such as A’ for which 
¢’=0, and any other point in space, such as T’. These 
are to be kept invariant in position during the formula- 
tion of the subsequent potential expressions. In this ZX 
transformation, the potential difference between T’ 
and A’ with ¢’=0, due to a positive whole charge,® 
+r’ at the grid center is represented by 


, 


r 
Er’ — Ea’ = —2(+17¢’) en (1) 
(Sc’— Rc’) 


Next, due to a negative whole charge, —7c’, at the 
filament kernel the potential difference is 


Er: — Ea = —2(—1¢’) 
[(r’)?+(Sx’)?—2r'Sx’ cos¢’ }} 
XIn 
Re’ — (Sc’— Sx’) 
where the numerator of the logarithm corresponds to 
T’P’ in Fig. 3a and the denominator thereof corre- 


sponds to A’ P’. Finally, due to a negative whole charge, 
—rtc’, at the filament center the potential difference is 


Er: — Eg = —2(—1¢’) 
[ (r’)?+ (S¢’)?— 2r'S¢’ cos¢’ }} 
XIn 
Ro’ 
where the numerator of the logarithm corresponds to 
T’Y in Fig. 3a and the denominator thereof corresponds 


to A’Y. As a consequence of the three whole charges 
the total potential difference is (observing that E4-= Ec’) 





» (2) 





» (3) 


r’ 


(Sc’—Re’) 
1 
Re’ { Re’ — (Sc’—Sx’)} 
x {L(r')?+-(Sx’)?—2r'Sx’ cose’) 
XL(r’)?+ (Se’)?—2r'Sc’ cos’ }}}. 


8 Line charges are referred to unless otherwise stated. 


Er —Ec= om 2(+ Tc’) In 


—2(—r¢’) In 
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Since, the path between the filament and plate planes 
is characterized by the value of ¢’ in Fig. 3a, it is 
appropriate, for illustrative purposes at least, to con- 
sider only two paths, i.e. centrally through and halfway 
between strands which correspond to ¢’=0 and 7, 
respectively. Furthermore, recalling that the charac- 
teristics of the composite Z plane entails only magni- 
tudes of distances, it follows that 


r’ 


Er — Ec= 2(- Tc’) in———_—_———_ 
(Se'— Rc’) 


—— Re'{ Re’ —(Sc’— Sx’)} 
+2(— 7’) In > 





|r’ +Sx’ | |r’ +Sc’ | 
Now, from potential theory (Sc’)?—1=(Rc’)? and 
Sx’=1/Sc’ so that, using half of the previous whole 
charges, the above expression assumes the form 
+ , (Ro’)? 
Er: — Ec=(—T¢’) In » (4) 
|Se’+1/r'||r’+Sc’| 
Next, the YX transformation, Fig. 3b will be analyzed. 
The potential difference between 7” and Sc with 
¢’=0, due to a positive whole charge, +7¢’’, at the 
grid kernel is represented by 
Ep —Ec=—2(+1¢") 
C(r’)2+ (Sx”’)?— 2v"Sn" cos¢”” }! 
X In . . (5) 
(Sx”’—Sc) 


The potential difference between T” and Sc, ¢”=0, 
due to a negative whole charge, —rp”, at the grid 
center is represented by 


Ere —Ec= —2(—rp”) 
C(r’’)2+ (S@’’)?—2r"Sg"" cos¢”” }} 
n . 
(Se —Se) 


The potential difference between J” and S¢ with 
¢’=0, due to half of the two whole aforementioned 
charges, is expressed as 











(6) 








Ep —Ec 
[(r’’)2+ (Sq@’’)?—2r’S¢" cos¢”’ }! 
=—(—rp”) In 
Sq” 
C(r’’)2+ (Sx")?— 29S" cos¢”” }? 
—(+7c¢”) In ° 
Sr” 


for case J, i.e., when Sc&Sx, Se. 

As previously discussed, paths for which ¢’=0 and x 
will only be considered. Hence, the above expression 
reduces to 


Ep — Eg= —(—rp") In| r”/S@"+1| 
—(+7¢”) In| Se”’r’’+1 | , (7) 


since Sx”’=1/Sg"’ as a consequence of potential theory. 
From the foregoing it should be noted that although 
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¢’ and r’ are derived from transformation ZX, whereas 
¢~” and r” are derived from transformation YX, the 
former variables are independent of the latter. Inci- 
dentally, when the chosen angles ¢’ and $” each exceed 
the secant angles of the transformed filament and grid 
circles, respectively, the coordinates are off the re- 
spective strands. 

Now, expressions (4) and (7) are next combined to 
obtain the over-all potential at any point, 7’ = T’’=T, i.e. 


(Re’)? 
n 
| Se’+1/r’| |r’ Sc’ | 
+rp" In| r"/Se"+1| —7¢" In| Se”r" +1], 
for 2Ec=0. Since 





Er=—te'l 


oo" { =ne(—70")} =a0'| =ne(—T0’)} 
(—71¢")=(nc/ne)(—Te’), 
the above expression for potential becomes 


(Rc’)? 





Er= (—re') ft +(nc/ne) 
|Se’+1/r'| |r’ +Sc’| 


X In| S6"r"+1| }+c+ tp’) In|r”’/Se”’+1|, (8) 


which represents the general potential function for the 
system of Fig. 1. It should be recalled that this ex- 
pression is valid only for case J, ie. when ScSg. It is 
necessary to return to the basic expressions (5) and (6) 
THIS REGION IS OUTSIDE OF THE 2x TRANS- 
FORMATION WHERE THE G PLANE IS THE 
BOUNORY 


GRID CENTER 
(r*91,¢'s o) 








PLATE PLANE 











Fic. 3b. The YX conformal transformation with a strand in 
plane 2 of Figs. 1 and 2b, e.g. the grid as center. 


in order to derive the general potential function for 
case IJ, i.e. when the aforementioned condition is not 
valid. 

Continuing the consideration of case J, particular 
solutions of (8) are required for the plate to cathode and 
grid to cathode potentials in terms of potential coeffi- 
cients and charges in the expression (8). These are 
obtained by first putting T on the plate, such that when 

¢’=0, 


/ 
7! = e2tnclots) 


and when ¢”=0, r’=e*"@ are substituted in ex- 
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Fic. 4. The actual W 
filament approximated by 
a parallel strand type, for 
\ case 7 of the analysis. 
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Fic. 5. Conventional static plate characteristic curves for a 
typical thermionic vacuum tube structure considered in the 
example. 


is readily obtained. Then, by putting T on the grid, 
such that when 

¢@=0, r’ =e, 
and when ¢”=0, r’=S¢’’"+R,” are substituted in ex- 
pression (8), the grid to filament potential, 


(Re’)? 





Eg—Ec= a Tc) In 
1 


S¢’- 





| e?rnce— So’ | 





e2™™ca 


nc 
+— In | Se” (Se"+Re"’) —_ 1 | 
NG 


+ ‘on tp’) In | +Re"/S¢" 





, (10) 
results. 


Relations (9) and (10) may be written in terms of the 
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potential coefficients, A1, Az, By, Bz, and the charges, 
tc’ and rp”, as 


Ep= A i(— tc )+Ao(+ we 


and (11) 


Ecg= By(—1¢')+ Bo(+7p’”) 
where the potential coefficients are 




















(Rc’)? 
A 1= In 
| 1 
| S;°/ —-—————— | e2™nc(atd) —S¢’| 
| e2™nc(atd) 
nc 
+— In| S¢@"e?*"@—1]. (12a) 
NG 
| e2*nab e2™nad 
A2= In| —1 in ’ (12b) 
| ud ff Sq” 
(Ro)? 
B,=In 
Se’— | e2enca_ S| 
e?*nca! 
+(nc/ng) In| Sq’"(Se""+Re")—1|, (12c) 
and 
By=I|n|+Re"/Se@’’|. (12d) 


Expression (8) can be solved with the aid of (11), from 
which the charges, rc’ and rp”, are derived in terms of 
specified values of Ep and Eg. 


Ill. THE AMPLIFICATION CONSTANT, u 


With an aspect to completeness it is now appropriate 
to formulate the relations for the amplification constant, 
u, for the geometric structure under consideration, 


Fig. 1. This constant is a measure of the grid and plate 


potentials in the determination of the electric field 
strength that exists just off the filament (or cathode) 
surface, in the absence of space charge. Hence, the 
amplification constant, u, can be expressed in terms of 
the potential coefficients, (12b) and (12d), that is 


Az In(e?*"4>/S¢"’) 
~~ InSg¢”/Re” ; 


Since S¢’”’=cosh2angRg and Rg’ =sinh2rn¢Re, 





u= (13) 


(2xngb—I\n cosh2rngRe) 
y= : (14) 
In coth2angRe 





as usual. 

This expression is also valid so long as the distance 
between the filament and grid is such as to justify that 
Sce«Sg¢ in the Z plane. When this condition is not 
fulfilled the coefficients in (13) must be derived di- 
rectly from expressions (5) and (6) as they stand. 

The amplification constant can also be represented in 
terms of the previously derived potential distribution 
function, that is 


w= —(Ep—Ec)/(Ec—Ec), 
at cut off. 
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IV. A PRACTICAL EXAMPLE 


To illustrate the foregoing theory and particularly to 
apply it, a thermionic vacuum tube is selected which 
approximately fulfills the essentials of the geometry 
portrayed by Fig. 1. The Western Electric type 104D 
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Fic. 6. Theoretically derived potential distribution between 
plates, along a path halfway between grid and filament strands 
of Fig. 1. 
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Fic. 7. Theoretically derived potential distribution between 
plates, along a path halfway between grid strands and centrally 
through filament strands of Fig. 1. 
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Fic. 8. Theoretically derived potential distribution between 
plates, along a path centrally through grid and filament strands 
of Fig. 1. 
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Fic. 9. Theoretically derived potential distribution between 
plates, along a path halfway between filament strands and 
centrally through the grid strands of Fig. 1. 
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thermionic vacuum tube is suitable for the immediate 
purposes. The pertinent geometric data obtained from 
the dismantled tube are: 


Grid strand radius, Re =0.00825 cm, 
Grid strand pitch, mg=2.53 strands/cm, 
Grid to plate distance, b=0.239 cm, 
Filament to grid distance, a=0.195 cm, 
Filament strand radius, Rc=0.0128 cm, 
Filament strand pitch, »c=2/cm. 


The filament wire pitch is derived from an equivalent 
model. The actual filament configuration in the chosen 
thermionic vacuum tube is an inverted W shape fila- 
ment, 13.0 cm long. To simplify the analysis it was 
considered valid, in the light of the assumption that 
ScKSg, to treat the actual filament geometry as sym- 
metrically disposed parallel strands orthogonal to the 
grid wires as shown in Fig. 4. 

The plate characteristics of the chosen thermionic 
vacuum tube are illustrated by Fig. 5. 

The potential function (8) is evaluated for this tube 
in the following manner: The values of the potential 
coefficients (12a) through (12d) in expressions (11) are 
computed by merely substituting the corresponding 
geometric constants and potential differences. Then, 
the charges rc’ and rp” in expressions (11) are solved 
for and substituted in the potential function, (8). From 
the foregoing geometric data it is evident that 


Sc’= cosh2rncRc= 1.013, 

So’ = cosh2rngRe¢= 1.008, 

Re’ =sinh2ancRc=0.1607, and 
Re” =sinh2rngRg=0.1304. 


The selected operating point with reference to Fig. 5 
is defined by 


Ep= 100 volts, 
Eg=—10 volts, and 
Ec=0. 


These are substituted in expressions (11) to yield the 
charges, 
(— rc’) = —.0138 statcoulomb 
and 
(+7rp”’)=+.0678 statcoulomb. 


Upon substitution of these charges in expression (8) 
with pertinent other factors the potential function, 


0.0258 





E= —0.0138} In- 


1| 
1.0134! |r’ +1.013) 
| r’| 


Ld 


r 





+0.790 In| 1.0087’’+1| | +0.0678 In 
. 1 





2) 


is obtained. 


This particular potential function is next evaluated, 


with the previously fixed values of Ep, Eg and Ec, for 
the following four important cases, between plates 
along a path (a) halfway between grid and filament 
strands, i.e. for ¢’=¢’’=7, (b) halfway between grid 
strands though centrally through the filament strands, 
i.e. for ¢’=0, ¢’’=7, (c) halfway between the filament 
strands though centrally through the grid strands, i.e. 
for ¢’=2, ¢’’=0, and (d) centrally through the grid 
and filament strands, i.e. for ¢’= ¢” =0. 

It should be recalled that the choice of sign in the 
absolute terms are dependent upon the chosen values 
of ¢’ and ¢”. As previously mentioned, only the two 
extreme values of these angles were selected. So that 
whenever these angles are equal to z, in the Z plane, 
the signs in the absolute terms are simply positive, 
whereas when zero the sign chosen is negative. 

The final results are illustrated by Fig. 6 through 9 
and are self-explanatory. 





Application of Thermodynamics to Chemical Problems Involving the Oxide Cathode 


Appison H. WHITE . 
Bell Telephone Laboratories, Murray Hill, New Jersey 


Received February 23, 1949 


Computations of the free energy of some chemical reactions involving dissociation or reduction of alkaline 
earth oxides have been made on the basis of rough estimates of the thermochemical quantities which have not 
been measured. The resulting values of free energy, equilibrium constant and equilibrium pressure of barium 
vapor derived for specified ideal conditions are tabulated. Several applications of these results are discussed. 


I. INTRODUCTION 


HE oxide-coated cathode is a complex system 
whose chemistry and physics are still the subject 

of extensive investigation.' Over the long history of re- 
search in this field there has been very little application 


1-J. P. Blewett, J. App. Phys. 17, 643 (1946). 
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of the powerful methods of thermodynamics to the 
chemistry of reactions involving these alkaline earth 
oxides.2 The primary reason for this is the lack of 
thermochemical data, particularly the high-temperature 
heat capacities of BaO and SrO. Nevertheless, it is 


2 For a review of such applications, see J. P. Blewett, J. App- 
Phys. 10, 831 (1939) ; see also reference (12). 
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possible to base illuminating calculations on reasonable 
estimates of such quantities. This paper is designed to 
provide a tabulation of the results of such calculations 
made several years ago which, it has since developed, 
are of interest to a number of workers in the field. It 
must be emphasized that the free energies of reaction 
thus derived are only approximations, but at a rough 
guess are correct to within five kilocalories per mole of 
BaO or SrO involved. 


Il. METHOD OF COMPUTATION 


The change of the Gibbs free energy accompanying a 
chemical reaction is 


AF°(T)=AH(T)—TAS(T), (1) 


where H and S are respectively the heat content and the 
entropy.* The superscript o indicates as usual that each 
substance involved in the reaction is conceived to be in 
its standard state. It is thus sufficient for the calculations 
below to know AH and AS for the reactions and tem- 
peratures in question. One convenient way to obtain 
these quantities is to start with the known values of 
AH(To) and AS(To) (To=298.1°K) in calculating the 
desired quantities at the higher temperatures of interest, 


according to 
T 


AC, dT; (2) 


AH(T)=AH(T).)+ 
To 
T 
AS(T)=AS(To)+ AC ,dlnT. (3) 
To 


AC, is the total molar heat capacity at constant pressure 
of the substances “produced” less that of those “‘con- 
sumed” in the reaction. The molar heat capacities C, of 
the substances involved below can all be expressed in the 
approximate analytical form* 


C,»=A+BT-—CT?+DT“ 

so that, for any of the reactions involved, 
AC,=AA+ABT—ACT~?+ADT-—. (4) 
Equations 1 to 4 yield, after the indicated integrations 


AF°(T) = AH(T»)—TAS(T») 
4+fsdA+frdB+fcAC+foAD (5) 
where 
fa=T—To—T In(T/T)), 
fe=TT.—-}(T°+T), 
fe=1/2T—(1/To)(A—T/2T»), 
fo=4T'—2TT,“'—2T,}. 


These functions of temperature are presented in Fig. 1 

with sufficient accuracy for calculations of this type. 
The thermochemical data used here are recorded in 

Table I, together with the sources from which they were 


*See e.g., G. N. Lewis and M. Randall, Thermodynamics 
(McGraw-Hill Book Company, Inc., New York, 1923). 
*See Table I, reference (d). 
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obtained. It will be noted below that some of the heat 
capacity data are employed beyond the temperature 
range of their validity as specified by the last column of 
the table. Where data are lacking, the basis of an esti- 
mate is briefly indicated in the footnotes appended to 
the table. In particular the high-temperature molar heat 
capacities of BaO and SrO are taken to be equal to that 
of CaO, which has been measured. This is empirically 
justified by the fact that at 298.1°K the three heat 
capacities fall within the range from 10.23 to 11.34 
calories per degree-mole, and that below room tempera- 
ture this range decreases with increasing temperature.‘ 

The values of free energy obtained as above for twelve 
chemical reactions at four elevated temperatures are 
recorded in Table II. The free.energy of formation of 
water in Row 1 is useful for obtaining the free energy for 
reactions not included here by addition or subtraction 
according to the standard procedure.’ Reaction 2 does 
not occur but is similarly useful in combination with 
other reactions. Reactions 3 to 12 are relevant to various 
aspects of the oxide cathode problem. 


Ill. THERMAL DISSOCIATION OF ALKALINE 
EARTH OXIDES 


In order to illustrate the application of the data to 
problems of interest the equilibrium constants of these 
reactions have been calculated according to 


AF°=—RT InK, 


and are recorded in Table III. One problem is to esti- 


LA 


6 





a 
» 

oO 
+ 
\ 








. 


Dall 








en 
SSS 





Nke 











wr 
-4 cP 
“ox 
-6 N 
2, 
2 
-2 ay 





-10 























wt ATT 


“12 
300 500 700 900 1100 1300 
TEMPERATURE IN DEGREES KELVIN 


Fic. 1. 





1500 


* Table I, reference (c). 





——— ers 


Caen a ee 


So 


" 
7 
Et 
. 
: 
i 


TABLE I. Thermochemical data employed in free energy calculations. 

















H(T»o)* S(To)e Cp =A +BT—CT-?+DT- (calories deg.-! mole-1). 
(Kilocalories/ Mole) (E. u /Mole) Ad BX10°4 CX10-54 De Valid range 

H: 0 31.23 5.797 1.003¢ 0 +14.63 273°-2500°K 
O: 0 49.03 9.629¢ 0.369¢ 0 — 42.83 300°-2500° 
H,0(g) — 57.798! 45.17 8.828¢ 1.917e 0 — 29.58 300°-2500° 
CHi(g) — 19.100 44.46 5.34 11.5 0 0 273°-1200° 
CO(g) — 26.428 47.32 7.74¢ 0.74¢ 0 —20 273°-2500° 
Ba(g) (+44.500)« 40.67 4.97 0 0 0 All 
BaO(s) — 133.100 16.8 10.0 4.845 1.0805 0 
Ni(s) 0 7.1 6.99 0.905 0 0 626°-1725° 
NiO(s) — 57.830 9.25: 13.693 0.833 2.915 0 273°-1373° 
C(gr) 0 1.36 2.673 2.617 1.169 0 273°-1373° 
Si(s) 0 4.50 5.74 0.617 1.01 0 273°-1174° 
Si0,.()) — 198.300 11.2 12.80 4.47 3.02 0 273°-1973° 
BaSiO;(s) — 356.040 28* 27.95! 2.056! 7.456! 0 
Mg(g) +35.907> 35.51 4.97 0 0 0 All 
Mg(s) — 145.760 6.4 10.86 1.197 2.087 0 273°-2073° 
Sr(g) (+39.2)™ 39.33 4.97 0 0 0 All 
SrO(s) — 140.700 13.0 10.05 4.845 1.080% 0 
Ca(g) +42.485> 37.00 4.97 0 0 0 All 
CaO(s) — 151.710 9.5 10.00 4.84 1.080 0 273°-1173° 

* International Critical Tables, Vol. V, p. 169 (1929), unless otherwise specified. 

>’ K. K. Kelley, U. S. Bureau of Mines Bulletin 383 (1935) 

eK. K. Kelley, U. S. Bureau of Mines Bulletin 394 (i936). unless otherwise specified. 

4K. K. Kelley, U. S. Bureau of Mines Bulletin 371 (1934), unless otherwise specified. 

* M. deK. Thompson, Total and Free Energies of Formation of the Oxides of Thirty-two Metals (Electrochemical Society, New York, 1942). 

'F. D. Rossini Research, Nat. Bur. Stand. 22, 407 (1939). 


® Heat of sublimation assumed independent of temperature and equal to Van Liempt’s value, 800°-1050°K 


Pay-Bas, 55, 468 (1936). 
» Assumed equal to values for CaO, 300° <T <1500°. 
i Seltz, DeWitt, and MacDonald, J. Am. Chem. Soc. 62, 3527 (1940). 


(J. A. M. Van Liempt, Rec. Trav. Chim. de 


3A. F. Kapustinski and K. A. Novosel'tsev, J. Phys. Chem. USSR 11, 61 (1938). 
k Assumed equal to the sum of S(T») for BaO and SiOs:, by analogy with the result observed in CaSiOs. 


! Assumed equal to values for Wollastonite, CaSiOs, 273° <T <1573°. 


™ Taking heat of vaporization of liquid strontium from (6) to be 37.046 kilocalories at 298°K 
1030°K of 2.190 kilocalories, assuming it to be equal to heat of fusion at 298°K. 


mate the pressure of gaseous Ba and Oz in equilibrium 
with solid BaO at 1273°K, as compared with its molecu- 
lar vapor pressure of 4X 10-7 mm.° For this we have toa 
good approximation® 


K= pps*po, (6) 


where the equilibrium partial pressures identified by the 
subscripts are expressed in atmospheres. For the ideal 
case with no source of Ba or O» except the dissociating 
BaO 


PBa= 2po,= (2K)}. 


This pressure as listed under #,, in the table is 2x 10-" 
mm of mercury, showing that at 1273°K the rate of loss 
of BaO by dissociation to gaseous products is negligible 
compared with the loss by evaporation of BaO mole- 
cules. In SrO and CaO however similar estimates yield 
“dissociation pressures” ,, in Table III which are 
smaller than the vapor pressures by a factor of less than 
ten. In practice this is of little significance when the 
pressure of oxygen due to contamination is higher than 
that due to dissociation of the alkaline earth oxide, in 
experiments such as the determination of the vapor 
pressure of this oxide. 

The equilibrium pressure of 8X10-'’ mm of Ba or 
4x 10~-” mm of O, obtained by extrapolating the data of 
Table III to 1000°K may seem in conflict with Dush- 


5A. Claassen and C. F. Veenemans, Zeits. f. Physik 80, 343 
(1933). 
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, and adding the observed heat of fusion of strontium at 


man’s® estimate of < 10~*° atmosphere for the dissocia- 
tion pressure of BaO at this temperature. However 
Dushman’s result presumably refers to the conventional 
meaning of the term “dissociation pressure;” i.e., the 
equilibrium pressure of oxygen above the oxide in the 
presence of the solid metal. This may be obtained from 
Eq. (6) by setting ps, equal to the vapor pressure of 
solid Ba, 1.2X10-? mm at 1000°K.’ The equilibrium 
pressure of oxygen thus derived is 2X10~* mm or 
3X10~-** atmosphere, in agreement with Dushman. 
Since no such pressures are attainable it is of course 
concluded that no solid barium nor for that matter any 
significant concentration of excess barium dissolved in 
BaO can be obtained by thermal dissociation of BaO at 
or near 1000°K. 


IV. PRODUCTION OF FREE BARIUM 


A widely accepted theory of the oxide cathode starts 
with the assumption that BaO is an excess electronic 
semiconductor by virtue of a slight stoichiometric excess 
of barium atoms dissolved in it.* The possibility of sup- 
plying such excess barium by the reaction of various 
reducing agents with BaO has been considered from this 
thermodynamic point of view. The results are shown, in 


6S. Dushman, Rev. Mod. Phys. 2, 433 (1930); vd. also J. A. 
Becker, Trans. Am. Electrochem. Soc. 59, 209 (1931) and Wagner, 
Electronics 1, 178 (1930). 

7 Table I, reference g. 

8 See e.g., J. P. Blewett, J. App. Phys. 10, 668 (1939) for a re 
view ; see also reference (10). 
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terms of the equilibrium pressure of barium vapor, by 
the data for Reactions 4 to 10 in Table III. The concen- 
tration of excess barium in the solid in equilibrium with 
this vapor can be estimated as to order of magnitude by 
assuming that Raoult’s law holds over the entire range 
of solid solubility. Then, since the saturated solution 
would be in equilibrium with solid or liquid barium, 


fes/fea°=ppa/PBa” (7) 


where fp, and fga° are, respectively, the equilibrium and 
the saturated mole fractions of dissolved barium, while 
psa’ is the vapor pressure of solid or liquid barium. The 
value of fp.° may be taken from Schriel to be of the 
order of 0.01 at 1400°K.® No estimates of fp, are in- 
cluded in Table ITI, since pp, is equally useful for com- 
paring the effectiveness of the various reducing agents. 
By way of example however, the pp, of 1 mm at 1473°K 
obtained from Reaction 5, together with pp,°=17 mm 
and fpa°=0.01, yields 6X 10~ for the concentration of 
dissolved barium in equilibrium with 10~° atmospheres 
of methane. 

Reactions 4 and 5 represent procedures that might be 
employed in a laboratory experiment designed to increase 
the excess barium content of BaO. The equilibrium 
pressures of barium vapor tabulated in the last column 
show that methane at a pressure of 10~° atmosphere is 
at least 100-fold more effective than hydrogen at one 
atmosphere in reducing BaO. The low pressure chosen 
for the methane calculation was that employed in an 
experiment which is described elsewhere.!° 

Reactions 6 to 10 are related to mechanisms by which 
BaO might be reduced in a commercial oxide cathode, 
where the (BaSr)O is in contact with a nickel core con- 
taining small concentrations of such impurities or 
additives as carbon, silicon, and magnesium. Reaction 6 
shows that nickel itself can have no effect, since the 
equilibrium pressure of barium vapor above BaO alone 
is higher than that above BaO in the presence of both 
nickel and NiO. In Reaction 7 carbon is shown to have a 
considerable effect, although it should be emphasized 
that Reactions 7 to 10 are written for saturated solid 
solutions of the reducing agent in nickel, whereas the 
commercial concentrations are probably much lower. 
Reactions 8 and 9 show that silicon would be much less 
effective than carbon if it were not for the side reaction 
between BaO and SiO, to form a silicate, here assumed 
for convenience in calculation to be BaSiO; rather than 
the experimentally observed Ba2SiO,." It will be noted 
that the equilibrium constants for Reaction 10 have 
been calculated for gaseous rather than liquid mag- 
nesium. However the barium pressures listed under pm 
are derived by taking the pressure of magnesium to be 
that in equilibrium with its liquid at each temperature. 





* M. Schriel, Zeits. f. Anorg. allgem. Chemie 231, 313 (1937). 

” Hannay, MacNair, and White, J. App. Phys. 20, 669 (1949) ; 
see also C. H. Prescott and J. Morrison, J. Am. Chem. Soc. 60, 
3047 (1938). 

"A. Fineman and A. Eisenstein, J. App. Phys. 17, 663 (1946) ; 
H. P. Rooksby, Nature 159, 609 (1947). 
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Taste II. Free energy changes accompanying certain reactions. 











AF®° (Kilocalories) 





Reaction T =873° 1073° 1273° 1473° 

1. 2H: +03 =2H:20(g) —95 —90 --84 -79 
2. 2CHs+02=2CO+4H2 —91 —110 -—128 —147 
3. 2BaO(s) =2Ba(g) +O: . 273 255 238 221 
4. 2SrO(s) =2Sr(g) +O2 273 254 236 218 
5. 2CaO(s) =2Ca(g) +O: 300 280 262 243 
6. BaO(s) +H: =Ba(g) +H20 89 82 77 71 
7. BaO(s)+CHs=Ba(g) +CO+2H:2 91 73 55 37 
8. BaO+Ni(s) = Ba(g) +NiO(s) 98 93 88 84 
9. BaO(s)+C(gr) = Ba(g) +CO 91 78 65 53 
10. 2BaO(s) +Si(s) =2Ba(g) +SiOa(/) 111 102 93 84 
11.. 3BaO(s) +Si(s) =2Ba(g) +BaSiO;(s) 87 77 69 60 
12. BaO+Mg(g) = Ba(g) +MgO(s) 1 2 4 6 








The results are thus to be compared with those of 
Reactions 6 to 9 and show magnesium to be thermo- 
dynamically the most effective of the reducing agents 
considered here. 

Tungsten is also of considerable interest as a reducing 
agent for alkaline earth oxides, but is not treated 
here in view of the careful work of Moore, Allison, and 
Morrison.!? 


Vv. CONCLUSION 


This is not intended to be an inclusive discussion of 
the application of thermodynamics to the chemistry of 


TaBLE III. Equilibrium constants of various chemical reactions. 











1 . 2 3 4 
Reaction T K Pm(mm) 
1. 2BaO(s) =2Ba(g) +O2 1073 1X10 6 X10715 
1273 2X10 2 x10-" 
1473 2X10" 1X10 
2. 2SrO(s) =2Sr(g) +Ox2 1073 2X10" 5 x<10715 
1273 3X10™ 3X10" 
1473 4X10 2x10 
3. 2CaO(s) =2Ca(g) +O2 1073 8xX10-% 9 x1071" 
1273 110-6 1 X1071% 
1473 8X10-" 9 X<10710 
4. BaO(s) +H2(g) = Ba(g) +H20(g) 1073 2X10" 3x10 
1273 7X10-“" 2x10 
1473 3X10™" 4x10 
5. BaO(s) +CHa(g) = Ba(g) +CO(g)+2H2(g) 1073 2X10" 6X10 
1273 4xX107 1 X1071* 
1473 4X10°* 1° 
6. BaO(s) +Ni(s) = Ba(g) +NiO(s) 1073 1107" 9 xX107" 
1273 7X10-*% 5X10" 
1473 4X10" 3 xX107% 
7. BaO(s) +C(graph) = Ba(g) +CO(g) 1073 1X10"* 8 x107% 
1273 6X10" 2X10 
1473 1X10 9X10"? 
8. 2BaO(s) +Si(s) =2Ba(g) +SiO2(l) 1073 2X10" 3X<107% 
1273 1X107* 8 xX107* 
1473 3X10-" 4X10 
9. 3BaO(s) +Si(s) =2Ba(g) + BaSiOs;(s) 1073 2X10" 1x10 
1273 2X107 1X10 
1473 1X10 3x10 
10. BaO(s) +Mg(g) = Ba(g) +MgO 1073 3X10! 7> 
1273 2X10" 4 X10> 
1473 1X10" 2X10*% 








® In equilibrium with CH, at 10~§ atmosphere. : 
>In equilibrium with magnesium gas at the vapor pressure of liquid 
magnesium. 


2 Moore, Allison, and Morrison, Phys. Rev. 75, abstract in 
rint; paper to be published; cf. also G. E. Moore, J. Chem. Phys. 
5 427 (1941). 
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the oxide cathode, but rather to illustrate the type of 
problems that can be attacked profitably even without 
complete thermochemical data. However the data of 
Table II may be combined to obtain the free energies of 
a number of reactions not explicitly considered above. 
For example, Reaction 1 or Reaction 2 may be added to 
Reaction 4 to obtain the free energy for reduction of 
SrO by hydrogen or methane respectively ; or Reactions 


3, 5, and 9 may be combined to yield the free energy of 
reduction. of CaO by graphitic carbon. 


Finally, it is to be emphasized that equilibrium condi- ' 


tions are seldom attained; hence the barium pressures 
listed represent upper limits not likely to be experi- 
mentally realized. 

The author isindebted to L. A. Wooten, C. Herring, and 
G. E. Moore for illuminating discussions of this subject. 





Time Constants for Vacuum Gage Systems 


S. A. SCHAAF AND R. R. Cyr 
University of California, Berkeley, California 
(Received January 20, 1949) 


A typical system for measuring pressures in a rarefied gas consists of a vacuum gage connected to the gas 
by a length of tubing, usually with an orifice of diameter smaller than that of the tubing at the end. The 
interior surface of this system will outgas at a fairly constant rate for a long period of time so that an equi- 
librium state will be reached with a definite non-constant pressure distribution along the tubing. The 
response time constant for the system to come to equilibrium and the difference in pressure between gas and 
gage are functions of the outgassing rate and the geometry of the system. An analysis is made assuming 
molecular flow, and formulas are obtained for the response time and the pressure difference. It is shown 
that there is an optimum value of the tubing diameter for a given orifice diameter, gage volume, and length 
of tubing. This optimum is rather critical; if the diameter is too small, response time constants of the order 
of hours may result. Experimental work has been done and the results compared with theory. This work 


was done under ONR contract. 


NOMENCLATURE 


A =inside surface area of gage. 
C=constant, }{1/(2*RT)#]~ 10-5 sec./cm at 20°C. 
L=length of tube. 

p= pressure. 
po= pressure in system (applied at time ‘=0). 
pi=initial pressure in gage. 
p2= equilibrium pressure in gage. 
Q=mass rate of flow. 

q= surface outgassing rate. 

q' =RTq. 

R=gas constant. 

ro= orifice radius. 

r,;= tube radius. 

opt = optimum tube radius. 

i= time. 
t,= time constant (lag coefficient). 
V =gage volume. 

x= distance along tube from gage. 
p=density. 

g=r/reL. 

O=ro/3L5/3/V, 


INTRODUCTION 


HE reading given by a vacuum gage connected to a 
vacuum system by a length of tubing (see Fig. 1) 

does not respond immediately to a change in the pres- 
sure level in the system. There is a certain time lag 
during which gas molecules flow in or out of the volume 
of the gage and tubing until an equilibrium condition 
is reached (see Fig. 2). The length of this time lag de- 
pends on the geometry of the gage and tubing system, 
and may be of the order of hours or more unless care is 
taken in choosing dimensions properly. There is, in fact, 
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an optimum size of tubing diameter for any given 
installation; too large a diameter increases the time 
lag because a larger volume must be brought into equi- 
librium, too small a diameter increases the time lag by 
increasing the resistance to flow through the tubing. 
In the present paper it is shown that the gage pressure 
varies exponentially with time, and a formula is ob- 
tained for the time constant or lag coefficient of the 
gage—i.e., the time within which the gage will have 
registered 1/e of its complete change following an 
instantaneous pressure change in the vacuum system. 
This formula is derived under the assumption of “free 
molecular flow,” i.e., that the pressure level is so low 
that all collisions between gas molecules may be neg- 
lected. The formula has been verified experimentally 
and shown-to be adequate for design purposes. 

The equilibrium state in general results in a pressure 
in the gage which is different from that in the system. 
This is apparently due to surface outgassing or in- 
gassing, i.e., the emission or adsorption of gas molecules 
at the inside surface of gage and tubing. The rate of 
the outgassing or ingassing also changes with time, 
resulting in a gage pressure change, but does so very 
slowly, being almost constant for hours or more. It is 
desirable to have an estimate of the initial equilibrium 
pressure difference between gage and system. Such an 
estimate is obtained. 


RESULTS 
The time constant is found to be 
tp=C(VL/r2+8V/3re+8ar2L/3re+4L?/nr,]. (1) 
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In dimensionless form this reduces to 





tire? 8r 4 

——=8/3+1/¢+ |= ones | (2) 

CV 3 rps 
where 0=79‘*L5?/V and g=r/re?L. Graphs of Eq. (2), 
which illustrate the rather critical optimum, are pre- 
sented in Fig. 3. The four terms appearing on the right 
side of Eq. (1) may be interpreted as the time constants 
associated with, respectively: (i) the flow of gas in the 
gage volume V through the resistance L/r; of the 
tube, (ii) the flow of the same gas through the resistance 
8/3ro? of the orifice, (iii) the flow of gas in the tube 
volume rr?Z through the resistance of the orifice, and 
(iv) the flow of gas in the tube through the distributed 
resistance 4L/x*r,3 of the tube itself. 

It is clear from (1) that ¢; may be decreased by de- 
creasing V and L and by increasing ro. In general, how- 
ever, there are physical limitations on how small V and 
L may be made, and on how large ro may be consistent 
with the purpose of the instrument. For any given ro, L, 
and V there is an optimum 7;. In most actual applica- 
tions the volume of the tube is small in comparison with 
that of the gage. The term 4L*/zr, in (1) may then be 
neglected and the optimum value of r; is determined by 


rope =[9re2V /169 }/5=[0.18r2V "5. (3) 


Experimental time constants have been obtained for 
three gage systems having the same values for ro and V, 
but with different values for r;, respectively smaller 
than, equal to, and larger than the optimum indicated 
in Eq. (3). The results are given under I, II, and III 
in Table I. 

The time constant for a gage system with no orifice 
other than the tube itself is found to be given by 


tp=CUVL/ré+4L2/ar;]. (4) 


Experimental values have been obtained for three such 
gage systems and are presented in Table I under IV, V, 
and VI. The theoretical and experimental pressure re- 
sponse as a function of time for the six gage systems is 
presented in Fig. 4 in dimensionless form, using the 
theoretical time constant appropriate to each gage 
system. 

The pressure difference at the initial equilibrium be- 
tween gage and system, due to outgassing flow, is 
found to be 


8A AL 16rxL xl’ 
Peree— Paystem=9'C —+—+ +— ° (5) 


3ro?>_—sar,® 310" r,? 





The outgassing rate g’ apparently depends upon the 
pressure and temperature history of the surface for 
several hours preceding a measurement, as well as 
upon the surface material. It may be positive (out- 
gassing) or negative (ingassing), and approaches zero 
with time if the surface is exposed to a constant pressure 
gas for several hours (see Fig. 2). From the measured 
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pressure differences and Eq. (5), we estimate that the 
maximum value of q’ is of the order of 10-? wcm/sec. 
for a Pyrex glass surface which has been exposed to 
pressures of the order of one micron for several hours. 
Further work on this phase of the problem is planned. 
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THEORETICAL DEVELOPMENT 


Under free molecular flow conditions, the flow rate at 
any point in a tube depends on the pressure and tem- 
perature distribution along the entire tube. For steady 
isothermal flow in a long tube with a constant pressure 
gradient, the flow rate is proportional to the over-all 
pressure gradient.' The basic simplifying assumption of 
this paper is to set the instantaneous flow rate propor- 
tional to the instantaneous local pressure gradient only ; 
the justification for this assumption is probably best 
found in the generally good agreement between theory 
and experiment indicated in Section 2. 

Under the above assumption, we have the axial mass 
rate of flow given by! 


Q=(—1r#/CRT)(dp/dx). (6) 


The radial mass rate of flow in from the sides due to 
surface outgassing is, for a ring of width “dx,” 


dQ=2nrqdx. (7) 


Using the gas law p=RpT, and assuming T and q are 
constant, a mass balance on the element of length dx 


| 


VERY SLOW STEADY CHANGE 
DUE TO CHANGE IN OUTGASSING RATE 
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Fic. 2. Response of pressure ‘gage. 


gives the partial differential equation 
Op/dt= (r./wC) (0? p/ Ax?) + 29'/r:. (8) 
The boundary condition at x=0, is obtained by writing 


1 Leonard Benedict Loeb, The Kinetic Theory of Gases (McGraw- 
Hill Book Company, Inc., New York, 1934), p. 293. 
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a mass balance for the gage volume. It is 
Op/dt=(r2/VC)(dp/dx)+Aq'/V atx=0. (9) 


For an infinitely thin orifice, the mass rate of flow is 
given by? 


Q= (4rCro?/3)(p— po). (10) 
The boundary condition at x= L is, hence, 
Op/ dt= (3r0?/8r?)(po— p) at x= L. (11) 


If the aperture consists in a short tube of length L’ and 
radius ro’, an effective orifice ro, given by 


3 L’\3 
ard /(1+- wert oF (12) 
8 ro 


should be used in Eq. (11). The case of no orifice other 
than the tube itself, i.e., p= po at x= L, may be obtained 
formally by letting re in Eq. (11). . 
The differential Eq. (8) is equivalent to a heat con- 
duction equation with a distributed heat source,’ and 
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Fic. 3. Variation of time constants with geometry of gage system. 
C=3/4(2rRT)$=2.61X 10-5 sec./in. at 530°R. 0=ro4/*L9/3/V, 


together with the boundary conditions (9) and (11), 
has the solution 


p= pot ¢'C[8A/3re+AL/r,3 


+ 167.9 L/3ro?+ 2L?/r,? | 
— ACq'x/rP—x2Cq'x?/r/?? 





Ea wr 
+> (A, cosd,x+B, sind,x) exp( - ="), (13) 
1 


2 See reference 1, p. 302. 
+See, for example, H. S. Carslaw and J. C. Jaeger, Conduc- 
tion of Heat in Solids (Oxford University Press, London, 1947). 
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where the eigenvalues \, are determined from the 
equation, 


3ar Pro? —8rAVX,? 
tandA,L= 7 (14) 
(3Vro?+8ar®)An 





From Eq. (13), we have at x«=0, 
p—p Ai*r 
~— ~exp( - *), (15) 
P2— pi rC 


where ), is the smallest eigenvalue. By definition, the 
time constant fz is given by 


tp=mC/APrt. (16) 





Solving Eq. (14) approximately for \1, we obtain finally, 
Eq. (1): 


trp~CLVL/r2+8V /3re+8ar?L/3re+4L*/ar; |. 


Eq. (4) is obtained by letting ro>— in Eq. (1). Eq. (5) 
is obtained directly from Eq. (13) by setting x=0 and 
letting tH. 


EXPERIMENTAL SET-UP 


The equipment used for making the experimental 
tests of time constants consisted of the small vacuum 
system in Fig. 5 which is ordinarily used for gage cali- 
bration. This system has a volume of approximately 
10 liters and a pumping speed of 25 liters per sec., 
yielding a time constant of 0.4 sec.—considerably faster 
than any of the gage systems tested. 

The gage systems (dimensions listed in Table I) were 
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Fic. 4. Experimental and theoretical vacuum gage responses. 
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TABLE I. Theoretical and experimental time constants 
for six gage systems. 











Tube I Tube II TubeIII TubeIV Tube V Tube VI 

V (cm) 164 164 164 52.5 52.5 52.5 
L (cm) 25.4 25.4 25.4 25.4 25.4 153 
rt (em) 0.11 0.41 1.6 0.13 0.05 0.41 
ro (cm) 0.02 0.02 0.02 — oo _— 
Time constant 

theory (sec.) 41 14 24 6.7 91 2.1 
Time constant 

exp. (sec.) 51 13 32 5.0 100 2.6 








constructed of Pyrex glass using standard commercially 
available ionization gages. Ionization gages were se- 
lected, since their pressure response is instantaneous in 
comparison with the Pirani and thermocouple gages 
which have a slow response due to the necessity of 
obtaining thermal equilibrium in the gage. The two 
ionization gages were driven by circuits which regulated 
the voltage on the grids and the emission current.‘ The 
ion currents were measured by a balanced vacuum tube 
voltmeter bridge circuit in which an Esterline Angus 
recording meter was used to measure the output. 

The measurements were made by pumping on the 
gage system until equilibrium had been reached, then a 
needle valve connected to the main vacuum system was 
opened, causing a sharp rise in the pressure. The varia- 
tion of pressure in the glass gage system was measured 


4L. N. Ridenour and C. W. Lampson, Rev. Sci. Inst. 8, 162 
(1937). 





Fic. 5. Experimental equipment. 


on the recorder as indicated in Fig. 2, and the data 
used in Figs. 3 and 4 were reduced from the recorder 
tape. 

This work was done under a contract between the 
University of California and the ONR. 





A Generalization of Nyquist’s Stability Criteria 


A. VAZSONYI 
Underwater Ordnance Department, Naval Ordnance Test Station, Pasadena, California 


(Received January 24, 1949) 5 


A generalization of Nyquist diagrams is developed. These new diagrams yield a lower limit to the values 
of “damping ratios” and so permit the evaluation of the stability and the “‘degree”’ of stability of systems 


characterized by ordinary linear differential equations. 


I. INTRODUCTION 


HE classical method to investigate ordinary linear 
differential equations consists of determining the 
roots of the characteristic equation and computing the 
response to some special input functions, such as a step 
or a pulse. Often one is interested only in the nature of 
the possible oscillations of the system. In this latter 
case from the roots of characteristic equation, the so- 
called “damping ratios” are computed. These dimen- 
sionless quantities can be considered as a measure of 
the oscillatory nature of the system. 
The classical procedure is very laborious for compli- 
cated systems. The frequency method developed by 
Nyquist* permits the stability investigation of very 





*H. Nyquist, Bell Sys. Tech. J. 11, 126 (1932). 
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complicated linear systems without an excessive amount 
of labor. From the Nyquist diagram of the system, it 
can be readily determined whether the system is stable 
or not. Furthermore, by special experience a certain 
insight about the “degree” of stability of the system 
can also be obtained. However, the classical “degree” 
of stability is not immediately correlated with the 
“degree” of stability of the system obtained from a 
Nyquist diagram. 

The purpose of this paper is to present a generaliza- 
tion of Nyquist’s theorem. It will be seen that the 
theorem in this more general form permits a quantita- 
tive correlation between the classical and the frequency 
(steady-state) method. More specifically, the general- 
ized Nyquist diagram yields a lower limit to the damp- 
ing ratios of the linear system. 
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Il. NYQUIST’S CRITERIA AS A MAPPING THEOREM 


Consider the general form of a homogeneous n’th- 
order linear differential equation with constant coeffi- 
cients 








d"0 d"—"@ d’6 dé 
Aw—+A,-1 +++++A2—+A;—+A)’= . (1) 

di” dy" dt? dt 

A, 

Zz Plane 
a 
+j 
-1 oO 41 
Az 








Fic. 1. A; and Az in complex plane. 


The auxiliary algebraic equation associated with this 
equation is 


A,N"+A,-1AN"!'+ eee +A od?+ AyA+Ap=0. (2) 


Let the roots of Eq. (2) be Au, As, ---, An. The general 
solution of the differential equation is 


0= ce" + coe*+- + + - +-6,€%*. (3) 


The differential equation is stable when each term of 
the general solution goes to zero with increasing time. 
Thus we have 

Lemma 1. The differential Eq. (1) is stable when all the 
roots of Eq. (2) have negative real parts. 

Let us consider now the complex function 


w=f(z)=Ap2"+Ani2™ +-°--+ArztAo. (4) 


This function maps the z plane into the w plane. In 
particular, assume that the positive half of the z plane 
(where the real part of z is non-negative) is mapped 
into the “interior” of the domain D of the w plane. 
‘According to our Lemma, Eq. (1) is stable when the 
domain D does not contain the origin of the w plane. 
The contour of the domain D is given by the mapping 
of the imaginary axis of the z plane and so is given by 


. f( jo) = An(jw)"+ An-i(jo)” "+ nica +A 1(jw)+Ao, (S) 


where w is a real number varying from —* to +. 
Thus the stability of Eq. (1) can be determined in 
three steps: 

1. Plot the f(jw) function in the w plane. 

2. Determine the “inside” of the domain limited 
by f(jw) (that is, determine on which side the positive 
half of the z plane is mapped). 
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3. If the domain contains the origin, the differential 
equation is unstable. If it does not contain the origin, 
the equation is stable. 


III. MEASURE OF STABILITY 


Consider the second-order differential equation 


d*6 dé 
A y+ 1i— + A o° = 0. (6) 
dt? dt 


This equation can be written in the form 


(d?0/dt)+-2{w(d0/dt)+w*=0, (7) 


INTERIOR 


om A we 
A3 ni wae iA ee 








Fic. 2. z plane. 


where 
w= (Ao/A2)! (8) 
and 
f=A 1/2(AoA 2)4. (9) 


w is the undamped natural frequency and ¢ is the dimen- 
sionless damping ratio. The quantity ¢ is a measure of 
the stability of the system. When ¢<0 the system is 
unstable, when ¢=0 the system is neutrally stable, 
when 0<¢<1 the system oscillates with decreasing 
amplitude, when ¢$1 the system is not capable of 
oscillations. Suppose now, that the system is oscilla- 
tory and let \; and 2 be the conjugate complex roots 
of the auxiliary equation of Eq. (6). By solving Eq. (7) 


M=[—f+j(1-—$) "Jo (10a) 
Aw=[—F—j(1— 5) * Jo. (10b) 


Figure 1 shows a diagram of \; and 22 in the complex 
plane. Let us call the angle between the vector \; and 
the imaginary axis the “‘excess’”’ angle u. Thus, 


w= Argr— 2/2. (11) 
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From Eqs. (10a) and (10b) follows, that 
sinu=¢. (12) 


_ Thus, the excess angle is a measure of the stability of 


the system. When y=0 the system is neutrally stable, 
when w= 7/2 the system is critically damped.'. 

Consider now the n’th-order equation (1). The roots 
of the auxiliary equation are real or conjugate complex. 
Each pair of conjugate complex roots corresponds to a 
second-order equation of the type of Eq. (6). Each pair 
of roots defines a damping ratio and an excess angle. 
Lemma 1 can now be generalized. 

Lemma 2. Assume that the roots of Eq. (2) in the upper 
half plane lie in the sector to the left of the excess angle u 
(Fig. 2). The damping ratios associated with Eq. (1) are 
larger than ¢=sinu. 


IV. GENERALIZATION OF NYQUIST’S CRITERIA 


Equation (4) maps the z plane into the w plane. 
Assume that the exterior of the sector of Lemma 2 in 
the z plane (see Fig. 2) is mapped into the domain D of 
the w plane. According to Lemma 2, if D does not con- 
tain the origin, the damping ratios of Eq. (1) are larger 
than sinu. Now, the upper half of the sector in the z 
plane can be defined by 


z= w(—sinu+j cosu)=wl—¢+j(1—¢7)#], (13) 


where w is a real number varying from 0 to +. The 
half of the contour of domain D in the w plane is 
given by 


flol—s+jA—s?)*]} = Anfol—f+j(i-—&)#}"} 
+Anafol—f+j(l—s?)*]}e +--+ +A. (14) 


Thus the following criteria can be set up: 


1. Plot the f{w[l—¢+ 7(1—¢?)*]} function for a given 
¢ in the w plane. 

2. Determine the “inside” of the domain limited by 
wl—¢+j(1—$?)#). 

3. If this domain does not contain the origin, the 
damping ratios of the system are all larger than 
the given ¢. 


V. NUMERICAL EXAMPLE 
Consider the torque amplifier 
I(d?0o/dt*) + B(d0,/dt)= M, (15) 


where J is the moment of inertia, 8 is the viscous damp- 
ing, 0 is the output, and M is the torque of the motor. 

Assume that the torque M is related to the error 
through the differential equation 


T(dM /dt)+M=Ke=K(0;—%), (16) 


where 6; is the input angle. Equation (16) expresses the 





‘ Equation (12) was communicated to the author by Frank J. 
ednarz. 
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fact that for a given error, ¢, the torque is developed 
through the time lag T. Write Eqs. (15) and (16) in 
operational notation 





p(Ip+8)0o=M (17) 
(Tp+1)M=K(0;—%). (18) 
By eliminating M, 
K K6; 
(19+ 89+ = . (19) 
Tp+1 Tp+1 


We are interested in the homogenous equation only. 
This can be written as 





K/B 
[1+ Joo. (20) 
p(I/Bp+1)(Tp+1) 


The performance of the system can be studied with the 


—7/¥ 
; 








Fic. 3a. Mapping of z plane. 


T 











Fic. 3b. w plane. 
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-1 
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¢-0 


Fic. 4. The transfer function G,(z) =1/z. 
aid of the function 
K/B 


+ . (21) 
2(1/Bz+1)(Tz+1) 





f(z)=1 


Actually, it is more practical and useful to study the 
transfer function 
K/B 
G(z)=—— , (22) 
2(I/B2+1)(Tz+1) 





In this case the ““—1” point of G(z) takes the place of 
the zero point of f(z). Assume the following numerical 
values: 


I= 2.00 slug-ft.? 
B=5.50 ft. lb./rad./sec. 
T=0.08 sec., 


and investigate K for the value giving damping ratios 
smaller than ¢=0.2. (Actually the system is of the 
third order and has only one second-order factor and 
so only one damping ratio.) Figure 3a shows the ex- 
terior of the z plane to be mapped and Fig. 3b shows 
the mapping in the w plane. For K, the value K=10 
ft. lb./rad. was assumed. It can be seen that for this 
value, the damping ratio is larger than 0.2 as the —1 
point is not contained in the domain. By changing the 
scale of the diagram, it can be shown that if K=30.2 


‘ft. Ib./rad., the domain will just contain the “—1” 


point, and so in this case the damping ratio ¢ is ex- 
actly 0.2. 
VI. BASIC TRANSFER FUNCTIONS 


It will be of interest to present plots of some of the 
basic transfer functions. The simplest transfer function 
is given by 


Gi(z)=1/2=1/ol—$+j(1—¢?)*]. (23) 


Figure 4 shows the transfer function for different values 
of ¢. When ¢=0 the standard transfer function is 
obtained. 
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Another important transfer function is given by 
G.(z)=1+7z=14+Tol—¢+j7(1-—¢)*). (24) 


This is a “phase lead’”’ transfer function. Figure 5 shows 
the plots for this function. 
A third important transfer function is given by 


G3(z)=1/1+Tz=1/14+To[—¢+ j(1—¢)#]. (25) 


This “phase lag” or “time lag” function is shown on 
Fig. 6. The transfer function loci are circular arcs 
through the “0” and “1” points. 

These three transfer functions provide the building 
blocks for transfer functions of more complicated 
systems. It can be seen that the usual theory of transfer 
functions can be easily generalized from the transfer 
functions introduced in this paper. 


w Plane 7 
G(Z)=1+TZ=L+wT(-€ + j vi- $2) 
6 4 2 ¢:0 





-1 0 7 
Fic. 5. The transfer function Go(z)=1+T7:. 
VII. GENERALIZED ATTENUATION PHASE ANALYSIS 


More recent methods? of servo-analysis plot the 
logarithm of the transfer function 


InG=In|G|+j ArgG=A+j¢. (26) 


Separate diagrams for the attenuation and for the 





Fic. 6. The transfer function 
1 1 
Gl) = (5TH “14 Tol tg 


2 See: James, Nichols, and Phillips, Theory of Servomechanism 
(McGraw-Hill Book Company, Inc., New York, 1947), p. 163. 
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argument are prepared, both plotted on semilog paper. 
The basic plot in this method is given by 


A=In|1+Tol—§+j(1—§*)*]]. (27) 


Figure 7 shows the attenuation curve of Eq. (27) for 
different values of ¢. It is important to note that the 
asymptotic characteristic of the attenuation diagram is 
independent of ¢. Figure 8 shows the corresponding 
phase diagram. Obviously, 


Tw(1—¢?)! 
1-—¢Tw 


tang= (28) 


The asymptotic characteristics of the phase diagram are 
of great importance. Obviously when Tw=0, g=0; 
when Tw=%, cosp=—f. Thus the value of ¢ for 
large Tw is affected by ¢ while the asymptotic character 
of the attenuation is independent of ¢. 


A Comment on Dr. Vazsonyi’s Paper, 
“A Generalization of Nyquist’s 
Stability Criteria” 


S. J. Mason 


Research Laboratory of Electronics, Massachusetts Institute of Technology, 
Cambridge, Massachusetts 


(May 20, 1949) 


HE Nyquist diagram of a function F(S), where S=r expj@, 
may be defined as the F plane map of the S plane curve 
|@|=2/2. Dr. Vazsonyi presents and discusses a generalization 
in which some fixed value of |@| between 2/2 and x is mapped. 
Such maps represent loci of constant Q.' If the appropriate 
critical point is not encircled by the F-plane locus, then the 
system not only is stable but the associated transient oscillation 
must exhibit a Q less than cot(2@—). Gain-phase-log frequency 
techniques, generalized to constant-Q loci by Dr. Vazsonyi, 
provide an approach to system design for a specified maximum 
allowable Q. 

In working with Nyquist Diagrams and gain-phase diagrams 
the writer has used a method of Q determination which may be 
of supplementary interest to readers of Dr. Vazsonyi’s paper. 
The method requires only the standard Nyquist plot of F(jw) and 
makes use of a graphical process of analytic continuation to 
locate the critical frequency in the left half of the § plane. The 
steps are as follows. The Nyquist diagram, F(jw) or logF (jw), 
is plotted. A uniform frequency scale is marked upon the curve 
in the region of closest approach to the critical point. A grid- 
work of curvilinear squares is then sketched in the neighborhood 
of the critical point, using the frequency marks on the Nyquist 
plot as starting points for the construction. The grid-work repre- 
sents the F plane map of rectangular coordinates in the S plane. 
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Te 


Fic. 7. Attenuation diagram for In|1+7:|. 





Fic. 8. Phase diagram for In| 1+T7z]|. 


The critical frequency S;= 01+ jw:, corresponding to the critical 
point, may then be read from the curvilinear coordinates of the 
sketch and Q; may be computed. 

If a logarithmic (rather than an arithmetic) frequency scale 
is marked upon the Nyquist plot, then the associated curvilinear 
squares represent the F plane map of the coordinates log|S| 
and @. A critical point which lies at @, in this grid-work indicates 
the value of Q,;=cot(20,—7). 

The graphical method described here is particularly simple in 
problems where the Q of the critical frequency is high, since the 
required curvilinear squares will then exhibit little distortion. 
Conversely, if the critical frequency lies far from the jw axis in 
the S plane and close to singularities of the function F(S), then 
elaborate sketching will be necessary in order to continue the 
function from the F(jw) plot to the critical point. 


1 The Q of a complex-frequency S: =¢1+ jw: is defined here as 
—les| on 
Q 2e1 


and the associated transient is of the form 


f(i)= exp(=Je""") coswit. 








The Propagation of Electromagnetic Waves in a Tube Containing a 
Coaxial D.C. Discharge* ** 


Pure Rosen*** 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 


(Received January 25, 1949) 


The boundary value problem involved in finding the propagation characteristics of electromagnetic waves 
in a cylindrical tube containing a coaxial d.c. discharge is solved. Ohm’s law is assumed to hold for the 
alternating current. This assumption has been verified theoretically for the case where the alternating 
electric field is small compared with the d.c. field. Curves which give the relationship between the complex 
propagation constant, y, and the complex dielectric coefficient, K.’, have been computed for the solution in 
which the TEM mode of the coaxial line is approached as K,’ becomes infinite. 





INTRODUCTION 


ECENTLY microwave techniques have been ap- 

plied to the study of d.c. gaseous discharges by 
D. E. Kerr'** of the M.I.T. Research Laboratory of 
Electronics. He used a cylindrical conducting tube 
which contained a d.c. coaxial mercury vapor dis- 
charge. He measured the propagation characteristics of 
microwaves as a function of d.c. discharge current by 
standing wave measurements. The details of this experi- 
ment will be given in C. 

The purpose of the present work is to solve the 
boundary value problem involved in this experiment, to 
find a relation between the dielectric properties of the 
discharge and the propagation of electromagnetic waves 
through the coaxial line. 

In the past, the propagation of the electromagnetic 
waves in a cylindrical tube containing two non-lossy 
coaxial dielectric materials has been investigated by 
Buchholz,‘ Pincherle,’ Frankel,’ and by Bruck and 





Fic. 1. Wave guide with lossy center conductor. 


° ~ ® This pay paper is part of a dissertation presented for the degree 
of Doctor of Philosophy in Yale University. 

“** This research has been assisted by the ONR under contract 
N7onr-288, Task Order IT. 

*** Now at Rensselaer Polytechnic Institute, Troy, New York. 

1 Kerr, Brown, and Kern, Phys. Rev. 71, 480 (1947). 

?—D. E. Kerr, Gas Discharge Seminar, M.I.T., Cambridge, 
Massachusetts, October 7, 1946. 

*D. E. Kerr and W. P. Allis, M.I.T., Cambridge, Massachusetts. 
Private communication. 

*H. Buchholz, Ann. d. Physik 43, 5, 313-368 (1943). 

5 L. Pincherle, Phys. Rev. 66, 118-130 (1944). 

*S. Frankel, J. App. Phys. 18, 650-655 (1947). 
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Wicker.’ These investigators only considered dielectric 
coefficients equal to or greater than unity. In the present 
work we deal with a lossy dielectric and furthermore the 
dielectric coefficient is permitted to be less than unity 
and even negative. A dielectric coefficient less than one 
is characteristic of media which contain a large density 
of electrons (see Appendix). 


A. BOUNDARY PROBLEM 


Our basic task is to solve the boundary problem of 
the wave guide shown in Fig. 1. For simplicity we have 
taken the outer conductor to be perfect. We also assume 
that the glass container for the ionized gas has a negli- 
gible effect. But before we can proceed we must consider 
a general set of Maxwell’s equations applicable to a 
medium containing charges and currents. Also, for sim- 
plicity we will assume our center conductor to be homo- 
geneous and isotropic. If we assume Ohm’s law to 
hold, i.e., 


J=<c.E, (1) 


where o, may be complex, o.= ¢r+i¢7, and also assume 
a time variation of the form e‘“‘, we can then use a 
set of Maxwell’s equations which contain no charges or 
currents but which is equivalent to that for a medium 
having a complex dielectric constant, K-,.* It can be 
shown that: K.=K—i(or/w), where K the effective 
dielectric constant, = e+ (¢7/w), where ¢ is the ordinary 
dielectric constant. 

Thus our problem reduces to that of a cylindrical 
wave guide with two dielectrics, which has been solved 
by Pincherle’ for real positive dielectric constants. For 
a circularly symmetric transverse magnetic wave we 
obtain: 


K. JIi(§i) 
Ci Jo(Fi) 
z. {Ji(Fo)/Jol (b/a) Fo ]} — {Ni(So)/N oL (b/a) fo J} 


"i { Jo(0)/JoL(b/a)Fo}} — { No(So)/ Nol (b/a)to]} 


( 54. G. Bruck and E. R. Wicker, J. App. Phys. 18, 766-769 
194 

® Leigh Page, Electrodynamics (D. Van Nostrand Company, 
Inc., New York, 1940), Chapter 6. 
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where: 
27 \? 


Ao/a 





pe—te= (K—0( (3) 


and 
So=aLly?+(2m/dro)*}!, Ss=aLy?+ Ke (2m/do)?}}. 


Here y=a+i§, is the complex propagation constant 
with a, the attenuation, and 6, the phase constant; 
K/= K;/Ko—i(o/wKo) = Kj/—i(a/wKo) (Ko is the di- 
electric constant of vacuum); Ao is the free space wave- 
length; 5, and a, are the outside and inside radii respec- 
tively; i refers to the inside or central region and 0 refers 
to the outside region or vacuum. 

In our work we will only consider a circularly sym- 
metric TM wave. This is necessary for simplicity since 


for the non-symmetric case a linear combination of TM * 


and TE waves must be used. From Eqs. (2) and (3), it 
can be seen that our problem is determined by a choice 
of the ratios Ao/a, and b/a. 


B. APPROXIMATIONS AND NUMERICAL SOLUTION 


Our problem is to solve Eq. (2) subject to Eq. (3)— 
an attempt to do this analytically is unsuccessful. 
Therefore, we are left with the methods of approxima- 
tion and of numerical solution. Both methods will be 
used so that one can be employed to check the other. 


Approximations 


If the complex dielectric coefficient, K.’, becomes very 
large, one would expect to obtain the solutions for a 
coaxial line. It can be shown’ that for a coaxial line: 


y=[nx/(b—a)—w?/c?}}, n=0, 1, Y *y 
or in terms of £0: 


So=nm/(b/a—1). (4) 


The TEM mode (m=0) therefore means that ¢o=0. 
Because of this result, it seems wise to try the approxi- 
mations to the Bessel functions for large ¢; and for 
large {o first, leaving the approximations for large ¢; 
and small ¢» for later. 

Large ¢;, Large {o.—Using the asymptotic expressions 
for the Bessel functions of the first and second kind!® 
in Eq. (2), and assuming that the imaginary part of ¢; 
is large we obtain: 








—(K.’)? 1 1 
7. a} 6 
o fol sin?{[(b/a)—1]fo} 








Ke’ Ju(Ss)_ 1 {Hi (o)/Ho[(6/a)f0]} — {Hi (S0)/HoL (b/a) fo} 


The imaginary part of ¢; would certainly be large for 
high o and high negative K;,’ since: 


fi=aly?+ Ki — (ie /wKo) )(24/do)*}*. 


Thus a high negative effective dielectric constant has 
the same effect as a high conductivity. Now as | K,’| >, 
[(b/a)—1]to=n, n=0, 1, 2 ---. However, {o>=0 does 
not satisfy our initial requirement for {» to be large. 
To make { as large as possible for the other solutions, 
b/a should be close to unity. 

Large ¢;, Small {o.—For small values of the arguments 
of the Bessel functions we have:!° 


(a) Jo(z)~1 (b) 1i@)~- 


(6) 





—2 2 —2 
(c) No(=)=— In( ) (d) Ni(z)=—. 
rg 1.782 mz 


Using these for fo, and the asymptotic expressions for 
i, in Eq. (2) we obtain: 


fo=73(2m4a/do)*LInb/a }-*(K.’)-, (7) 


provided that {o?<1 and | K.’|>>1. We see from Eq. (7) 
that {>—0 (TEM mode) as | K.’|—>~. The usefulness 
of Eq. (7) lies in guiding our numerical work. It is not 
accurate enough in itself to be used for calculation 
except if one is satisfied with error of 5 percent or more. 

The following important conclusions can be drawn 
from Eq. (7). First if K-’ is real and highly negative, 
fo lies along the positive imaginary axis near the origin 
and moves downward as K;,’ decreases. Secondly if K,’ 
is a highly negative imaginary, (|¢/wKoK;,’/|>1), the 
argument of {» is 67.5°. 


Dielectric Rod Type of Propagation 


From Eq. (7) it follows that for real and highly nega- 
tive K.’, fo lies along the positive imaginary axis near 
the origin and approaches the origin as K,’ decreases 
further. This suggests that for low negative and real 
K.’, {9 may be a positive imaginary but large-n absolute 
value. To verify this guess, we will take { to be a large 
positive imaginary. Then from (3) ¢; is also a large © 
positive imaginary. 

Now, the Bessel functions, J, and N,, can be ex- 
pressed in terms of the Hankel functions.'! If we make 
use of this on the right of Eq. (2) we obtain: 








ci Jolfs) Fo {Ho (So)/Ho™(b/a) to ]} — {Ho (So) /Ho((b/a) fo ]} , 


From the asymptotic expressions for the Hankel func- 


*Sarbacher and Edson, Hyper and Ultra High Frequency Eng. 
(John Wiley & Sons, Inc., New York, 1943), Chapter 7. 

Stratton, Electromagnetic Theory (McGraw-Hill Book Com- 
pany, Inc., New York, 1941), p. 357ff. 
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(8) 





tions, it can be seen that the ratios involving the H, 
functions become much smaller than the ratios involv- 


1 Jahnke and Emde, Tables of Functions (Dover Publications, 
New York, 1943), p. 134. 
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ing the H,°” function as {> becomes a large imaginary. 
We therefore have: 


Ki Ii(5)_ 1 Hi) 
fs Jolt) fo Ho (to) 


Equation (9) is identical with the equation obtained by 
Hondros and Debye" for dielectric guides. This means 
that for large positive imaginary fo, the fields attenuate 
rapidly as one goes radially outward from the center 
into the vacuum, i.e., most of the energy is confined in 
the central dielectric rod. This is just the reverse of the 
case where {>=0 (TEM mode) where the energy is 
confined to the annular space. For large £» and large ¢; 
Eq. (9) can be put into the form: 


(2ma/Xo)? 
Sher eee ° ‘ (10) 
K/+1 








(9) 


This relation could have been obtained from Eq. (5), 


since: 
b 
sin| (~)- 1 fo} =, as {fp i@. 
a 


From Eq: (10) we conclude that {>=i« for K.’=—1 
and as K,’ becomes more negative, {o travels down the 
positive imaginary axis toward the origin. For a small 
enough imaginary part of K.’ (> nearly a positive 
imaginary), {> does not become infinite at K.’=—1, 
and we have a resonance effect with the function rising 
steeper the smaller the loss. 


1S 4) 
acc 


0.9 


0.8 








0.7 -o 3 i 


. 


Fic. 2. {> and aa plotted as functions of K;’ for the 
lossless case when fp is real. 


“” Hondros and Debye, Ann. d. Physik 32, 465 (1910). 
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We have not examined what occurs for K,;’>—1. But 
we will discuss this below. 


Numerical Solution 
We now put Eq. (2) into the form: 
Ke Ii(§i)Fox(b/a, £0) =F: Jo(FiW(b/a, fo), (11) 


where 
x(b/a, £0)=Jo(Fo) Nol (b/a)fo] 
—Jol(b/a)foJNo(So), (12a) 
and 
¥(b/a, £0) =Ji(Fo) Nol (6/a) fo] 
—Jof(b/a)FoJNi(So). (12b) 


The auxiliary Eq. (3) remains the same. 

The values of the ratios b/a, and \»/a were taken to 
be 3, and 52, respectively. These are the values used 
by Kerr,’ in his experiment. Thus Eq. (3) becomes: 


0.16(K.’—1) =¢2—0?; (13) 
also 
fo= (a*y?+0.16)*, :=(a*y?+0.16K.’)'. (14) 


Because x and wy are functions of {> only for constant 
b/a, it seems wise to use a trial and error method with 
fo as our independent variable. Also, since £o is limited 
to the first quadrant as is easily seen, our task is thereby 
lessened. Therefore x and y were computed as func- 
tions'* 4 of |f| for constant arguments of fo. The 


08 


0.6 ag 
\ 


a a a 3.6 7 6-3 ON 12 
. 








Fic. 3. {> and af plotted as functions of K,’ for the 
lossless case when fo is imaginary. 


13 Table of the Bessel Functions Jo(z) and J,(z) for Complex 
Arguments, Math. Tables Proj. (Columbia University Press, New 
York, 1943). 

4 Table of Bessel Functions of the Second Kind Yo(z) and Y\(s) 
for Complex Arguments, Math. Tables Proj. (microfilm copy). 
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Fic. 4. (a) p2—p: plotted against trial values of | {o| for various 
values of y, the argument of fo; (b) @2—86, plotted against trial 
values of | {| for various values of ¢; (c) the curves of p2—pi:=0 
and 6:.—6,=0 plotted in the {o plane, the intersection of which is 
the solution. 


arguments of {9 were chosen in the range 0°-90°, and 
were spaced at 5° intervals. As for the absolute values, 
|fo|, the range was chosen from | ¢»| =0 to | ¢o| = 1.60 
and intervals of 0.10 were used. 
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Fic. 5. The argument of {> is plotted as a function of ¢/(wKo) 
holding K;’ constant for various values of K;,’. 
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The reason for choosing this range of |{o| is as fol- 
lows. In our work we shall be interested only in the 
solutions of Eqs. (11) and (13) yielding the lowest at- 
tenuation. If no discharge is taking place in the central 
region, the dielectric coefficient of the center is equal to 
unity and we have essentially a cylindrical wave guide. 
The mode for a cylindrical wave guide which is both 
circularly symmetric and is attenuated least is the 
TMo: mode for which {9>=2.405/3~0.8.° If one starts 
with this mode and changes K,’ continuously, it is 
possible to follow the solutions and they should be of 
lowest attenuation. If now the center conductor has 
| K.’| =, the mode of the cylindrical guide will either 
change to the TMo; mode of the coaxial line with 
to= 2/2, or become the TEM mode ({)=0). It seems 


_more probable that a cylindrical guide excited with the 


TMo: mode will tend to go into the mode having a field 
configuration most nearly like it, if a center conductor 
is added with | K.’| = ©. In as much as the TMo; mode 
of the cylindrical guide has the same field configuration 
as the TEM mode of the “coax” (except for the axial 
electric field) we expect the TEM mode of the “coax” 
as |K,.’|—>+. But our work is definitely incomplete in 
its restriction to this transition. 


Real K.’ 


We first consider the conductivity equal to zero. Since 
y can then either be a pure positive imaginary or a pure 
positive real if c=0, {> must be either a pure positive 
imaginary or a pure positive real. 

Let us consider fo, real. For a particular K,’ we choose 
a trial {> and substitute this {> into Eq. (13) to find the 
corresponding {;. With this particular pair, fo, {; we 
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Fic. 6. The absolute value of {> is shown as a function of K,’ 
with o/(wKo) taking on various constant values. 
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Fic. 7. The absolute value of {o as a function of K;,’ 
when @/(wKo) is small. 


find both the left and right hand sides of Eq. (11), and 
when the difference between these is plotted against 
the trial {), we obtain a zero difference for the correct 
value of {o. In this manner Fig. 2 is obtained. Values of 
Kj from 0.5 to —11.5 inclusive with intervals of 1 were 
used here. Figure 2 shows that when no discharge is 
taking place in the center region (K;,’=1), | o|~0.8 

(TMo: mode of cylindrical guide), and when K;,’ be- 
comes extremely negative, |{o|—>/2 (TMo: mode of 
coaxial line). Also plotted in Fig. 2 is the product aa 
(a is the real part of y), obtained from the {o curve. 
Because of the values used here for b/a and \»o/a these 
solutions yield completely attenuated waves. 


Imaginary {» 


With {> imaginary and with the same method as 
above, Fig. 3 is obtained. This figure verifies what was 
found from Eqs. (7) and (10). Namely, for K,/=—1, 
fo=ie and that ¢,—0 (TEM mode) as K,/—>— ~. Also 
for K;,’ slightly below —1, the propagation resembles 
that in a dielectric rod. One other important result can 
be obtained from Fig. 3; it is that above K,/=—1, 
(—1<K,/ <1) no solution can be found for imaginary fo. 

Also shown in Fig. 3 is the product a8 plotted against 
K;’. Solutions with {9 a positive imaginary can yield 
only completely unattenuated waves. 

It was stated previously that the TMo; mode of the 
cylindrical guide can either go into the TEM mode of 

_ the “coax” or into the TMo; mode of the “coax” 
as |K.’|-+«. If the conductivity is zero and one 
started out with a cylindrical guide and makes K,’ pro- 
gressively less than 1, continuing into negative values, 
we can either arrive at the TMo: mode of the “coax” 
({o= 2/2) as shown in Fig. 2 or we can make a discon- 
tinuous jump in {) at K,/=—1 and then continue 
downward along the positive imaginary axis to the 
TEM mode of the “‘coax” ({)»=0) as shown in Fig. 3. 
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The latter alternative means that {o in Fig. 3 is really 
a continuation of the part of {> in Fig. 2 above K,’= —1. 
This latter alternative, as will be seen below, gives rise 
to continuous curves when the conductivity is not zero. 
It is this case that will be investigated, that is, the case 
where the TEM mode of the “coax” is approached 
as | K.'|-@. 


Complex K.’ 


We are now ready to deal with the case where the 
dielectric coefficient is complex. The method used for 
real or imaginary fo is not useful here since the differ- 
ence between the right- and left-hand sides of Eq. (11) 
is a complex number and cannot be plotted readily 
against the complex trial ¢». However, if the left-hand 
side of Eq. (11) is put into the form p,e and the right- 
hand side into p2e“, the following method can be used. 

















{ 3 3 At =) aT 13 


Fic. 8. The argument of {> as a function of K;’ when 
a/(wKe) is held constant. 


Plot p2—p; as a function of |{o| for various constant 
arguments of {o(¢1, ¢2, g3) and also plot @2.—4@; in a 
similar manner (see Fig. 4(a), (b)). Then in Fig. 4(a), 
points A, B, C are points where p2—pi=0, and in Fig. 
4(b) points D, E, F are points where @.—6,=0. When the 
curves of p2s—pi=0 (determined by A, B, C) and of 
62.—6,=0 (determined by D, E, F) are plotted on the 
fo plane (Fig. 4(c)) we get an intersection where pi= 2 
and where 6,= 62. This is our solution, that is the value 
of {> which satisfies both Eqs. (11) and (13) for a 
particular K.’. 

The values of the imaginary part of K.’ were taken in 
the range 0 to —24 inclusive with steps of 2 and the 
range of the real part of K.’ was the same as before. 
Thus —11.5<K/ $0.5, 0€ (¢/wKo) £24. 
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At the start it was found best to begin with the 
imaginary part of K,’ small, so that the solution would 
lie not far off the real or imaginary axis, as the case may 
be. Then we progressively decrease the imaginary part 
of K.’ and follow our solution along. Thus we could 
narrow down the region in the {» plane where our solu- 
tion occurred. We could then use 2 arguments of £o 
spaced 5° apart for our trial solution and 2 straight lines 
for the curves of pe—p,;=0 and 6.—6,=0. In this man- 
ner Figs. 5, 6, 7 were obtained. 

Figure 5 which is a plot of the argfo against ¢/(wKo) 
for constant K,’ shows that as the imaginary part of 
K,' decreases the argument {> approaches 67.5°. This 
was found previously from Eq. (7). It is also seen from 
Fig. 5 that curves of constant K,’ for K;’ 2—1 start 
from 0° and for K,/<—1 they start from 90°. At 
K,/=—1 the curve goes down to 0° and then passes 
discontinuously up along the axis of ordinates to 90°. 
This is better seen from Fig. 8 which is a cross plot of 
Fig. 5. In this figure we note the discontinuous jump 
at K/=—1 for ¢/(wK»o)=0. However, if (¢/wKo)+0, 
we obtain a continuous transition. This was mentioned 
previously when it was stated that the TMo: mode of 
the cylindrical guide (K,’=1) will either become the 
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Fic. 9. The product a@ is plotted against K,;’ when o/(wKo) 

takes on various constant values. 


= 


TEM mode or the TMo; mode of the coaxial line as K;’ 
approaches — , 

Figures 6 and 7 which are plots of | {o| against K,’ for 
constant ¢/(wKo) show resonance type curves. That is, 
the peaks are steeper and narrower as the losses become 
smaller. Also the curves in Fig. 8 resemble the phase 
angle curves obtained in radio resonance measurements. 
This resonance behavior was hinted at in Eq. (10). 
Figure 6 also shows us that | {o| approaches zero (TEM 
mode) as |K,.’| goes to . This is just the solution 
we wanted. 
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The solution in which the TMo: mode of the “‘coax”’ 
is approached as | K,’|—> has not been investigated. 
This solution will yield waves of greater attenuation. 
Since no anomalies were found in the results here ob- 
tained, this solution either does not intersect with the 
solutions investigated or it is tangential at regions where 
the solutions become one and the same. 

Figure 7 shows | {| plotted against K,’ for low values 
of o/(wKo). The part of the curve for ¢/(wK»o)=0 which 
approaches {9=2/2 (TMo mode of “coax”) is to be 
ignored for the solution of interest in the region K// < —1 
because of what was said previously. The extreme steep- 
ness of the curve for o/(wK»)=1 indicates that nu- 
merical work in the low conductivity region is not as 
precise as for the higher conductivity regions. However, 
an over-all probable error in both |f{o| and argo is 
estimated to be about 1 percent and the maximum error 
in each is about 2 percent. 

Figures 9-13 show plots of aa and a8 which were 
obtained from the curves of {o. In Fig. 11 the curve for 
a/(wKo)=0, K,/ 2 —1 lies along the axis of the abscissae, 
ie., a8=0. In Fig. 12 the curve for ¢/(wKo)=0 for the 
solution of interest is the one indicated by the arrows. 
This curve has a discontinuity at K,/=—1 and then 
lies along the axis of the abscissae (aa=0) for K,/ <—1. 
The other solution for ¢/(wK»o)=0 yields a curve that 
approaches da=1.52 at K/=—. 

In Fig. 13, we have drawn a map of ay as a function 
of K,. That is, curves of constant K,’ and constant 
a/(wKo) which intersect are plotted in the ay plane. 
(For some unknown reason) these curves seem to be 
nearly orthogonal to each other. This figure should be 
most useful for obtaining information at a glance. 

We should mention that for K.’=1 (cylindrical wave 
guide) we could have any of the higher circularly sym- 
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Fic. 10. The product aa is plotted against K,’ when o/(wKo) 
takes on various constant values. 
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metric modes. That is, we could have started with any 
of the TMo, modes of the cylindrical guide at K.’=1. 
We have chosen to start with the TMo; mode because it 
has the lowest attenuation. To investigate these other 
solutions would involve tremendous labor and is at 
present not feasible. Besides, tables of the Bessel func- 
tions for complex arguments are not available for abso- 
lute values of the arguments of the Bessel functions 
beyond 10. 


C. COMPARISON WITH EXPERIMENT 


Kerr obtained standing waves by placing a shorting 
bar at one end of the wave guide and used a traveling 
E-probe and a standing wave detector to measure the 
dependence of the radial electric field E, on the axial 
distance from the short. He intended to measure the 
attenuation, a and the phase constant, 8, by this 
method. If the short is placed at Z=L, the theoretical 
dependence of EZ, on the distance from the shorted 
end is:** 


E,~([cosh(2ax)—cos(28x) }!, x=Z—L.t 


By comparing this with the experimental dependence 
of E, on x, Kerr was able to determine a and 8. By 
varying the direct current in the arc he obtained in- 
formation on how a and £ vary with arc current. 


a 
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i 
Fic. 11. The product a8 as a function of K;’ when ¢/(wKo) is small. 


It is the purpose of this work to find y as a function 
of the complex dielectric coefficient, K.’, of the dis- 
charge. Thus if a and 8 are known from experiment the 
value of K.’ would be known. Therefore, knowing K-’, 
it may be possible to relate certain parameters of the 


t This is seen as follows. Since we now have waves traveling in 
both directions along the axis of the cylinder we have: 


E,~(Bye~™” — Bye") f(r). 


At Z=L, E,=0 for all r and therefore B,/B,=e~?7". If we choose 
B,=4, we have: E,~sinh(yx). Now since E-probe measures 
(E,E,*)* we have (E,)x dependent part ~(cosh2ax—cos2fx)}. 
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Fic. 12. The product ga as a function of K;’ when ¢/(wKo) is small. 


discharge, such as electron mean free path, collision 
cross sections and the number of electrons per cc, with 
its measurable properties (see Appendix). 

Available data are not as good as is desirable. Diffi- 
culties arise from leakage of the field through a hole 
made by the arc and because of the fluctuation of certain 
discharge parameters. The latter reason is most impor- 
tant since, as will be seen, the data indicate ¢/(wKo) is 
small and good data are difficult to obtain in a region 
where the curves of Figs. 6 and 7 are steep. However, a 
definite trend is observed. 

For zero arc current (K,/=1), {00.7 and as the 
current is increased {») remains substantially real, in- 
creasing to approximately unity at a current near 1 
ampere; {9 may have a very small imaginary com- 
ponent, but this is very difficult to measure. Upon in- 
creasing the arc current further, a discontinuity occurs 
with {») becoming nearly a pure imaginary (the argfo 
remaining at 86°) and | {o| decreases from about 1.3 at 
1 ampere to about 0.6 at 1.4 ampere. 

In terms of a and 8 this can be described as follows. 
At zero current @ is zero and the wave is attenuated. As 
the current increases 6 remains either zero or is very 
small while a increases rapidly. Then we have (at the 
discontinuity) a sharp drop in a, while 8 becomes high 
suddenly. Finally, a becomes smaller still while 6 drops 
and approaches the free space value (TEM mode) at 
the highest current. 

The facts can be explained by our results. Initially at 
zero arc current we have cylindrical guide with ¢o~0.8. 
With the discharge present, at low currents, we have a 
low pressure discharge through which high frequency 
waves are passing. Therefore, we should expect K.’ to 
be real or nearly so (see Appendix). We can therefore 
explain the low current data with the part of Fig. 2 
for K/>—1, i.e., {o rises to approximately unity as K/’ 


JOURNAL OF APPLIED PHYSICS 








—_ hte on wa fh OUftlClU CUCU 


> a tee tet 


‘2s Pwcen 2a 








SS a Ve ee CL” 


nw eos & - 








decreases. The decrease in K,’ is due to the fact that 
the electron density is increasing as the arc current goes 
up (see Appendix). 

The discontinuity between the low and high current 
regions indicates the conductivity is very low. This is so 
because as was seen previously, a violent jump in the 
argto from 0° to about 90° takes place at K,/= —1 when 
the conductivity is low. The fact that ¢/(wKo) is small 
and that K,’~—1 is seen from the fact that just after 
the discontinuity | |=1.3 and the argt»>=86° (see 
Figs. 5, 7, 8). Then, finally as the current increases 
further | {o| is decreasing and approaching f>=0 (TEM 
mode) since K,’ is going to high negative values. This is 
seen from the data at 1.4 amp. (| {o| =0.6, argéo= 86°), 
which means K,;’~—8 and a/(wKo)~2 (see Figs. 5 
and 6). 

The reason the conductivity, which was previously 
infinitesimal, has increased is because the pressure of 
the arc and the d.c. field maintaining it have increased 
with arc current (see Appendix). 

Professor Allis* of M.I.T. has estimated the data at 
1.4 ampere (|{o|=0.6, argfo=86°) to mean that 
Kj’ = —8, o/(4Ko)~2. He used a power series approxi- 
mation in the high current region. These results check 
ours semi-quantitatively. 

One sees many practical aspects of the wave guide 
with a coaxial d.c. discharge. By controlling the dis- 
charge current one can vary the impedance of wave 
guide continuously. Also because of the fact that at a 
particular arc current, a rapid transition occurs between 
a completely attenuated wave and a propagated wave 
with small attenuation whose phase constant is a con- 
tinuous function of arc current, this wave guide can be 
used as a radio frequency switch and as a modulator. 


APPENDIX 


Above we had assumed that Ohm’s law held. Now 
we wish to justify this assumption theoretically using 
the methods of kinetic theory. We will show that 
Ja.c.=o-Ea.c., where o- is a function of the d.c. field. 

Consider the Boltzmann transfer equation :'*—"” 


df/dt+v-Vi+a-Vf=6f/at\|., (15) 


where f is the distribution function of electron velo- 
cities= f(v, r, ¢); v is the velocity, a is the acceleration, 
r is the position vector in configuration space, V and V 
denote the gradients in configuration and velocity space 
respectively, and 5f/é/|, denotes the rate of change in 
f due to encounters. 

Let 


f=fO+(1/r)v-£, (16) 


This is the start of an expansion of f in spherical har- 
monics;!® those spherical harmonics higher than the 


*S. Chapman and T. G. Cowling, The Mathematical Theory of 
Non-Uniform Gases (The Cambridge University Press, New York, 
1939), Chapter 3. 
®H. Margenau, Phys. Rev. 73, 297-308 (1948). 
 H. Margenau, Phys. Rev. 69, 508-513 (1946) 
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second are neglected. Then if Eq. (16) is substituted 
into Eq. (15) provided that a= (eE/m) where e, m are 
the electron charge and mass, respectively, we obtain: 








af y afm vf, cE 
$v Vf -4y-0/ )+— vy 
v ot m 


at . 
eE VAM) Sf) 8 sy-£M 
+—y/ )- +—( ) - (17) 
m v dt |. dX 9 s 


If we only consider encounters with molecules, then it 
can be shown that :16 18 


5(v-£/v) 
ét 














yf 
Te ’ (18) 


e v 








where / is the gas kinetic mean free path. Upon using 
Eq. (18), Eq. (17) yields two equations; one contains 
terms linear in v,, vy, v, and the other contains constant 
terms and terms containing v,”, v,7, v,7; vz0y, V2, etc. 











"3 4 oe 7 > Ke) 


Fic. 13. Curves of constant K;,’ and constant o/(wKo) 
plotted in the ay plane. 


The equation linear in vz, v,, 2, is: 


v of eE af vi vf 
--——+ y-Vf(4— -—+ =0, (19a) 


v ot m ov vv 





and the remainder: 


af v-VfM-v eE yf 
+=" +" .y(—)- 


ot v m v 











sf 
ot 





. (19b) 
e 
18 See reference 15, pp. 384 ff. 
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If our distribution is homogeneous (neglecting diffusion) 
Vf=0 and V-f“=0. Equation (19a) becomes after 
factoring out the v-(v is arbitrary) : 


10f 1 af € £(™ 
~ ——-++-- ——_ —E= ——. (20) 
v Of v ov m l 





If we average Eq. (20) over all directions of v we get 
no new result. However, averaging Eq. (19b) over all 
directions of v yields something different. We have: 


v VEM)-y v 
(——* =-V-f=0, 
v ay 3 


vf 1 @ 
o()L-se 
v mw 30" Ov 


Now in the steady state df /d/=0, if the electron 
density is to remain constant. Thus, Eq. (19b), when 
averaged, finally becomes: 


(0) 

a “(—e. f »)- =. 

3v? dv bt 
Now suppose E= Ey.c.+Es.c.= Eo+E,e'. The field 
will induce the corresponding frequencies in f“ so that: 


f= FY4 FP, Yeiet 


F, is complex because we might have a phase change. 
If we substitute the real parts of f and E into Eq. (21) 
we obtain: 


and 





(21) 





where F,‘°=u+iw. 


1 dfv’e 
—— — Real part (E,+E,e‘**‘). 
3v? dv 


(0) 





Real part (Fo Five) |= (22) 





6t 


e 


We must use the real parts here since { is real. Also 
Eq. (20) yields after equating static and time depend- 
ent parts: 











L af eEy 
F,“) = —- —_ — (23a) 
v ov m 
iwl+1 L af eE, 
( )rw=--———. (23b) 
v v ov m 
Instead of Eq. (23b) we have: 
vl (df /dv)(eE,/m) wl?(df / dv) (eE,/m) 
u=-— => . 
v?+ wl? v?+w?l? 


By virtue of these results Eq. (22), when averaged 
over ¢ (f is independent of time) becomes: 


1 9 (vie? af LEY? 6f 
lo al ieeoll” ol 
1+ (w/*/v?) bt 


30? Ov 








(24) 





m*> dv 
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If E,;«KE whether w//v is large or small compared with 
unity, we can neglect the term containing (£;)*. This 
assumption is certainly true for a weak signal passing 
through a discharge maintained by a high d.c. poten- 
tial. Also since the discharge is maintained by a high 
d.c. field, with collisions mostly inelastic, one can as- 
sume with some assurance that the rate of change of 
f{ with encounters is independent of E,..., the alter- 
nating field. Thus from Eq. (24), we see that f® is 
independent of Ea... 

Now the current density J= nev where n is the elec- 
tron density and V is the average electron velocity. We 
find ¥ as follows: 


1 r) ” 2e 
v=- J f f vfv"dv sin6déd ¢, 
See “oe 0 


where @, ¢ are the polar angle and azimuth in velocity 
space. However, since f= f+ (v/v)f™, we obtain 


v= ax/3n J £edv, 
0 
Thus: 


J = (4re)/3 J (Fo + F, Me") v8do, 
0 


and using Eqs. (23a) and (23b) we have: 


dre p* le afr Exe’ 
=—-— — Ey pe 
1+70l/v 





3 So vm av L 


We see therefore that the menaing current density 
Jac. =OcEac. where: 


4re p® le afr 1 
3 -" onan 


3 Jo vom dv Lit (iwl/o 





[ae (25) 


Thus Ohm’s law holds, since o, is independent of Ea.c.. 

Since we have assumed that the discharge is main- 
tained solely by the high d.c. field such that the ratio 
of inelastic to elastic collisions is high it should prove 
interesting to substitute an expression for f which is 
suited to such.a case, into Eq. (25). An expression for 
{© which is Maxwellian, holds'® under the condition 
that the probability of an inelastic encounter is inversely 
proportional to the energy. This Maxwellian formula: 
{=A exp(—/v’) is a special case of a more compli- 
cated expression, and it is used here for the sake of 
simplicity. Here, h°=[$m?/1;,/(eE)? ], where li, is the 
inelastic mean free path, and A is a normalization con- 
stant found to be n(h/x)!. If the above Maxwellian 
expression is substituted into Eq. (43), we obtain an 
expression exactly the same as was found by H. Mar- 
genau,” except for the fact that the temperature of the 
electron gas is different. This expression is: 


4 e'ln 


o.= 3 3 Gama ent) Ranta], (26) 
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where k is the Boltzmann constant, T is the absolute 
temperature, x«;=[m(wl)?/2kT ] and where: 


© xPe-# mv? 
dx x«x=—. 


K,(x1)= 
0 %+x 2kT 





These K’s have been shown” to be expressible in 
terms of the exponential integral, E;(—x,), and the 
error function erf{[_(x:)']. The only difference between 
our case and that of Margenau is that our tempera- 
ture T= (Il;,/6)*(eEo)/k as can easily be seen since 
h=m/(2kT) for a Maxwellian distribution. 

Now since x;=[m(wl)?/2kT ], we have 


(6)! mw? 7? 
2 eEs (Ilin)® 





(27) 


With moderate d.c. fields, for high frequency waves 
traveling through a low pressure gas, x;—>. Again 
using Margenau’s result for this limiting case, we find 
for our problem 


_16 en seEo\  fllin\? 
—(— ) ( ) E, coswl 
3 mw*l\2xrm 6 
e*n 


+—E, sinwt. (28) 
mw 
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For moderate d.c. fields, and extremely high fre- 
quency waves passing through a very low pressure gas 
only the second term on the right of Eq. (28) is of 
importance. 

Thus the imaginary part of conductivity o7= 
—(e’n/mw). The effective dielectric constant, K=e« 
+(¢7/w), becomes: 


K=e—(e?*n/mw?). (29) 


But for most gases that are not composed of permanent 
dipoles and especially if the gas is of low density, e has 
the value Ko, the dielectric constant of vacuum.” Thus 
the dielectric coefficient 


K’ =1—(e?n/mw?Ko). (30) 


’ We see therefore that at high electron densities, the 
dielectric coefficient is less than unity and may be 
negative. 

The author wishes to express his appreciation to Pro- 
fessor Henry Margenau for suggesting this problem and 
for the indispensable assistance, advice and encourage- 
ment he has given him. 

The author also wishes to express his gratitude to 
Professor W. P. Allis and D. E. Kerr of M.LT. for 
their stimulating discussions of this problem. 


9 Leigh Page, Introduction to Theoretical Physics (D. Van Nos- 
trand Company, Inc., New York, 1935), p. 383. 
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It has been observed that the microwave power loss of a section of wave guide filled with iron powder 
material could be varied by subjecting the iron powder to an externally applied magnetic field. This paper 
describes some experiments performed with various kinds of iron powders. Measurements have been made 
with the magnetic field applied parallel to, and perpendicular to the plane of the magnetic component of 
the electromagnetic wave inside the guide. The possibility of ferromagnetic resonance is considered. A theo- 
retical formula for the power loss exhibited by low conductivity iron powder attenuators is developed and 


compared with experimental measurements. 





INTRODUCTION 


ARLY in 1947 a series of experiments were under- 
taken at the Westinghouse Research Laboratories 
to determine the effectiveness of iron powder materials 
as X band wave guide attenuators. While performing 
certain tests, it was observed that the insertion loss of 
such attenuators could be varied by subjecting the iron 
powder to an external magnetic field. In the hope of 
developing practical magnetically controlled attenu- 
ators and wave guide amplitude modulators, a consider- 
able amount of experimental data was collected on the 
effect of a magnetic field upon the microwave trans- 
mission properties of various iron powder slugs com- 
pletely filling a section of X band wave guide excited in 
the fundamental TE», mode. This paper will present 
and discuss some of this data. A theoretical formula for 
the attenuation of iron powder materials inside wave 
guides will be derived and compared with experimental 
measurements. 

A search of the literature revealed a communication 
to Nature magazine by J. H. E. Griffith' in which is 
reported the effect of a magnetic field applied to a 
ferromagnetic film electroplated on the inside wall of a 
microwave cavity. In a paper by Charles Kittel? inci- 
dental reference is made to an M.I.T. Radiation Labora- 
tory report describing the characteristics of iron powder 
cores at microwave frequencies. 

More recently W. A. Yager and R. M. Bozorth’ have 
published experimental data on the ferromagnetic 
resonant effect observed with a supermalloy foil in a 
microwave cavity operating at 24,050 megacycles and 
subjected to a magnetic biasing field of 4920 gauss. 
Their results confirm Kittel’s* theoretical deductions 
concerning the ferromagnetic resonance effect at micro- 
wave frequencies. 

‘Although the experimental procedures reported in 


1 J. H. E. Griffith, Nature 158, 670 (1946). 

2 Charles Kittel, “Theory of the dispersion of magnetic perme- 
ability in ferromagnetic materials at microwave frequencies,” 
Phys. Rev. 70, 281 (1946). The report in reference is Rad. Lab. 
Report No. 906 which to date has not been available to the writer. 

*W. A. Yager and R. M. Bozorth, “Ferromagnetic resonance 
at microwave frequencies,” Phys. Rev. 72, 80 (1947). 

‘Charles Kittel, “Interpretation of anomalous Larmor fre- 

uenciés in the ferromagnetic resonance experiment,” Phys. Rev. 
iL, 270 (1947). 
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this paper are in practice somewhat similar to that 
reported in the foregoing references (i.e. one observes 
the effect of a magnetic field on a ferromagnetic material 
at microwave frequencies), it is believed that the opera- 
tion of the magnetically controlled iron powder attenu- 
ators does not depend primarily upon the ferromagnetic 
resonance effect. Instead it is thought that the applica- 
tion of the magnetic field varies the incremental perme- 
ability of the iron powders, and in this manner the 
hysteresis and eddy current losses at the microwave 
frequencies may be controlled by the magnetic field. 
However in certain cases, possible resonance effects do 
appear and are described in a subsequent section of this 


paper. 
PREPARATION OF THE ATTENUATORS 


The iron powder attenuators were made by mixing 
in discreet proportions various iron powders with poly- 
styrene insulating material. Before mixing, the poly- 
styrene was first dissolved in 75 cubic centimeters of 
toluene. After mixing, the solution was spread on a large 
flat dish, inserted into an oven held at a temperature of 
100°C and dried for about 12 hours. The resulting 
brittle mixture was removed from the plate, inserted 
into a one-inch diameter mold, heated at a temperature 
of 140°C until the enclosed polystyrene mixture began 
to melt, and allowed to cool to room temperature under 
a pressure of 3000 pounds per square inch. After re- 


. moval from the mold, the iron powder cylindrical ingot 


was cut into small slabs which could be inserted into a 
section of wave guide through narrow vertical slots cut 
in the side walls of the guide. 


EXPERIMENTAL APPARATUS 


The block diagram of Fig. 1 represents the wave 
guide apparatus and associated equipment used in the 
various experiments to be described. The variable mag- 
netic field required for the experiments was conveniently 
obtained from a magnetron electromagnet. This magnet 
is wound with 6000 turns of number 18 wire. The pole 
faces are 0.75 in diameter, providing an air gap area of 
1.727 square inches. The pole faces were set one inch 
apart and the air gap flux measured as a function of the 
coil current. The residual flux remaining in the air gap 
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TABLE I. Physical characteristics of carbonyle iron powders. 











To : 
retained Weight Par- 





Tron Carbon Oxygen Nitrogen by 325 average ticle 

content content content content mesh dia den- 

Grade %) (%) (%) (%) screen microns sity 
L 99.7-99.9 0.005-0.03 0.1-0.2 0.005-0.05 up to30 20 7.85 
HP 99.7-99.8  0.01-0.3 0.1-0.3 0.01-0.1 0.5-2.0 10 7.86 
Cc 99.4-99.8  0.03-0.12 0.1-0.3 0.01-0.1 0.5-2.0 10 7.86 
E 97.9-98.3 0.65-0.80 0.45-0.60 0.6-0.7 Trace 8 7.77 
TH 98.1-98.5 0.5-0.6 0.5-0.7 0.5-0.6 Trace 5 7.79 
SF 98.0-98.3 0.5-0.6 0.7-0.8 0.5-0.6 Trace 3 «7.81 








after the coil current was reduced to zero was approxi- 
mately 35 gauss. 

The experiments were performed in standard X band 
wave guide with the magnetic field applied in one or 
the other of two directions. The “b—0b” direction will 
refer to the situation in which the magnetic field is 
applied parallel to the broad dimension of the guide. 
The term “a—a”’ direction will be used to designate 
application of the magnetic field perpendicular to the 
broad dimension of the guide. Unless otherwise specified 
all measurements were made at a free space wave- 
length of 3.2 centimeters and an incident: power level 
of 0.5 milliwatt. 


EXPERIMENTS WITH CARBONYLE POWDERS 


The first experiments consisted of measuring the 
microwave transmission characteristics of samples made 
with various carbonyle iron powders supplied by 
General Analine and Film Corporation. The physical 
characteristics of these powders as furnished by the 
manufacturer are shown in Table I. The samples were 
all 0.063 inch thick and were cut from ingots made by 
mixing 1.23 ounces of iron powder with 0.18 ounce of 
polystyrene. The attenuation loss and reflection coeffi- 
cient of these samples were measured as a function of 
the magnetic field applied in the “b—b” direction. The 
initial reflection (i.e., reflection at zero magnetizing 
field) from each of the samples was “matched out” with 
a three stub tuner situated directly in front of the 
sample. 

The curves of Fig. 2a indicate the performance of 
the various iron powders as magnetically controlled 
attenuators. The samples made with the smaller particle 
size iron powders apparently make the most effective 




















attenuators. The curves of Fig. 2b illustrate the re- 
flection characteristics of these samples. — 

Several samples were made from an insulated iron 
powder supplied by the Stackpole Company and desig- 
nated by them as G-1 powder. This powder is the same 
as the carbonyle type E powder except for the particle 
insulation applied by an unknown Stackpole process. 
In order not to damage this insulation, no toluene was 
used in making an ingot from this iron powder. Instead 
the polystyrene was added to the iron in the form of 
finely crushed crystals. The two ingredients were dry 
mixed, and a test ingot made in the manner pre- 
viously described. 

The data in Fig. 3a show the attenuation character- 
istics of three different size samples made from an 
ingot consisting of 1.23 ounces of G-1 Stackpole iron 
powder and 0.18 ounce of polystyrene crystals. It is 
interesting to note that the initial insertion loss is 
almost directly proportional to the thickness of the 
sample, but as magnetic saturation is approached this 
proportionality no longer holds. The reflection proper- 
ties of these samples are displayed in Fig. 3b. No 
explanation can be offered for the large values of 
P.S.W.R. exhibited by sample number seven. 


FERROMAGNETIC RESONANCE 


Possible ferromagnetic resonance effects were ob- 
served with an iron powder sample made from G-1 iron 
powder mixed with a thermo-setting resin. This powder 
was also obtained from the Stackpole Company and is 
designated as G-1 MIX. Since a resin binder was 
already available in this powder, it was not necessary 
to add any polystyrene in making the ingot. The 
attenuation and reflection characteristic of a test sample 
(0.063” thick) was measured at two different wave- 
lengths, and with the magnetic field applied in the 
“a—a” and the “b—b” direction. The performance of 
this sample when the magnetic field is applied in the 
“a—a” direction is shown in Fig. 4. For this direction 
of the magnetic field a pronounced peak can be ob- 
served in the attenuation characteristics. Such peaks 
are indicative of resonance absorption of microwave 
energy. 

Kittel’ has shown that for ferromagnetic resonance, 
the magnetostatic field applied normal to the magnetic 
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* Charles Kittel, “Ferromagnetic resonance absorption,” Phys. Rev. 73, 155 (1948). 


VOLUME 20, SEPTEMBER, 1949 








component of the microwave field must satisfy the 
following condition: 


wo= 7 { [H.+ (Ny- N,)M.)X [H.+ (V.- N.)M,)} , 


where wo=angular microwave frequency; H.=mag- 
netostatic field strength ; M.=component of magnetiza- 
tion vector parallel to magnetostatic field; V., NV,, N; 
=demagnetizing factors of sample in guide; and 
y= electron spin gyromagnetic ratio (y/24r= 2.80 mega- 
cycles/oersted). Because of the difficulties involved in 
calculating proper values for the demagnetizing factors, 
and the uncertain assumptions that would have to be 
made in assigning a value to M,, it is impractical to 
attempt to compute the resonant frequency condition 
for the iron powder sample. However from Fig. 4 it may 
be seen that as the microwave frequency increases the 
peaks of the attenuation curves occur at greater mag- 
netostatic field strengths. This is at least in qualitative 
agreement with theory. 

The attenuation and reflection characteristics of 
sample number 16 with the magnetic field applied in 
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Fic. 2a. Insertion loss of samples made with carbonyle powders, 


the “b—b” direction are shown in Fig. 5. No peak or 
resonance attenuation occurs for this direction of the 
magnetic field. In comparing Fig. 4 with Fig. 5 one is 
impressed with the tremendous difference in the reflec- 
tion characteristics. In designing a magnetically con- 
trolled attenuator for which a very low reflection 
coefficient is desirable, it apparently would be advan- 
tageous to apply the magnetic field in the “a—a” 
direction. 

The static permeability of sample number 16 was 
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measured by the ballistic galvanometer method® and 
found to be equal to 1.80. 


EXPERIMENTS WITH PERMALLOY POWDER 


Several samples were made from 3 micron diameter 
Permalloy powder obtained from the Inter Chemical 
Corporation of New York. These samples were cut from 
ingots made by introducing the iron powder into dis- 
solved polystyrene as previously described. The attenu- 
ation and reflection characteristic of two of these 
samples are shown plotted in Fig. 6. The magnetic 
field was applied in the “b—b” direction. Sample 
number 10 was made from an ingot consisting of 1.23 
ounces of Permalloy and 0.18 ounce of polystyrene; 
sample number 11 was cut from an ingot consisting of 
0.50 ounce of Permalloy and 0.15 ounce of polystyrene. 
Because of the comparatively steep slopes of the 
attenuation curves, these samples appear to make the 
most effective magnetically controlled attenuators. This 
probably can be attributed to the high initial perme- 
ability of Permalloy powder which saturates at rela- 
tively small magnetizing field strengths. 


THEORETICAL CONSIDERATIONS 


At microwave frequencies a low conductivity di- 
electric loaded with ferromagnetic materials will exhibit 
a complex dielectric constant and a complex perme- 
ability.” This section of the paper will consider the 
theoretical power loss experienced by a TE, wave 
propagating through a section of wave guide filled with 
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Fic. 2b. Reflection characteristics of samples made 
with carbonyle powders. 


6 Page and Adams, Principles of Electricity (D. Van Nostrand 
and Company, Inc., New York, 1931). 

7J. B. Birks, “Measurement of the permeability of magnetic 
— at centimeter wave-lengths,” Nature 158, 671 
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a loaded dielectric material. In view of the assumptions 
that must be made, acceptable correlation is obtained 
between the theoretical and the measured loss values. 

Maxwell’s equations for a wave traveling in a medium 
exhibiting complex electrical constants become in m.k.s. 
units 


VXE=—(u;+ine)dH/dt 
V: (uitins)H =0 
V- (e:+ie,)H=0 


where the subscripts 1 and 2 denote, respectively, the 
real and imaginary parts of the permeability and di- 
electric constant. Because of the extremely low con- 
ductivity of the medium under consideration, the 
conduction current term in Maxwell’s equations may 
be neglected. : 

A solution of the wave equation, and application of 
the boundary conditions for a wave guide excited in the 
TE: fundamental mode, yields the following definitions 
for the electric and magnetic field components of a wave 
traveling along the Z direction 


Ey= Ke~® sin(2x/b) cos(wt— BZ) 


(ui—ipe)e*% =x 
H,=-—K 





S 
(u41°+ wo”) bw b 
X |B cos(wt— BZ)+ a sin(wi— BZ)} 


—Kr(ui-—ipe)e 2 x 
ony’ sin(wt— BZ) 





Hz= 


(417+ po”) bw 
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Fic. 3a. Insertion loss of three different size samples made 
with stackpole G-1 iron powder. 
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where 6 is the wave guide dimension along the x direc- 
tion (broad dimension). The propagation constant 
y=a+iB evaluates to 


= (9?/b) — w?(urt+ipe) (€1+tee) ]}. 


Separation of the real and imaginary parts leads to the 
following expressions for the attenuation constant and 
phase constant. 


a=[({2/b}?—w*{ wie1— ures} )* 

+ (w*{ u2€1+p1€2} )? }* cos(8/2), 
B= C( 1/b}*—w{ w1€1— wee} )? 
+ (w?{ oer t ures} )* }* sin(6/2), 
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Fic. 3b. Reflection characteristics of samples made 
with stackpole G-1 iron powder. 


The total instantaneous power flowing through the 


guide at any point may be obtained by integrating the 
Poynting vector over the cross-sectional area of the 


guide. Thus 
z= y=a dl 
Pi= f f (EXH)dxdy. 
0 0 


Substituting Ey and Hx for FE and H, and averaging by 
integrating over a complete time period, the average 
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power flow across a section of guide evaluates to 


a K*e-?42 4 Ba: b 


4 (17+ 2”) 





“q” is the wave guide dimension in the y direction. 
Establishing the ratio of the average power at point Z 
to the average power at point (Z+1) meters 


P2/(P 24» ]=e*. 
Converting this power ratio to decibels 
db= (10/2.3) log.e**. 


Thus the decibel power loss per meter length of guide 
may finally be written as 


db= 20/2.3[ ({ r/b}?—w?{ wre: — ree} )? 
+ (w*{ w2ert wi€2} )* }' cos(8/2), (3) 


where @ is given by Eq. (2). 



































wr ATTENUATION at 3.075 CMs. 
” ™*\ crema ath anon 

/ i iat e 
2 ‘Zs l “ \ had 
8) 7 \ ; 
z\|7 \ \ 5 
> | PSWR. AT A=3.075 - : 
e eS all \ 1.6 & 
3 2 PW.S.R. ATL =3.325 a * 8 

12 wre om, 12 

, 8 

2 3 a 














U 4 : 
MAGNET AIR GAP FIELD INTENSITY~KILOGAUSS 


Fic. 4. Attenuation and reflection characteristics of sample No. 16 
with magnetic field in a—a direction. 


It is now necessary to determine the effective perme- 
ability and dielectric constants of the iron powder 
samples. Kittel? has shown that the effective perme- 
ability of iron at X band frequencies is complex and is 
given approximately by 


p= 9e-**” in c.g.s. unit. 


Separating real and imaginary parts and converting to 
m.k.s. units. 


pi= 1.96 10-6 henry/meter 
we= 11.1 10-6 henry/meter. 


A typical handbook value for the dielectric constant of 
polystyrene at X band frequencies is 


€:= 21.7X10-" farad/ meter, 
€2= 1.06 10-" farad/meter. 
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From theoretical considerations based on scattering 
of an electromagnetic wave from spherical particles 
imbedded in a non-conducting medium, Lewin® derives 
the following formulas for the bulk constants of a loaded 
dielectric. 


ém= eal 1+ (3f/1—f) ] 





3f 
tama 1+ ) 
(up— 2ua)/(up— Ma) —f. 


where ¢€,=dielectric constant of loaded dielectric; 
€a= dielectric constant of dielectric only; u»=perme- 
ability of loaded dielectric; ya=permeability of di- 
electric only; 4,»=permeability of iron particles at 
frequency under consideration; and f=ratio of volume 
of particles to total volume of mixture=0.504 for 
samples shown in Fig. 2a. Substituting the given values 
for the parameters in Lewin’s formulas, and treating 
the real and imaginary parts separately, there are 
obtained the following values for the bulk constants of 
the loaded dielectric. 


€:=91.8X10-” farad/meter, 
€2= 4.48 X 10-4 farad/meter, 
i= 1.59X10~* henry/meter, 
po= 3.54 10-* henry /meter. 


Substituting these constants into Eq. (3) and neglect- 
ing the possibility of ferromagnetic resonance, the loss 
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Fic. 5. Attenuation and reflection characteristics of sample No. 16 
with magnetic field in )—6 direction. 


8 L. Lewin, “The electrical constants of a material loaded with 
spherical particles,” J. Inst. Elec. Eng. 94, 65 (1947). 
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in the loaded dielectric becomes 


Power loss= 7.91 decibels per millimeter 
(at zero magnetizing field). 


Assuming that for magnetic saturation the real and 
imaginary parts of the effective permeability decrease 
to that corresponding to free space, and that the di- 
electric constant of the loaded dielectric is not a function 
of the magnetic field, the power loss may be computed 
under conditions of magnetic saturation and compared 
to some of the experimental curves shown in Fig. 2a. 
The loss at magnetic saturation evaluates to 5.89 
decibels per millimeter. The sample thickness equals 
1.6 millimeters. Consequently the computed loss with- 


out a magnetostatic field is 12.65 db and the com- , 


puted loss at magnetic saturation is 9.32 db. These 
loss values may be compared to the measured values 
for sample number 3 in Fig. 2a which has an initial in- 
sertion loss of 11.8 db and a magnetic saturation loss 
of 8.4 db. 

It is interesting to note that in most cases the theory 
predicts rather accurately the difference between the 
initial and minimum insertion loss values for the 
samples shown in Fig. 2a. The most striking deviation 
from the value predicted by the theory occurs with 
sample No. 2. This may be caused by ferromagnetic 
resonance which might sustain the permeability at 
large magnetostatic field strengths. 

The inability of the theory to determine the initial 
insertion loss for all of the samples, stems from the fact 
that the actual radii of the iron particles does not enter 
into Lewin’s formulas—only the volume concentra- 
tion.” If sufficient data were available it might be 
possible to compute an empirical constant which would 
take into consideration the variation in particle size. 
However such factors as sample homogeneity, particle 
geometry, and particle isolation would also have to be 
considered. 


®The measurements made with sample number 5 cannot be 
compared with the theory because for this sample the iron par- 
ticle size is several times greater than the skin depth at the 
frequency under consideration, and therefore only a small portion 
of the iron is effective in controlling the attenuation. 

10 Another factor not considered in the theory is the power loss 
associated with the circulating currents in the matching section. 
However measurements indicate that the maximum matching loss 
for any of the samples considered in Fig. 2a is only about 0.4 db. 
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CONCLUSION 


The attenuators described in this paper require very 
large magnetizing forces of the order of five thousand 
ampere turns, and for this reason are not too well suited 
for practical applications. However, some work has 
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Fic. 6. Attenuation and reflection characteristics of samples 
made with Permalloy powder. 


been done toward reducing the strength of the required 
magnetic field, and several attenuators and amplitude 
modulators have been developed which require less 
than 1000 ampere turns of magnetizing force. It is 
planned to report on these attenuators in the near 
future. 
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Aniline and Film Company in supplying many of the 
carbonyle powders used in the experiments. Acknowl- 
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experimental data presented here. Thanks are also due 
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N a recent article, Mahlman!' encounters some difficulties in 
explaining his results on the basis of the assumption that an 
oxide-coated cathode can be treated as an excess-semiconductor. 
Some time ago we proposed? a conduction mechanism, of the oxide- 
coated cathode, which, we think, offers an explanation of the fact 
that according to his measurements 


¢thermion. + 4E > ¢photoelectr., 


4E being the activation energy of the conductivity. 

In our theory particular stress is laid on the fact that an oxide 
coating is an electron-emitting semiconducting solid with a 
porous character. Because of its electron emission at high tempera- 
tures (in particular>800°K), each grain of oxide in the porous 
coating is surrounded by a cloud of electrons. This means that all 
pores can be considered as filled with an electron gas the density 
of which can easily be calculated. On the assumption that the 
mean free path of these electrons is of the order of magnitude of 
the average pore diameter (a few microns), a conductivity of an 
oxide coating can be calculated which is of the same order of 
magnitude as that actually found in experiments. 

It will be understood that the temperature dependence of the 
density of the electron gas and consequently of the conductivity 
will be the same as that of the thermionic emission. A calculation 
showed this to be correct in a first approximation. 

Extensive experiments showed that the loge against1/7-curves 
always had the form reproduced in Fig. 1. The bent is situated at 
a temperature of about 750°K. An obvious explanation of this 
form of the curves is that there are two conduction mechanisms 
in the oxide coating, connected in parallel, one with a low activa- 
tion energy preponderating under 750°K, the other with a higher 
activation energy preponderating above 750°K. Because of this 
connection in parallel it is not allowed to calculate activation 
energies directly from the slope of one of the parts of the curve in 
Fig. 1. The true temperature dependence of each conduction 
mechanism can only be found if one bears in mind that the con- 
ductivity experimentally found is the sum of two separate con- 
ductivities. A typical example of the activation energies computed 
on this basis is given in Table I. 

The rather good agreement between the values in columns 2 and 
4 (compare with those of column 3) is a confirmation of our theory. 
We suppose the activation energies preponderating under 750°K 
(column 1) to be those of the conducting grains themselves. 

Our explanation of the results of Mahlman is: 
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Fic, 1, Loge, 1/T—curve for an oxide coating, 
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TaBLe I. Activation energies and work functions of one oxide-coated 
cathode in various degrees of activation. 











1 2 3 4 
Appararent 
— activation 
Activation Activation energy calcu- 
energy — energy lated directly Thermionic 
preponderating preponderating from the slope work 
under 750°K above 750°K above 750°K function 
Volts Volts Volts Volts 
we 0.22 0.98 0.84 1.07 
Es 0.14 1.14 0.60 —_ 
ac 0.12 1.24 0.54 1.10 
os 0.11 0.96 0.56 0.90 
28 0.10 0.96 0.54 1.02 
= 0.09 0.94 0.48 0.87 








(1) That his activation energy of the conductivity is an ap- 
parent activation energy. The true activation energy will be some 
tenths of a volt higher as can be seen in Table I. It cannot be esti- 
mated from his measurements as he did not extend them below 
750°K. 

(2) That according to our opinion the true value of 4£ is not 
1.2 v but 0.2-0.4 v, as could be estimated from the degree of ac- 
tivation of his cathode and our experimental results. If no allow- 
ance is made for the Franck-Condon principle 


¥photoelectr. = ¢thermion. + sE. 


As Mahlman found ¥thermion. = 1.5 V, ¢photoelectr. would be 
1.7 to 1.9 v which is in striking agreement with the 1.8 v he ob- 
served. A paper with full details of our investigations is in prepara- 
tion.® 

1G. W. Mahiman, J. App. Phys. 20, 197 (1949). 

2H. J. Vink, Het geleidingsmechanisme in de oxyde-kathode, thesis, Leiden 


University, 1948. 
3 To be published in Philips Res. Rep. 





Note on Some Plastic Flow Effects in Steel 


Orrin H. CLARK AND ERNEST M. LyMAN* 


Socony-Vacuum Laboratories (A Division of Socony-Vacuum Oil 
Company, Inc.), Research and Development Department, 
Paulsboro, New Jersey 


June 13, 1949 


N application of interference fringe measurements to the 
study of plastic flow in steel has given some interesting 
results. In the course of studies of the lubrication of bearing metals, 
it became important to know the approximate magnitude of the 
depression existing at the center of a contact circle between two 
dissimilar steel pins pressed together under load sufficient to 
induce plastic flow in the softer one. The solution obtained by 
interference fringe measurements illustrates a method of studying 
deflections occurring in simple cases of plastic flow and is therefore 
described briefly. 

A photograph of the interference fringe contour map of a resid- 
ual pit is shown in Fig. 1. The pit was made in a soft (Rockwell 
B=90) steel pin of }-in. diameter. Another }-in. steel pin with 
axis at right angles was pressed into the soft pin under 25 lb. 
The loading pin was hard (Rockwell C=64) and suffered no 
permanent deformation. The picture was taken with sodium light 
by laying an optical flat on the pit. A standard metallurgical 
microscope with vertical illuminator was used. 

This contour map of course does not represent conditions under 
load, because of the elastic recovery. It is convenient in correcting 
for elastic recovery to consider it as taking place in two stages. 
The first stage is the shrinkage of the contact circle to the size of 
the residual pit. That is, under load, the contact circle exceeds 
the measured pit size by an amount consisting of an annular ring 
of elastic contact outside this region of permanent deformation. 
When unloading, the contact area first “peels off” the annular 
ring of elastic contact, both pins incurring vertical rise at the 
center in this stage. Further vertical rise takes place during the 
second stage of removal of the rest of the load. 


JOURNAL OF APPLIED PHYSICS 








FIG 


Th 
betwe 
howe 
since 
of thi 
sagitt 
defor 
The 1 
fact, 
pin’s 
carrie 
of the 

Th 
sagitt 

Co 
from 
fringe 
of for 
meas' 
amet 
being 
Henc 
curva 

It: 
diam: 
was ( 
ment: 
conta 
this r 
or 0. 
4.45> 
to be 
cm, ¢ 
the s 
the h 
displ. 
havin 
in the 
elasti 
cover 
Unde 
3.40> 

Th 
differ 
An e 
feren 
case 
the c 


Fic 


VOL 





i- 


he 
ng 
ls, 
he 
vO 
to 
by 
ng 
re 


ell 
ith 
Ib. 


ght 
cal 


der 
ing 
res. 
> of 
eds 
ing 
on. 
ilar 
the 
the 


ics 











Fic, 1. Fringe pattern in pit left in soft pin after loading with hard pin. 


The interference fringes give directly the difference in depth 
between the center and each edge, after unloading. The edges, 
however, suffer no vertical recovery during the second stage, 
since they are under no pressure. During the first stage the rise 
of the transverse edges can be calculated from the conventional 
sagittal formula, provided the curvature of the loading pin is not 
deformed appreciably, so that it may be used to give the arc. 
The fringes show that the edge depth in the residual pit does, in 
fact, correspond to the sagittal value given by using the loading 
pin’s radius. Since the elastic contact annulus is very small and 
carries very little load, it is a good assumption that the curvature 
of the loading pin at the transverse edges is unchanged throughout. 

The figures showing the correspondence of edge depth to the 
sagittal calculation are as follows. 

Counting fringes, it is seen that there is a rise of six fringes 
from the transverse edges to the center, and a further rise of eight 
fringes from the center to the axial edges, making a total depth 
of fourteen fringes. Calculation of the depth of the arc from the 
measured radius of the loadung pin (0.310 cm) and the pit di- 
ameter (0.0315 cm) gives the result 4.0010 cm. Each fringe 
being a depth of 2.95 10-5 cm, this corresponds to 13.6 fringes. 
Hence, the curvature of the pit after unpeeling is the same as the 
curvature of the loading pin. 

It remains to calculate the edge depth under load. For this the 
diameter of the contact area under load is needed. This diameter 
was obtained by using Prussian blue as an indicator. Measure- 
ments on fresh pins gave a value of 1.051 for the ratio of total 
contact diameter to the diameter of the residual pit. Applying 
this ratio, the contact diameter under load was 1.051 0.0315 cm, 
or 0.0332 cm. The edge depth was then, by saggital formula, 
4.45X 10-4 cm. The center depth, after unloading, has been seen 
to be eight fringes, or 2.36 10~* cm. The difference, 2.09 10-4 
cm, consists of two parts, namely, the central elastic recovery of 
the soft pin when unloaded plus the central elastic deflection of 
the hard pin under load. These two parts are the vertical elastic 
displacements of two similarly shaped pins under the same load, 
having the same elastic properties, differing somewhat, however, 
in the stress distributions away from the contact surface. If these 
elastic displacements are assumed equal, the desired elastic re- 
covery of the center of the soft pin is obtained as 1.05 10~ cm. 
Under load the depth of the central depression was therefore 
3.40 10~4 cm. 

The necessity of assuming equal elastic deflections for pins of 
different hardness does not arise in cases where like pins are used. 
An example of this is illustrated in Fig. 2, which shows the inter- 
ference fringe pattern obtained when both pins were soft. In this 
case there is a rise of eleven fringes from the transverse edge to 
the center and a further rise of four fringes to the axial edge. It is 
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Fic. 2. Fringe pattern in pit left in soft pin after loading with soft pin. 


VOLUME 20, SEPTEMBER, 1949 


obvious that the seven-fringe asymmetry was equally caused by 
elastic recovery of each pin, and that therefore the recovery in 
the second stage of unloading was 3} fringes, or 1.0310-* cm. 
For this combination the pit diameter was 0.0328 cm, and the 
factor given by Prussian blue measurements for increase due to 
elastic contact was 1.057. Applying the sagittal formula gives an 
additional elastic recovery for each pin of 0.28X10~ cm, with a 
resulting total elastic recovery of 1.3110~* cm. 

In the case of two hard pins, all deformation is elastic and the 
total elastic deflections are obtained by again applying the sagittal 
formula. The pit diameter was calculated to be 0.0266 cm from 
the standard formula of elasticity theory, and three checks with 
Prussian blue gave measurements of 0.0269, 0.0270 and 0.0270 cm. 
Using the calculated value, the central elastic deflection turns out 
to be 1.43 10~ cm for each pin. 

Comparing the combined elastic deflections of both pins in the 
three cases, hard-soft, soft-soft, and hard-hard, we get the values 
2.10 10~* cm, 2.62 10~* cm, and 2.86 10~ cm, respectively. 
Since the hard-soft case would be expected to be intermediate, 
it appears that the measuring error in the value obtained for 
elastic deflection in a single pin is of the order of 0.26 10~* cm. 
This is about the amount of one fringe and is ascribed to uncer- 
tainty in matching fringes to the pit edges. While the Prussian 
blue measurements have little precision, their total contribution 
is too small to make this a major source of error. 

* University of Illinois, Urbana, Illinois. 

1 It is interesting to observe that virtually the identical problem arose in the 
recent work of Holm, Holm, and Shobert, J. App. Phys. 20, 319-327 (1949). 
In a study of the theory of hardness, they needed to know the deflections in 
the case of a ball pressed into a plane. They arrived at a calculated result 


by assuming that the flattening of the ball and of the plate could each be 
considered separately and then combined. 





Interpretation of X-Ray Patterns of 
Cold-Worked Metal* 


B. L. AVERBACH AND B. E. WARREN 
Massachusetts Institute of Technology, Cambridge, Massachusetts 
(July 7, 1949) 


ITH modern experimental technique it is possible to 
measure the shape of a powder pattern line with sufficient 
accuracy to justify an interpretation based on the shape of the line 
rather than using just the line breadth. By making full use of the 
shape of the line it is possible to get a great deal of information 
about the nature of cold work without making any @ priori as- 
sumptions. 
Starting with the same formulation used by Stokes and Wilson! 
we represent the position of cell mmom; by 


Rmymym3= M10, + M2d2+ mM303+ dmimem3 
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Fic. 1. Relative displacements vs. crystal distance from Fourier coefficients 
for powder lines of cold-worked a-brass. 
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where 6 is an arbitrary displacement. The intensity from one crys- 
tal is then 


I=PDVVVZI exp| 774(s—s.)-(Ra— Rw). 
mi m,;’ 

Let s—so=A(hbi+h2b2+hsb3) where bjb2b3 are the reciprocal 
vectors and /,/2h3 are continuous variables. With cubic crystals it 
is always possible to choose a new set of orthogonal axes so as to 
make any reflection of the type (00/). The intensity distribution in 
the (00%) line of a powder pattern is then given by the cosine 
Fourier series 


+2 


Pa=K = Am cos2amh; (1) 


m 
-—2 


where the variable /; is given by 
sin@ 
hy = 2a 
and the coefficients A, are given by 


A m= (COS2 mlz m)aymi =m1’ (2) 
m2=mz2’. 


If the 2x/z,, values are small, or if the probability of the z,, values 
is given by an error curve (2) becomes 


—2n7P(Zn,7)mi a 


m2 =m’ 


An=exp| 


(3) 


Using crystal monochromated radiation and a Geiger-counter 
spectrometer, the intensity distribution in the line (00/) is measured 
for a cold-worked sample and for an annealed sample. The Fourier 
coefficients A,,, relating to cold-work broadening only, are then 
readily obtained from the two curves by the method of Stokes.? 
From each coefficient A m, a value of Z»,? is computed using Eq. (3). 
Consider two unit cells mymym; and m,'m,'m;' with m,=my’, 
m2= m2’ and |m;—m;'|=m. The distance between these cells is 
ma; and the component along a; of the relative displacement of the 
two cells is zma@3. The quantity (Z,,?)m: =m)’ is an average over all 

m2=my2’ 
such pairs of cells. If strains in the material are uniform over large 
regions, a plot of (Z,*)4 against m should give a straight line whose 
slope is the root-mean-square strain. Figure 1 shows results for 
70-30 a-brass filings. The initial slopes give the true root-mean- 
square strains in the corresponding directions. It is evident that 
they are very closely in an inverse ratio to the corresponding values 
of Young’s modulus, indicating that the root-mean-square stress is 
approximately isotropic.* The average initial slope is very close to 
the macroscopic yield strain for 70-30 brass. The departure of the 


curves from linearity offers a promising method for studying the 
non-uniformity of strains over regions of appreciable size and 
demonstrating the highly localized nature of the strains in cold- 
worked metal. 

A treatment of powder pattern lines in terms of the Fourier 
coefficients A,» offers a very general and powerful method for 
handling cold-work broadening, particle-size broadening, and 
many other effects. A plot of A» against m will distinguish between 
cold-work broadening and particle-size broadening. For a sub- 
stance in which there is particle-size broadening only, the initial 
slope of the A», vs. m curve gives an average particle size and the 
second derivative of the curve gives directly a size distribution. 
The integral breadth of a line 8 is given by 


1 
a 


z2Am 


and hence a measurement of line breadth gives only a small part 
of the information contained in the full set of Fourier coefficients. 


* Experimental work supported by the ONR under Contract NS 
ori-07832. 

1A. R. Stokes and A. J. C. Wilson, Proc. Camb. Phil. Soc. 38, 313 (1942); 
Proc. Phys. Soc. 56, 174 (1944). 

2A. R. Stokes, Proc. Phys. Soc. 61, 382 (1948). 

7C. S. Smith and E. E. Stickley, Phys. Rev. 64, 191 (1943). 

4 Stokes, Pascoe, and Lipson, Nature 151, 137 (1943). 





Zirconium Carbide as a Thermionic Emitter 


R. E. Happap, D. L. Go_LDWATER, AND F. H. MORGAN 


Cathode Research Group, Bartol Research Foundation of the 
Franklin Institute, Swarthmore, Pennsylvania 


June 6, 1949 


S part of a broad program of investigation, the thermionic 
emitting properties of zirconium carbide are being studied. 
The zirconium carbide is cataphoretically coated onto tungsten 
wire and the cathode is then mounted in a guard-ringed diode. 
Values of A and 6b are obtained for the Richardson equation 
I,=AT*e~/T, (b6=11,605¢). First results yield values of ¢=2.1(6) v, 


A=0.2—0.5 amp./deg. K?/cm?. 


Two interesting values of the space-charge limited emission 
available are: at 1500° C,, at least 0.6 amp./cm? d.c., at 1700°C,, 
6 amp./cm? pulsed. 

In a subsidiary experiment, the spectral emissivity of zirconium 
carbide has been obtained for temperatures up to 1800°C, and at 
a wave-length of 0.6554. The value obtained is 0.96+0.04. 
Brightness temperatures given above can hence be regarded as 
true values. 
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